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Zusammenfassung 



Das Standardmodell der Teilchenphysik, welches in den spaten sechziger Jahren in 
seiner heutigen Form eingefiihrt wurde, hat sich iiber die letzten 40 Jahre hinweg 
in zahllosen Experimenten als eine aufierordentlich genaue Beschreibung der Phy- 
sik auf kleinsten Skalen erwiesen. Alle von dem Modell vorhergesagten Fermionen 
und Vektorbosonen wurden mittlerweile entdeckt, ihre Eigenschaften experimen- 
tell vermessen, und die Vorhersagen des Modells auf Basis dieser Daten haben sich 
auf Schleifen- und sogar Mehrschleifenniveau als gultig erwiesen. 

Der letzte noch nicht entdeckte Baustein des Modells ist das Higgs-Boson. Da 
von den bisher durchgefiihrten Experimenten im wesentlichen der Skalenbereich 
unterhalb von 100 GeV abgedeckt wurde, ist dies jedoch noch kein Problem und 
immer noch mit der Theorie konsistent. Der LHC (der in der nahen Zukunft anlau- 
fen sollte) wird den experimentell zugangliche Skalenbereich voraussichtlich um 
eine Grofienordnung auf mehrere TeV erweitern. Sollte das Higgs auch bei diesen 
Energien nicht gefunden werden, so ware das Standardmodell falsifiziert und aus 
den Grundlagen der Quantenmechanik wiirde f olgen, dafi entweder bei etwa 1 TeV 
die Storungstheorie zusammenbrechen oder neue Physik in Erscheinung treten 
mufite, welche die Perturbativitat der Theorie gewahrleistet — anderenfalls wiirde 
bei hohen Energien die Unitaritat des Zeitentwicklungsoperators verlorengehen 
und die statistische Interpretation der Quantenmechanik zusammenbrechen. 

Obwohl viele Modelle neuer Physik das Konzept der Symmetriebrechung durch 
fundamentale Skalarfelder vom Standardmodell iibernehmen und lediglich die De- 
tails des entsprechenden Sektors modifizieren, gibt es auch eine grofie Gruppe von 
Szenarien, welche keine derartigen Felder enthalten. Unter diesen Higgslosen Mo- 
dellen (und moglicherweise auch durch die beriihmte AdS/CFT-Korrespondenz 
verkniipft) sind die wohl bekanntesten Beispiele Technicolor-Modelle sowie extra- 
dimensionale Higgslose Szenarien. Beiden Klassen von Modellen gemeinsam ist 
das Auftreten neuer Resonanzen im Spektrum oberhalb von etwa 100 GeV, deren 
Austausch perturbative Unitaritat bei hohen Skalen gewahrleistet. Diese neuen Re- 
sonanzen stellen jedoch ebenfalls eine ernstzunehmende Gefahr fur solche Szena- 
rien dar, da ihr Austausch zusatzliche Beitrage zu den bei LEP / LEP-II sehr ge- 
nau vermessenen Prazisionsobservablen liefert. Falls keine speziellen Vorkehrun- 
gen getroffen werden, haben diese Beitrage die Tendenz, derartige Modelle bereits 
auf der Basis existierender Daten auszuschliefien. 

In den letzten Jahren wurden extradimensionale Modelle vorgeschlagen, in wel- 
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chen die Fermionen des Standardmodells in der Extradimension delokalisiert sind, 
wodurch die Prazisionsobservablen durch Justage der Kopplungen an die neu- 
en Resonanzen korrekt reproduziert werden konnen. Derartige Modelle sind ein 
gangbarer Weg, die elektroschwache Symmetrie zu brechen und die perturbative 
Unitaritat an der TeV-Skala aufrecht zu erhalten, ohne ein fundamentales Higgs- 
Feld zu postulieren. Allerdings sind extradimensionale Modelle (von Trivialfallen 
abgesehen) nicht renormierbar und nur unterhalb einer Cutoff-Skala gultig, und 
die meisten neuen Resonanzen liegen oberhalb dieses Cutoffs. Eine „ehrliche" Er- 
weiterung des Standardmodells sollte lediglich die Struktur unterhalb dieses Cu- 
toffs enthalten und den extradimensionalen Mechanismus zur Symmetriebrechung 
und Verzogerung der Unitaritatsverletzung implementieren, ohne Annahmen iiber 
die Physik jenseits der Cutoff-Skala zu machen. 

Das „Three-Site Higgsless Model" ist eine minimale Implementation dieser Idee. 
Obwohl dieses Modell durch extradimensionale Higgslose Modelle motiviert wer- 
den kann, enthalt es lediglich eine Generation zusatzlicher Resonanzen, welche 
vollstandig unterhalb des Cutoffs liegt und die Verletzung der Unitaritat auf 2 - 
3TeV hinausschiebt. Der nicht mit dem Standardmodell ubereinstimmende Teil 
des Spektrums besteht aus einem kompletten Satz von Partnern fur alle Standard- 
modellteilchen mit Ausnahme des Gluons und des Photons. Eine Analyse der be- 
stehenden experimentellen Einschrankungen zeigt, dafi das Modell die Prazisi- 
onsobservablen korrekt reproduzieren kann, falls die Kopplungen zwischen den 
schweren Eichbosonpartnern und den Fermionen des Standardmodells sehr klein 
(etwa 1% des Isospinkopplung) und die Fermionpartner mit Massen > 1.8 TeV 
verhaltnismafiig schwer sind. 

In dieser Doktorarbeit wurde die LHC-Phanomenologie dieses Szenarios unter- 
sucht. Zu diesem Zwecke wurden die Kopplungen und Breiten aller neuen Teilchen 
berechnet und das Modell in den Monte-Carlo-Generator WHIZARD / O'Mega 
implementiert. Diese Implementation wurde verwendet, um die Produktion der 
Fermion- und Eichbosonpartner auf Partonniveau in verschiedenen Kanalen zu si- 
mulieren, welche sich fur die Entdeckung am LHC eignen konnten. In dieser Arbeit 
werden die Ergebnisse zusammen mit ein Einfiihrung in das Modell sowie einer 
Diskussion der Modelleigenschaften prasentiert. 

Obwohl ihre fermiophobe Natur die Entdeckung der schweren Eichbosonen 
an Teilchenbeschleunigern grundsatzlich erschwert, zeigt sich, dafi der LHC die 
entsprechenden Resonanzen finden kann sowie sogar einige Ruckschlusse auf die 
Starke der fermiophoben Kopplungen (was ein wesentlicher Test der Konsistenz 
eines solchen Szenarios ware) zulassen sollte. Bei der Berechnung der Breite der 
schweren Fermionen stellt sich heraus, dafi zu der grofien Masse auch relative Brei- 
ten von 10% und mehr kommen, so dafi diese Teilchen sich eher schlecht fur eine 
direkte Entdeckung am LHC eignen. Trotzdem zeigen die Simulationen dafi, hin- 
reichend viel Zeit, Geduld sowie ein gutes Verstandnis von Detektor und Hinter- 
grund vorausgesetzt, eine direkte Entdeckung zumindest in einem Teil des Para- 
meterraums moglich ist. 



Abstract 



The Standard Model of particle physics, conceived in the late 1960s, has been con- 
firmed as an extremely accurate description of microscopic physics in a multitude 
of experiments conducted over the last 40 years. Over time, all fermions and vec- 
tor bosons predicted by the model have been discovered, their properties have 
been measured, and the predictions of the model based upon these properties have 
shown to be accurate to the one loop and even multiloop level. 

The last piece of the Standard Model which hasn't been discovered yet is the 
Higgs boson. As experiments have essentially only been probing scales below 
100 GeV, this is not necessarily a problem and still consistent with the theory. How- 
ever, the LHC (which should be commencing operation anytime soon) will be hope- 
fully extending the energy scale accessible to experiments by a an order of magni- 
tude to several TeVs. If the Higgs is not discovered in this energy range, then the 
Standard Model would be falsified and the very fundamentals of quantum mechan- 
ics would tell us that either perturbation theory must break down at about 1 TeV 
or that new physics must enter the stage in order for the theory to remain pertur- 
bative — otherwise, unitarity would be lost at high energies and the probabilistic 
interpretation of quantum mechanics would break down. 

While many models of new physics retain the concept of fundamental scalars 
in the spectrum being responsible for the symmetry breaking and just modify the 
details of the implementation, there is also a large group of scenarios which do not 
contain any such fields at all. Among these Higgsless models (and probably also 
connected by the celebrated AdS/CFT correspondence), the arguably most promi- 
nent examples are technicolor and extra dimensional Higgsless models. In both 
classes of models, new resonances appear in the spectrum above 100 GeV, the ex- 
change of which retains perturbative unitarity at high scales. However, at the same 
time, the presence of these new resonances also turns out to be a severe threat to 
such scenarios, their exchange leading to additional contributions to the precision 
observables measured in the LEP / LEP-II experiments to very high accuracy and 
tending to exclude such models if no special care is taken. 

In the last years, extra dimensional models have been proposed which can 
evade these constraints by delocalizing the Standard Model fermions within the 
extra dimension, thus allowing to tune the couplings to the new resonances in 
order to avoid these constraints. This way, such models are a viable method of 
breaking the electroweak symmetry and retaining perturbative TeV scale unitar- 
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ity without introducing a fundamental Higgs field. However, extra dimensional 
models (excluding trivial cases) are intrinsically nonrenormalizable and valid only 
below a cutoff scale, with most of the new resonances lying in fact above the cut- 
off. Conceptionally a honest extension of the Standard Model should only contain 
the structure below this cutoff, incorporating the extra dimensional mechanism of 
breaking the symmetry and delaying unitarity violation without making assump- 
tions on the high energy physics above the cutoff scale. 

The Three-Site Higgsless Model is a minimal implementation of this idea. While 
it can be motivated by extra dimensional Higgsless models of electroweak symme- 
try breaking, it in fact contains only one set of extra resonances which lies below 
the cutoff, delaying unitarity violation to « 2 - 3 TeV. The non-Standard Model 
part of the spectrum consists of a set of heavy partners for all Standard Model par- 
ticles with the exception of photon and gluon. The analysis of the experimental 
constraints reveals that, while the model is consistent with the precision observ- 
ables, the couplings between the new heavy gauge bosons and the Standard Model 
fermions have to be exceedingly small 1% of the isospin gauge coupling) while 
the new fermions are constrained to be rather heavy with masses above 1.8 TeV. 

In this thesis, we explored the LHC phenomenology of this scenario. To this 
end, we calculated the couplings and widths of all the new particles and imple- 
mented the model into the Monte-Carlo eventgenerator and WHIZARD / O'Mega. 
With this implementation, we simulated the parton-level production of the gauge 
boson and fermion partners in different channels possibly suitable for their discov- 
ery at the LHC. The results are presented in this work together with an introduction 
to the model and a discussion of the properties and couplings of the model. 

We find that, while the fermiophobic nature of the new heavy gauge bosons 
does make them intrinsically difficult to observe at a collider, the LHC should be 
able to establish the existence of both resonances and even give some hints about 
the properties of their couplings which would be a vital test of the consistency of 
such a scenario. For the heavy fermions, we find that their large mass is accom- 
panied by relative widths of more than 10%, making them ill-suited for a direct 
discovery at the LHC. Nevertheless, our simulations reveal that there is a part of 
parameter space where, given enough time, patience and a good understanding of 
detector and backgrounds, a direct discovery might be possible. 
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Chapter 1 

Introduction 



Of course our model has too many arbitrary features for these predictions to be 
taken very seriously... 

(Steven Weinberg, "A theory of leptons") 

This quote taken from U on what would later become the Standard Model of el- 
ementary particle physics may well be one of the most spectacular cases of false 
modesty in the history of science. Only six years after Steven Weinberg wrote these 
memorable lines on the then speculative unification of electromagnetic and charged 
current interactions in a spontaneously broken SU(2)l x U(1)y gauge theory, the 
discovery of the neutral current interaction at the Gargamelle bubble chamber in 
1973 at CERN suggested that there was more to this model than its authors origi- 
nally had dared to hope. 

Attempting to include the three known quarks into this model of weak interac- 
tions and explain the smallness of flavor violation, Glashow, Iliopoulos and Maiani 
postulated an additional quark (the charm) in 1970 GO, while in 1972 Kobayashi 
and Maskawa suggested another doublet of quarks (top and bottom) in order to 
explain the violation of CP symmetry observed in the hadron sector |3]|. Again, 
nature chose to prove the theorists right, and in 1974, the charm quark was dis- 
covered in form of the J/ ^ charmonium state, while the T bottomonium state was 
discovered in 1977. 

By the end of the 1970s, the Nobel prize committee in Stockholm deemed the 
evidence for the unified theory of electroweak interactions to be convincing enough 
to win Glashow, Salam and Weinberg the 1979 Nobel prize in physics, even before 
the postulated W and Z bosons were discovered as resonances at the UA1 and UA2 
detectors in 1983, with the 1984 physics Nobel prize going to Rubbia and van der 
Meer for this discovery. 

Up to this day, the success story of the Standard Model has continued with the 
impressive confirmation of the one loop structure of the model by the LEP / LEP-II 
experiments, the discovery of the top quark in 1995 at the Tevatron and the dis- 
covery of tau neutrino in 2000. Indeed, quite contrary to Weinbergs initial feelings, 
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every single prediction of the Standard Model has been confirmed in experiment, 
and all significant discoveries which are usually called "new physics" (e.g. neu- 
trino oscillations) can be easily accommodated within the model. 

The last remaining specimen from the Standard Model particle zoo is the Higgs 
boson. If it exists, then it must have managed to escape detection both at LEP / 
LEP-II and Tevatron, and together with indirect bounds on its mass coming from 
a global fit to the Standard Model, this implies that it most probably is lurking 
around 120 GeV and in any case must be lighter than several hundred GeV where 
the LHC together with the ATLAS and CMS experiments would surely find it. After 
all the past success of the model, it is arguably not the discovery but the absence 
of a Higgs resonance at the LHC which would be exciting news from the land of 
particle physics. 

Still, there is much truth in Weinbergs statement (after all, it comes from a No- 
bel laureate), and the Standard Model does exhibit many features that deserve to be 
called arbitrary and seem to be unfitting for a deep model of nature, e.g. the large 
number of free parameters (especially the Yukawa sector), the seemingly magical 
cancellation of gauge anomalies which is essential for the celebrated renormaliz- 
ability of the model or the three generations of fermions, each of which is an exact 
copy of the others, differing only in mass. 

Also, at the very least, the Standard Model does not describe dark matter and is 
incompatible with general relativity, and this fact alone implies that it can be only 
an effective description of nature. If we accept this as fact and still want to keep the 
Standard Model as-is up to the Planck scale, then we must answer the question why 
nature would adjust the Higgs mass at the matching scale to ridiculous precision 
in order to achieve the moderate value we observe. This is the famous Hierarchy 
Problem which would leave an ugly aftertaste if the Higgs were the only "new" 
piece of nature waiting to be discovered at the LHC. 

So, what happens if the LHC fails to discover the Higgs and puts an end to 
the success of the Standard Model? If there is no Higgs boson, then (unless our 
understanding of quantum field theory is fundamentally flawed) consistency con- 
ditions require that either new physics enters the game below 1 TeV or that the 
model becomes strongly interacting at this scale. In this sense, the LHC cannot fail 
in its mission: Higgs or not, something yet undiscovered most certainly has to lurk 
within its reach. 

Over the last decades, particle physicists have invested a lot of time and thought 
into finding out what such a discovery might look like. The resulting models can 
be divided into two categories: renormalizable models which (even if they may not 
be candidates for a fundamental theory of nature due to their physics content) can 
be valid to arbitrarily high scales, and effective field theories which are only valid 
below a cut-off scale and which do not claim to cover physics at scales significantly 
higher than those probed by the LHC. 

Models which fall into the first category are often motivated in a top-down 
fashion by theoretical concepts of how nature might organize itself at hight scales, 
for example supersymmetry or technicolor, while effective models are usually built 
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in a more pragmatic way as an extension of the low-energy Standard Model physics 
with little or no assumptions on the high energy structure of nature. 

Either way, it turns out that building a model which is in agreement with all 
experimental constraints collected over the last decades is quite challenging. In 
particular, any model of electroweak symmetry breaking different from the usual 
Higgs mechanism tends to change the structure of the leptonic current-current in- 
teractions which have been very precisely mapped out in the LEP / LEP-II exper- 
iments. The major part of the resulting constraints can be summarized in only 
three numbers which, however, have a devastating effect on many models of elec- 
troweak symmetry breaking. In particular, technicolor models and models of Hig- 
gsless symmetry breaking from compact extra dimension tend to violate these pre- 
cision constraints badly, and extra care has to be taken if models of this kind are to 
remain viable candidates for a description of nature^. 

In the specific case of compact extra dimensions, all particles usually come with 
a tower of partners of increasing mass which correspond to the higher excitations of 
the 5th momentum component (exceptions are particles which are explicitly local- 
ized on four-dimensional submanifolds by constuction), the so called Kaluza-Klein 
(KK) modes. Depending on the conditions enforced on the fields on the boundary 
of the extra dimensions, massless modes can be forbidden for the towers, thus facil- 
itating the breaking of electroweak gauge symmetry without introducing a Higgs. 

As stated above, if a model of new physics is to remain perturbative at the 
TeV scale, then it must contain some agent which retains perturbative unitarity at 
high scales, this part being played by the Higgs in the Standard Model. In Hig- 
gsless extra-dimensional models, this is achieved via the exchange of higher KK 
modes, which requires nonzero couplings between the Standard Model particles 
and their KK partners. However, these couplings also imply new contributions to 
the current-current interactions which lead to the aforementioned conflict with the 
LEP observables. Luckily, these constraints can be evaded by a moderate tuning of 
the wavefunctions of the Standard Model fermions, and this way, extra dimensions 
remain a potential candidate for the mechanism of electroweak symmetry break- 
ing- 

However, if we accept this fine-tuning of the parameters in order to comply 
with the precision observables, extra dimensional models contain another concep- 
tional flaw. Although they are built in a top-down way with the high energy struc- 
ture of spacetime in mind, the resulting models are in fact nonrenormalizable with 
a UV cutoff scale which is typically of the order of 5 — 10 TeV. Therefore, most of 
the new physics in these models is meaningless as it lies above the cutoff, and an 
effective field theory which only describes the KK modes below the cutoff would 
be a more honest way to describe the remaining sector which is not dominated by 
other high energy physics. 

The Three-Site Higgsless Model |4], the LHC phenomenology of which is the 

1 Of course, if the celebrated AdS/CFT correspondence is correct, then the difficulties technicolor 
and 5D theories are faced with are connected by this duality, and solutions in one formalism induce 
solutions in the other framework. 
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topic of this thesis, is an implementation of this idea. While it is clearly inspired 
by such Higgsless extra dimensional models of electroweak symmetry breaking, 
it only describes the Standard Model particles and the first generation of KK exci- 
tations, making no assumptions on the underlying high energy physics. The pre- 
cision constraints severely restrict the couplings between the new heavy particles, 
and we will find that, albeit the structure of the model is simple compared to other, 
more ambitious extensions of the Standard Model, these restrictions lead to a quite 
interesting and rich LHC phenomenology. 

In chapter [2J we motivate the model in a bottom-up approach, starting with the 
assumption of unitarity being maintained by W and Z partners in the absence of a 
Higgs boson and show that, while the model doesn't cure the arbitrariness in the 
Standard Model, most of the new structure is a straightforward consequence of this 
assumption. The explicit construction of the model using a more ideological top- 
down approach inspired by extra dimensions is shown in chapter [3J and in chapter 
H]we will move on to discussing the bounds on the model, its parameter space, the 
couplings and the widths of the new particles. Chapter|5]elaborates on the tools we 
have been using for the phenomenology study, the results of which are presented 
in chapters |6]-|8l 



Chapter 2 



The Model 
Approach 



Bottom-up 



I have a cunning plan. 



(Archetypical quote from Baldrick in "Blackadder") 



The Three Site Model is only one example out of a multitude of models for physics 
beyond the Standard Model that have been proposed. Before examining the phe- 
nomenology of any such model in detail, we should have a clear notion why look- 
ing at this particular kind of new physics is worthwhile. 

The purpose of this chapter is to give a bottom-up type motivation for the Three 
Site Model as a straightforward option for pushing the unitarity bounds of the Stan- 
dard Model once the notion of a fundamental Higgs field that unitarizes scattering 
amplitudes is discarded. 

Although this chapter contains some general remarks on quantum field the- 
ories, familiarity of the reader with the formalism and with the Standard Model 
ist assumed. These topics can be found covered in-depth in many textbooks, e.g. 

2.1 Restrictions on Quantum Field Theories 

Setting the stage 

According to our current understanding of nature, the world at very small scales is 
properly described in the language of relativistic quantum field theories. To define 
such a theory, we need to specify the field content of the theory as well as the 
Lagrangian. 

The field content tells us what different kinds of particles are described by the 
theory. The fields transform in different representations of the Poincare group 
which are distinguished by spin and mass of the field. 
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2.1 RESTRICTIONS ON QUANTUM FIELD THEORIES 



The Lagrangian C of the theory is a Lorentz invariant local function of the fields 
and their derivatives and determines the time evolution of the states described by 
the theory. It is commonly assumed that C either is a polynomial of the fields and 
their derivatives or at least can be expanded in a series in a meaningful way. 

The most useful known way to extract predictions out of this kind of description 
of nature is the use of perturbation theory to describe scattering experiments. Such 
experiments are modeled in an idealized way as a world that only consists of two 
particles with well-defined momenta in the asymptotic past, which then interact 
and scatter to finally evolve into a set of particles with well defined momenta in the 
distant future. This description relies on the hypothesis that there are such asymp- 
totically free states and that all interactions between the particles are negligible at 
asymptotic times ("adiabatic switching"). The time evolution operator that medi- 
ates the transition from the free particle states in the asymptotic past ("in states") 
to those in the asymptotic future ("out states") is called the S matrix operator. 

To describe the scattering process at finite times, the Lagrangian is split into a 
quadratic propagation part £f ree and an interaction part £i nt that consists of mono- 
mials of higher order in the fields^ • It is then assumed that the interaction part is 
just a small perturbation on top of the free propagation, and time dependent per- 
turbation theory is applied to calculate the transition amplitudes order for order 
as an expansion in the coupling constants appearing as prefactors at the monomi- 
als in £i nt . This formalism then leads to the well-known expansion of scattering 
amplitudes in terms of Feynman diagrams. 

Although this setup suggests a lot of freedom in the construction of a quantum 
field theory, there are some properties any sensible theory has to fulfill. The severity 
of these constraints depends on whether we want our theory to be acceptable as a 
fundamental description of nature or whether we just want it to be an effective 
description of our world that breaks down at some energy scale. 

Renormalization and renormalizability 

Any physical model comes with a number of free parameters that have to be fixed 
before the theory can be used to quantitatively predict experimental results. One 
criterion for a sensible model is the existence of predictions that can be checked to 
potentially falsify the theory, and therefore, these free parameters must be such that 
they can be determined with a finite number of measurements. 

When we set out to calculate higher orders in the perturbation expansion, we 
encounter loop diagrams that translate to unbounded integrals over four-momenta. 
Quite contrary to the assumption that higher order contributions are small com- 
pared to the leading order, the vast majority of these integrals is infinite. To study 
these infinities, it is useful to bound the integration domain and cut off the high 
energy contributions at some energy scale using a regularization scheme (e.g. a 



1 Because of the Lorentz invariance of the Lagrangian, linear terms in the fields are only allowed 
for scalars and generate a shift in the field vacuum expectation value in this case. 
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naive cutoff, Pauli-Villars, dimensional regularization etc.). Calculations in the reg- 
ularized theory then give finite predictions for observables which now carry an 
explicit dependence on the parameters and on the cutoff scale (in which they are 
divergent). 

Interestingly, there is a class of theories in which the reexpression of the free 
parameters in terms of a finite set of measurable defining observable^ leads to a 
cancellation of the divergences. In this case, the cutoff scale can be taken to infinity, 
leading to well-defined, finite observables that now only depend on a finite set of 
numbers that must be measured to completely fix the theory. This process is called 
renormalization, and theories which can be treated this way are called renormaliz- 
able. 

Of course, there is a lot of freedom in the choice of defining observables. Dif- 
ferent observables will lead to different convergence properties of the perturbation 
series. In particular, for a given observable, a clever choice of the scale \i at which 
the defining observables are measured can significantly reduce the contributions of 
higher orders; this corresponds to a resummation of the perturbation series. 

In the limit of infinite cutoff scale, the parameters which are expressed through 
the defining observables are divergent. Therefore, the above renormalization trick 
can be equivalently achieved by splitting the parameters g into a finite part §r and 
an infinite part 5g which is called a counterterm. The finite part is fixed by measure- 
ment while the counterterm is chosen according to some prescription such that it 
cancels the infinities. If the splitting prescription depends on some renormalization 
scale // like in MS, then a change of this scale shifts contributions between g^ and 
5g, changing the value of gR. This dependence of the gR on fj, is called the renor- 
malization group flow of the parameters, and a suitable choice of [i again allows to 
resum the perturbation series and reduce the contributions from higher orders. 

The whole process of renormalization can be interpreted as absorbing the effects 
of physics far above the renormalization scale into an (albeit infinite) redefinition 
of the couplings. 

Effective field theories 

Renormalizable field theories can be taken to be valid at arbitrarily high energies 
and therefore are candidates for a fundamental theory of nature. However, with 
any theory, the claim of it being fundamental is a very ambitious one and, further- 
more, the renormalizability condition turns out to be very restrictive on the types 
of allowed Lagrangians. 

However, if we drop the requirement of the theory being fundamental and in- 
stead think of it in terms of an effective field theory that is only valid below some 
scale A, then there is nothing bad in the idea of having a cutoff scale to make our 
integrals finite. Quite contrary, as the theory is only valid up to A, we should only 
include intermediate states below A anyway. Effects that come from high energy 

2 An introduction of a finite number of new free couplings that undergo the same treatment might 
also be necessary. 
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physics should be included as effective couplings in our low energy theory. The 
cutoff explicitly enters our calculations and again influences the convergence of the 
perturbation series. 

If we change the cutoff scale A, we can calculate how we have to adjust the 
couplings in front of the operators to retain the predictions. This leads to a renor- 
malization group flow of the couplings that can again be used to resum parts of the 
perturbation series and this way allow for a meaningful perturbative calculation of 
observables at different scales. The downside is that, as opposed to renormalizable 
theories, we must allow for infinitely many couplings as they will be introduced by 
the renormalization group flow anyway. 

To keep our theory predictive, we must assume some kind of ordering scheme 
where the higher order operators are suppressed by some scale A^p and therefore 
are negligible at low energies. If we construct an effective theory e.g. by integrating 
out heavy particles from a renormalizable theory, the new operators describing the 
effects of the degrees of freedom that have been removed obey exactly such an 
ordering scheme, so this is not an unnatural assumption. This way the theory can 
be predictive at low energies even though the Lagrangian describing it contains 
arbitrarily many operators. 

If we want to calculate observables at some scale A, we must evolve the cou- 
plings to this scale in order to perform the perturbative calculation. The higher 
order interactions increase if we evolve to higher scales, and at some scale Ao, the 
ordering scheme will break down. At this scale the very latest, we have no con- 
trol over the contributions from higher order operators anymore, and the effective 
theory is not predictive anymore. It then has to be replaced by another theory that 
might contain new degrees of freedom and which can either be a fundamental, 
renormalizable theory or again an effective theory. 

This way, if we don't require our theory to be fundamental, renormalizability 
is not necessary for the theory to make sense. We only have to be aware that the 
energy range in which the theory is valid is limited and that we have to replace it 
by something else above Ao- Not insisting on renormalizability allows for a much 
larger class of theories which can be effective descriptions of nature even if quan- 
tum field theories should turn out not to be the suitable language for the "theory of 
everything". Even better, the scale Ao as well as a self-consistent ordering scheme 
for the higher dimension operators can be estimated from the parameters of the 
theory by simple dimensional arguments without requiring any knowledge about 
the underlying high energy physics; this procedure is called NDA [9J ("naive di- 
mensional analysis"). 

The remainder of this chapter will be dealing with effective field theories, so 
renormalizability will not be an issue as long as we stay in the energy range where 
these theories are valid. 
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Unitarity — limits on the validity of the perturbation expansion 

Unrelated to the issue of renormalizability, there is another criterion for the con- 
sistency of field theories. Because it is the limit of a time evolution operator, the 
axioms of quantum mechanics require the S matrix to be unitary. This is the pre- 
requisite for the probabilistic interpretation of quantum mechanics. 

The diagonal elements of the S matrix contain contributions from the case of 
direct propagation of the initial into the final state without any scattering. As we 
are usually not interested in these contributions, the cross section is conveniently 
defined as a flux normalization factor times the modulus of an T matrix element, T 
being defined as 

S = iT + I (2.1) 
The unitarity condition on S then implies a relation for T: 

SS^ = I — ► i (V - T) = TT ] (2.2) 

Taking the matrix element between two states \a) and \b) and inserting a complete 
set of basis vectors, we obtain 

i ((b\ T \a)* - (a\ T \b)) = £ <a| T \c) (b\ T \c)* (2.3) 



From (I2.3D , important consequences for matrix elements can be derived. Re- 
stricting the equation to the forward scattering case |o) = \b) yields 

2%(a\T\a) = J2MT\c}\ 2 (2.4) 

\c) 

This is the famous optical theorem which basically states that the imaginary part 
of the forward scattering amplitude of some state is equal to the probability for this 
state scattering into anything. If |o) is an eigenstate of T with eigenvalue T\a) = 
1 1 a) and norm (a \ a) = 1, (|2.4|) further reduces to 

29t = |t| 2 

which is solved b)H 

t = i + e i5 (2.5) 

(I2.5D restricts the eigenvalue t to a "Argand circle" in the complex plane with center 
i and radius 1. 

Due to rotational invariance, two particle angular momentum eigenstates be- 
low the inelastic threshold (the regime in which the total energy is not sufficient to 
create any new particles) are eigenstates of the T matrix. If we examine the cross 



3 The same result can also be directly obtained from the definition of T in terms of the unitary S 
matrix J2.ll . 
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section between incident momentum eigenstates, then these states can be decom- 
posed in a relativistic partial wave expansion flQ|, and (|2.5|) implies relations for 
the partial wave amplitudes. If we cross the inelastic threshold (I2.5D changes into 
an inequality, and the amplitudes are constrained to lie within the Argand circle 
rather than on it. For a more complete discussion of unitarity and partial wave 
amplitudes see lfT0ll6fl 

If the Lagrangian is hermitian, then the S matrix will be unitary by definitional 
However, as (12.211 - (12.511 are nonlinear in the T matrix elements, they don't hold 
order by order if the T matrix elements are expanded in a perturbation series^. 
Therefore, at a fixed order in perturbation theory, the partial wave amplitudes can 
violate the unitarity criteria. Worse, if an amplitude grows with energy in pertur- 
bation theory, then it will leave the Argand circle at some scale and higher orders 
must become increasingly large to compensate and restore consistency to the the- 
ory. 

In particular, if the tree level amplitude for some scattering process grows with 
energy, then there will be a scale where the tree level scattering probability ex- 
ceeds one and at which perturbation theory stops to be a trustworthy calculational 
instrument and breaks down. If we write down an effective field theory, this uni- 
tarity cutoff will limit the range of applicability of the theory independently from 
the NDA cutoff. 



2.2 The Fermi Model 

At "low" energies much smaller than the W mass, our world can be well described 
by the Fermi model. At this scale, the observable fundamental degrees of freedom 
are|§ 

• Three generations of leptons e~ , r~ which differ only by their mass and 
which are described by Dirac spinors Lj . 

• The associated three generations of neutrinos v e , v ^ , v T which we will approx- 
imate to be massless and which are described by Dirac spinors Nj . 

• The photon which mediates the electromagnetic interaction and which is de- 
scribed by the vector field A^. 

The Fermi Lagrangian ^Fermi can be decomposed as 

^Fermi = £qED + ^4 

4 Exceptions can arise if the quantization is done incorrectly and states are inconsistently removed 
from the Hilbert space, e.g. due to improper quantization of gauge theories. 

5 Of course, the equations themselves can be expanded to hold order by order, but in this case, 
different orders in the expansion of the matrix elements will mix. 

6 As color and QCD are not related in any way to the line of reasoning aimed at in this chapter, all 
colored particles will be ignored for the remainder. 
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The first part is the usual QED Lagrangian £qed given by 

3 — 1 

£qed = ( L 3 W + m i) L o + iNtfNj) - eA„J% - -F^F^ (2.6) 

i=i 

with the electromagnetic gauge coupling e, the electromagnetic current 

3 

and the electromagnetic field strength tensor 

Defining the left- and right-handed fermion fields * i and via l|A.3|l , we can 
introduce the charged and neutral currents 

1 3 - 

^ = Jp (2-8) 
J ° = E (sm 2 ^ - ^n_) Lj + ^A/ .r^.V,., 

(with the Weinberg angle 9w) an d write the second piece £4 as 

£4 = -4^/2G F (J+ + J°J 0/1 ) (2.9) 

where the Fermi coupling Gp has mass dimension [Gp] = —2. 

(|2.6|l contains the bilinear parts of the Lagrangian that generate the propagators 
as well as the three point gauge couplings of the photon to two leptons, while (I2.9D 
encodes the charged and neutral current four point couplings among four fermions. 
Apart from the lepton masses, the only free parameters are the electromagnetic 
gauge coupling e, the Fermi coupling Gf and the Weinberg angle 9w- Once these 
parameters are fixed, the Fermi model does a very good job at describing the low 
energy phenomenology of our world. 

However, while pure QED is renormalizable, the four point interactions make 
the Fermi model a nonrenormalizable model that is only valid up to a UV cutoff 
scale Auv- To estimate this scale, consider the one loop correction to the charged 
current four point interactior0: 





*>o<* 



7 The t channel contribution can be avoided by considering the four point function for two differ- 
ent flavors (e.g. e~v e — > liTVp). 
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Naive powercounting reveals a quadratic divergence in the one loop contribution, 
and we can estimate the leading contribution by 

/d A n "M 4 Muv 2 

— ^ TV VL = G F 2 N f / dp p = -G F 2 Al v N f 
(2tt) 4 P 4 tt 2 Jo vr 2 

where Nj is the number of flavors which can run in the loop. Considering only 
the leptons we have Nf = 3, and demanding that the one loop correction must not 
exceed the tree level contribution if perturbation theory is to make sense then leads 
to a estimate for the cutoff scale 



Of course, if we had also included the quarks into our version of Fermi theory, 
they would appear within the loop, and the cutoff would be lowered further to 

A uv « 315 GeV. 

However, even when setting aside all issues related to renormalization, a straight- 
forward calculation of the amplitude for // — > ff type processes reveals terms 
induced by the four point couplings which grow quadratically with energy. Partial 
wave analysis yields a scale of 600 GeV [TTl [12| at which the s wave amplitude 
exceeds the unitarity bound. At this scale the very latest, perturbation theory can't 
be trusted anymore. If our world is to be described by a perturbative quantum 
field theory, then new physics must come into play below this scale and unitarize 
the scattering amplitudes or at least delay the scale of unitarity violation. 



2.3 The Standard Model 



Indeed, experimentalists didn't have to go to scales of 600 GeV to discover new 
physics that supersedes the Fermi model: in the early eighties, the and Z 
bosons were discovered at the Super Proton Synchrotron SPS at CERN and at Fer- 
milab with masses of mw ~ 80 GeV and mz ~ 91 GeV IIT31 . 

Together with the photon, the and Z are understood as the gauge bosons 
of a SU(2) l x U(l)y gauge group which unifies electromagnetism and weak inter- 
actions. In the following, we will denote the SU(2)l gauge fields as = W^t^ 
(with the SU(2) generators r^, see dA.21) ) and the U(l)y gauge field as -B M ; the 
gauge couplings will be called g and g' respectively. Introducing the field strength 
tensors 

F% = d^W u - d u W^ - ig [W, W u ] , = d»B v - d u B» 

the gauge kinetic Lagrangian can be written as 

£gauge = — xTV Fy^ FyVfiu — T-^B ^Bfiv 
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As the W and Z are massive, the gauge symmetry must be broken, and the 
only known mechanism which achieves this while still allowing for a consistent 
quantization is spontaneous symmetry breaking. Without any assumptions on the 
underlying dynamics, the breaking can be parameterized by the introduction of a 
SU(2) valued fielcfi 

E = ve l v Tk 

where v is the vacuum expectation valu^l of E 

(E) = vl 

This field is assigned the transformation behavior 

s SU(2) £ xU(l) y; ^ k r k ^X Y r, 

with the gauge transformation parameters \L,k, Ay. With this charge assignment, 
the covariant derivative of E reads (c.f. ((AT)) 

D^T, = d^E - igW^T k E + ig'B^n 

The minimal gauge invariant Lagrangian required to make E a dynamic field is 

£ s = Tr (£PE) (£> M E) f (2.10) 

Expanding (|2.10|) in the component fields 4>i, we find that the lowest order (aka the 
vacuum expectation value of E) leads to mass terms for the gauge bosons which 
can be written as 

t (» 2 t k*» + w,. *,> (4 ( ^ ) ) (2.1D 

(I2.11D is diagonalized by the mass eigenstates 

± = ^T^ 2 ,, gWs^g, g^+gg, 

and yields the Z and masses as eigenvalues 



m w = gv , m z = \J g 2 + g' 2 v (2.13) 

Classically, the E field can be put to vl at every point in spacetime by a gauge 
transformation, absorbing the three corresponding degrees of freedom into the 
gauge fields. This identifies the three 4>k fields as the Goldstone boson fields which 
are required by the Goldstone theorem [14J and which are "eaten" by the heavy 



8 More accurately, it is the rescaled field which is valued in SU(2). 

9 Note that this definition of v differs by a factor of \[2 from that usually used found textbooks. 
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gauge bosons to constitute the longitudinal degrees of freedom. As the 0k trans- 
form in a nonlinear realization of the gauge group, £ is called a nonlinear sigma 
field. 

To reproduce the physics of the Fermi model, we have to specify the represen- 
tation in which the leptons transform. The left-handed leptons of each generation 
are grouped in isospin doublets 



L L,j 

that transform in the fundamental representation of the SU(2) l, the right-handed 
leptons are singlets under this group factor. The hypercharge transformation prop- 
erties are derived from the electric charge of the fermions: 

with the electric charge Q. These assignments completely fix the covariant deriva- 
tive and the conserved currents which couple to the gauge bosons. Defining the 
electromagnetic gauge coupling e and the Weinberg angle 9\y as 

■ 1<J smO w = -=jL= (2.14) 



the resulting fermionic Lagrangian including the kinetic terms can be written as 

3 

^fermion = « ^ ( L k$ L k + N k 0N k ) + 

i=i 

~ ~n + W-J^) + -^—Z»J Q , - eA»J Qfl (2.15) 

The conserved currents appearing in (12.151) are the very same currents already de- 
fined in (23 and tZM . 

In the Fermi model, the fermion mass terms can be inserted into the Lagrangian 
directly. However, with the above SU(2)l x U(l)y charge assignments, such terms 
would now violate gauge invariance. Fortunately, the nonlinear sigma field £ al- 
lows to generate the necessary terms in a gauge invariant way from Yukawa cou- 
plings© 

Yukawa = E^ S (o (2 - 16) 

Putting together all pieces of the puzzle, the full Lagrangian of the theory reads 

£ = -£gauge + ^fermion + Yukawa + (2.17) 



1 For the quark masses, the mass matrix in J2.161 is replaced by ^diag (m u ,m d ) 
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The resulting model is the Standard Model, but without any assumptions on the 
dynamics of the electroweak symmetry breaking and therefore without the Higgs 
field. 

The Standard Model without a Higgs is nonrenormalizable due to the dimen- 
sion 6+ operators appearing in the expansion of the £ fields and therefore must be 
treated as an effective field theor y. E ven without knowing the physics responsible 
for the symmetry breaking, NDAlj can be used to get an upper bound on the UV 
cutoff 0: 

A NDA = 8ttv = 8vr— ra 3.1 TeV (2.18) 
9 

How does the transition from the Fermi model to the Higgsless Standard Model 
affect tree level unitarity? In the Fermi model, the divergent parts violating uni- 
tarity in the // — > ff amplitude stem from the four point fermion coupling. 
However, in the Higgsless Standard Model, these couplings are replaced by the 
exchange of the W and Z bosons: 

x->-<*X 

At low energies, the W ± and Z can be integrated out by replacing the propagators 
by — -^j to obtain the Fermi theory as an effective theory. However, above the W' 
mass scale, the propagators and the amplitude fall off as -K. This damping removes 
the quadratic growth in energy from the partial wave amplitudes and eliminates 
the unitarity violating terms. 

The entrance of the and Z bosons on the stage restores perturbative unitar- 
ity to the amplitudes for four fermion processes. However, the exorcism of unitar- 
ity violation is not complete — perturbative unitarity is still violated by amplitudes 
containing external heavy vectors. To understand this, recall the high energy limit 
of the longitudinal polarization vector 

a E—>oo 

^ 



At high energies, e l behaves like the four-momentum and therefore, the amplitude 
for the scattering of longitudinal vector bosons 




might develop a high energy behavior as bad as E 4 . Explicit calculation |15| shows 
that, while the E A parts cancel, a quadratic divergence remains which leads to a 



12.181 differs from |9| by a factor of 2 due to a difference in convention. 
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violation of tree level unitarity at approximately 1 TeV, well below the NDA cut- 
off. Therefore, new physics that unitarizes the scattering amplitudes should be 
expected below 1 TeV if the theory is to remain perturbative. 



2.4 The Higgs 

The question of unitarity at the TeV scale is not the only terra incongnita not cov- 
ered by the Higgsless version of the Standard Model. In particular, we didn't make 
any assumptions on the physics responsible for the dynamics of electroweak sym- 
metry breaking. As simplicity always has been an effective guiding principle in 
physics, it is not unreasonable to assume that these questions have the same an- 
swer and that the physics responsible for breaking the electroweak symmetry also 
unitarizes the scattering amplitudes. This is exactly what happens in the Standard 
Model: a scalar Higgs field is introduced which dynamically develops the vacuum 
expectation value which breaks the electroweak symmetry and which at the same 
time unitarizes the scattering amplitudes. 

The Higgs can be incorporated into the Higgsless Standard Model presented 
in the last section by replacing the vacuum expectation value of the £ field with a 
dynamical field H(x) 

S = ve^ k(x)Tk — > S = H(x)e^ k{x)Tk 
and adding the potential 

C H = -V H (x) = ^Tr SEt - jTr (sijt) * = 2 - ±H A (2.19) 

with the "mass" [i and the coupling constant A to the Lagrangian. This potential is 
the famous Mexican Hat potential whose nontrivial minimum constitutes the tree 
level vacuum expectation of the Higgs 

Shifting H(x) in order to perform the perturbation expansion around the ground 
state, we finally obtain Cr and £ in terms of the physical Higgs field 

£ = (v + h{x)) e^ttW 1 * 

£h = ^-t- ^ 2 h 2 - 2fiV\h 3 - -h 4 (2.20) 
2A 2 

Inserting the modified S into (12.171) and adding d2.20[) finally leads to the usual 

Standard Model. The nonlinear representation of the Higgs and Goldstone fields 

may look unfamiliar compared to the linear doublet representation usually used in 

textbooks, but both formulations can be transformed into each other by means of a 

nonlinear field redefinition and therefore are equivalent |16|. 
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What has changed with the introduction of the Higgs? In the Higgsless ver- 
sion of the Standard Model, the symmetry breaking was just parameterized by the 
introduction of the Goldstone bosons, but the physics responsible for the break- 
ing was not included and the vacuum expectation value v was fixed by hand. In 
the version of the model including the Higgs, the vacuum expectation value arises 
dynamically from the Higgs potential, and the observable consequence is the ap- 
pearance of a physical scalar in the particle spectrum. 

The Higgs couples to every massive field and, in particular, new contributions 
to the scattering amplitude for longitudinal gauge bosons arise 



+ 1 

As already mentioned above, explicit calculation shows that these new diagrams 
completely cancel the divergent pieces and, if the Higgs is not too heavy, perturba- 
tive unitarity is restored at arbitrarily high scales [15 J. 

In addition to resolving the unitarity problem, the inclusion of the Higgs makes 
the Standard Model a renormalizable theory [17J, removing the remaining obstacle 
in interpreting it as a fundamental theory valid at arbitrarily high energies^- From 
this, one might draw the conclusion that the Higgs is the preferred candidate for 
the mechanism of electroweak symmetry breaking. 

However, the Standard Model does not include gravity, and therefore, it cannot 
be a fundamental theory of nature but must be superseded by new physics at some 
(although possibly high) scale. Once we have accepted this fact, the Higgs looses 
much of its appeal through the so-called Hierarchy Problem. Consider the one loop 
corrections to the Higgs mass 

-O- + -0- +-■£-+-■&-- +--<>- ("I) 

Apart from the first two diagrams, all pieces of (12.211) contain quadratic diver- 
gences inducing a quadratic running of the Higgs mass. Therefore, if the Higgs is 
to have a mass mg at low scales, then the mass must be tuned with an accuracy of 
approximately at the scale Auv where the Standard Model is matched to some 
more fundamental description of nature. 

In order to preserve perturbative unitarity, the Higgs mass must be consider- 
ably smaller then 1 TeV. If we insert the Planck scale as Auv/ this implies a fine 
tuning to a rather ridiculous precision of roughly 1CT 14 %! While this extreme sen- 
sitivity to the underlying high energy physics is not technically inconsistent, it ar- 
guably seems unfitting for a theory of nature which is deep enough to be valid over 



12 Even if the model is renormalizable and tree-level unitary at all energies, Landau poles in the 
renormalization group flow might limit the range of applicability of the model, so this statement is a 
bit oversimplified. 
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19 orders of magnitude. Apart from the Hierarchy Problem, there is also a more 
mundane issue concerning the Higgs boson: no trace of it has been yet discovered 
in collider experiments and, looking at current exclusion plots (see e.g. |18|) for the 
Higgs mass, the air is growing thin. 

Considering the effective nature of the Standard Model in conjunction with the 
Hierarchy Problem and the reluctance of the Higgs to show up in experiments, the 
concept of a fundamental Higgs certainly drops in grace, and other mechanisms 
for electroweak symmetry breaking and/ or maintaining perturbative unitarity be- 
come attractive. 

2.5 An Alternative 

If the amplitudes for the scattering of longitudinal gauge bosons are to fulfill the 
tree level unitarity requirement, then their high energy behavior must be modified. 
This can be achieved either through modifying the propagators of the particles or 
by introducing new particles into the model, the exchange of which cancels the 
dangerous growth with energy. 

Of course, a scalar is not the only possible particle suitable for achieving this 
goal; an alternative ansatz would be the introduction of new vector bosons. These 
new vectors would have to be massive and rather heavy in order to have accom- 
plished the escape from detection until now. However, the only known case of 
massive spin 1 particles that can be consistently quantized is that of gauge bosons 
of a spontaneously broken symmetry. So, let us introduce an additional SU(2) 
gauge group with new heavy gauge bosons which we will call W' and Z' . To pa- 
rameterize the symmetry breaking, we will also need a new nonlinear sigma field 
Tl which describes the three new Goldstone bosons that arise from the symmetry 
breaking. 

If the new particles are to unitarize the scattering amplitudes, then they must 
couple to the and Z bosons. However, the couplings of gauge bosons are com- 
pletely determined by gauge invariance, and if we have only kinetic ancf^l mass 
terms, then we will only get couplings of the types 



(assigning double lines to the new heavy particles), but no couplings that mix the 
Standard Model gauge bosons and the new heavy ones. 

A way to lift this constraint is the introduction of terms into the mass matrix 
which mix the Standard Model gauge bosons and the new ones. As a results, the 

3 As this section focuses on the ideas behind the model, the discussion is qualitative and the details 
are postponed to the next chapter. 
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physical fields arise as mixtures of the old and the new fields 




) 




V W 3ll 

In 



with orthonormal mixing matrices U, V. This mixing induces the desired couplings, 
resulting in contributions to the longitudinal gauge boson scattering amplitude 



which can potentially restore perturbative unitarity. 

The mixing between the gauge bosons also induces couplings of the W'/ Z' to 
the standard model fermions, leading to contributions to the charged and neutral 
current four point function 



As new physics would also appear in loop contributions to the charged and neutral 
current correlators, this kind of processes was studied very precisely by the LEP-II 
experiments, resulting in an astonishing agreement with the one loop predictions 
of the Standard Model and leading to very stringent bounds on any kind of new 
physics which contributes to this kind of processes. The major part of this so-called 
electroweak precision data can be distilled into the famous aS, aT and all pa- 
rameters defined in |l9ll20| or equivalently the three ei/2/3 parameters defined in 



If no special care is taken to avoid these constraints, most models of new physics 
that introduce copies of the W and Z which couple to Standard Model fermions are 
excluded. In our construction so far, these couplings are completely determined by 
the mixing matrices U and V, and as we have no further influence on the result- 
ing couplings, we will almost certainly violate the electroweak precision bounds. 
Therefore, once we have introduced the W and Z' this way, we have to think of a 
trick to gain control over the undesired couplings. 

Of course, we could try to charge the Standard Model fermions under the new 
SU(2) gauge group, but the resulting couplings are completely fixed by gauge in- 
variance, giving us no additional freedom that could be used to avoid the con- 
straints. What could work, however, this the introduction of a partner fermion for 
each Standard Model fermion which is heavy, charged under the new SU(2) and 
mixes with its Standard Model partner. The price would be new set of new heavy 





mum. 
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fermions but we could use the mixing parameters to control the unwanted fer- 
mion couplings and tune them away. 

This way, we are rather naturally led to a model which contains a heavy copy 
of every Standard Model particle (with the exception of the gluon and the photon). 
There is no physical scalar in the spectrum and consequently there is no hierarchy 
problem. The task of unitarizing the scattering amplitudes is taken over by the new 
set of heavy vector bosons. 

This structure of massive replica of the spectrum is a well-known feature of 
theories with one or more compact extra dimensions. In such theories, gauge sym- 
metry can be broken without a Higgs by the introduction of suitable boundary 
conditions, and the massive gauge bosons that arise this way are known to de- 
lay the scale of unitarity violation II23II . Reducing the extra dimension to a finite 
point lattice, this can be exploited to build a model which has exactly the spectrum 
and features discussed above, the Three-Site Higgsless model |4|. The next chapter 
will be devoted to a detailed discussion of the construction and properties of this 
model. 



Chapter 3 

The Model — Top-down Approach 



Percy: Only this morning in the courtyard I saw a horse with two heads and two 
bodies. 

Blackadder: Two horses standing next to each other? 

("Blackadder I — Witchsmeller Pursuviant") 

In the last chapter the Three-Site Model was demonstrated to arise in a rather 
natural way when attempting to unitarize scattering amplitudes with heavy copies 
of the W and Z bosons instead of a Higgs. The actual construction was postponed 
and will now be carried out using a convenient framework. 

The first two sections are devoted to introducing the concepts of extra dimen- 
sions and dimensional deconstruction, while the third then continues with the ac- 
tual construction of the model. 

3.1 Extra Dimensions 

Manifolds and boundary conditions 

The first ingredient to an extra dimensional field theory is, obviously, the spacetime 
manifold on which it is formulated. As the theory must match onto our known 4D 
physics at low energies, this manifold must contain the usual Minkowski space- 
time (we are not going to consider general relativity), extended by one or more ex- 
tra dimensions. For the same matching reason, these additional dimensions must 
decouple from the low energy phenomenology and are therefore usually chosen to 
be compact with some compactification scale R which suppresses the new physics 
arising from them. 

In the following, we will specialize to a flat 5th dimension which (together 
with the case of "warped" extra dimensions with Randall-Sundrum metric 1 24 1) 
is arguably the best-studied case in particle phenomenology. To complete the ba- 
sic setup, we also need to specify the topology of the extra dimension. Common 
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choices are a circle with radius R or a compact interval [0; 2Rtt}. Out of these two 
examples, we now specialize to the latter. The boundaries of such an interval are 
also sometimes called "branes", while the open interval is called "bulk". 

The second ingredient to the theory is the set of fields that propagate on the 
manifold. In analogy to the 4D case, these form representations of the 5D Poincare 
grourU In 4D field theories, we require that the fields decay "fast enough" for 
asymptotic times and distances, and in a 5D theory, we need to specify boundary 
conditions on the boundaries of the compact extra dimension. On the interval, 
these are usually chosen either as Dirichlet or as Neumann boundary conditions^ 

Neumann: d y ${x, l/)|„ e{0)2Jfer} = 

Dirichlet: l/)| w6{0)2Jfer} = 

It is also possible to mix these conditions and e.g. impose a Dirichlet boundary 
condition at one end of the interval and a Neumann one on the other one. 



5D scalar field and Kaluza-Klein expansion 

The simplest possible 5D theory is that of a real scalar field on an interval: 

$ : (a?, t/)£R 4 x [0; 2Rir] — ► $(x, y) 6 R 

As an example for how to treat such theories, let us consider the Lagrangian for the 
flat 5D variant of <1> 4 theory (see appendix [A] for the conventions regarding Lorentz 
indices) 

A.* = l -{d a <$>){d a $>) - ht 2 ^ - |d> 4 (3.2) 

The action functional is then defined as the usual integral of the Lagrangian over 
spacetime, the y part of which is bounded 

Sm = J d*x J 2 *dy£& = J d^f^dv Q^.*)^*) - \m*& ~ f* 4 ) 

'(3.3) 

It is noteworthy that, in order for the action (I3.3D to be classically scale invariant, 
the mass dimensions of the scalar and of the coupling must be 

[*] = \ . b] = -i 

Variation of the action to obtain the equations of motion proceeds as usual. How- 
ever, the result 



SS = I d A x{ d y $5<i> 



y=2Rir r 2 ^ / a 

- dy I d a d a <S> + M 2 $ + ) 5$ 
y=o Jo ^ 3! 



1 Usually, the compactification breaks the 5D Poincare group down to the 4D one, resulting in a 
violation of 5D momentum conservation. 

2 Starting from a fifth dimension compactified on a circle, the case of compactification on an inter- 
val with Dirichlet or Neumann boundary conditions can also be obtained in a very elegant way by 
orbifolding (25). 
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contains a boundary term that results from partial integration. With the boundary 
conditions (|3.1[) or mixed ones, this term vanishes, and the resulting Euler-Lagrange 
equation is 

(d a d a + M 2 ) $ = (d 2 -d 2 y + M 2 )<5> = - J^ 3 (3.4) 

It is convenient to partially decompose the solutions to (13.41) as a Fourier series 
with the basis functions f k {y) chosen such that they satisfy the boundary conditions 
and are eigenf unctions of d 2 : 



1 ( k 

Neumann: $(x,y) = V cf> k {x)f k {y) = V^(x) - cos — y 



k=0 k=0 

oo oo 



oo oo 1 / k \ 

Dirichlet: $(x,y) = ^(t>k{x)fk{y) = ^ (j>k(x)-= sin ( — y J 

k = l fc=l V-K7T \ / 



(3.5) 

Inserting the expansion into the Lagrangian (I3.2D and defining the effective 4D 
action functional and Lagrangian as 

5[^>] = y <i 4 x £$4 , £$4 = j dy £$4 

we obtain 

£® 4 = \Y1 (( d M 2 - (m 2 + 4>l \ - »fci,fca,*s,*40*i0fca0*s0fc4 

(3.6) 

with the dimensionless coupling constants 

9 f 2Rn 

9ki,k2,k^,k4 — "77 / dy fk 1 {y)fk 2 (y)fk 3 (y)fk 4 (y) (3.7) 

This expansion is the so-called Kaluza-Klein (KK) expansion of the theory. 

What is the virtue of the KK expansion and of the effective Lagrangian (I3.6D ? 
While we started with a 5D theory of a single scalar field, we find that the physics 
can be equivalently described by an ordinary 4D field theory with an infinite num- 
ber of 4D scalar fields <j) k . Each of these fields has an effective mass of 

making them arbitrarily heavy for k > 0, the mass increasing with k. The 5D La- 
grangian £$4 contains a single four point interaction of the <!> whereas the effective 
4D version contains four point interactions involving all the different (f> k fields (al- 
though many of the coupling constants vanish). 
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The physics behind this is easy to understand: in the 5D version of the theory 
the 5-momentum satisfies the mass shell condition 

k a k a = k^ - (k 5 ) 2 = M 2 

with the fifth component of the momentum being discret^l rather than continuous 

h 5 - — 
n ~ 2R 

due to the compactification of the extra dimension. The straightforward interpre- 
tation is to think of <I> as a 5D field with a discrete fifth momentum component and 
mass M 2 . However, it can be equivalently regarded as a tower of 4D fields with 
masses given by M 2 + (k 5 ) . In the same way we could think about a 4D particle 
as a continuum of 3D particles, but as the four-momentum is a continuous degree 
of freedom, this is not a very useful intuition. 

At tree level, we have the liberty of setting the 5D mass M to zero. In the case 
of Neumann boundary conditions, there is a k = mode in the tower of scalars 
which is massless. In the case of Dirichlet or mixed boundary conditions the flat 
wavefunction fo(k) is forbidden, and all the 0& are massive, even if no tree level 
mass is put into the Lagrangian. 

Applying naive powercounting to the 5D version of the Lagrangian leads to 
the conclusion that it must describe a nonrenormalizable theory as can be easily 
understood by noting that the 5D coupling g has negative mass dimension. In the 
4D theory, the coupling constants at the four point interactions are dimensionless 
and therefore, the theory might appear to be renormalizable at a casual glance. 
However, in higher order diagrams, all the infinitely many (\>k can run in the loops 
leading to a possibly divergent series over k. This enhances the divergences and 
recovers the nonrenormalizable character of the theory. 

Although we have been examining just the simple case of a 5D <I> 4 theory, we 
have found a lot of interesting structure that generalizes to more complicated 5D 
theories. These theories can always be described by an effective 4D theory by virtue 
of the KK expansion. In this process, every 5D field gets replaced by a tower of 4D 
fields of increasing mass. Asymptotically, the masses of the modes in the tower are 
equidistantly spacec|^|. The presence of a zero mode which is massless in absence of 
an explicit mass term in the Lagrangian depends on the boundary conditions. Extra 
dimensional theories are usualhj! nonrenormalizable and valid only as effective 
field theories with an ultraviolet cutoff. Therefore, it is arguably inconsistent to 
include all modes in the KK towers in calculations. Instead, the towers should be 
cut off in some way at high scales, and if gauge fields are involved, extra care has 
to be taken to preserve gauge invariance in this step. 

3 To be a bit more precise, the fk are not eigenfunctions of id y but of d\, so we should rather be 
talking about the square of the five momentum and not about the momentum itself. 

4 This changes if the metric of is chosen nontrivially, e.g. in Randall-Sundrum type models. 

5 Trivial counterexamples exist, e.g. <f> 3 theory in 5 or 6 dimension (which, however, is a slightly 
pathological example as the potential is not bounded from below). 
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Gauge bosons 

A 5D gauge field transforms in the vector representation of the five dimensional 
Poincare group and hence has five components A a . However, boundary condi- 
tions break this group down to its 4D subgroup, under which transforms as a 
vector and A 5 transforms as a scalar. Once we perform the KK expansion to get the 
4D effective theory, the A^ will give us a tower of 4D vector fields, while A 5 will 
become a tower of 4D scalar s. In the following, we will restrict the discussion to 
the case of an abelian field for simplicity's sake. 

Under a gauge transformation with parameter field A, the fields A^ and A 5 
transform as 

A^x,y)^A^x,y) + -d^X(x,y) , A 5 (x, y) — ► A 5 (x, y) - -d y X(x, y) (3.8) 

9 9 

with the gauge coupling g. Similar to the scalar case, we have to choose boundary 
conditions for the fields. Let us choose Neumann boundary conditions for A^ . We 
can then KK expand the 4D vector field as 

00 1 / k \ 

A»(x, y) = Y A>j(x) - cos — y 



Looking at H3.8D , it is easy to see that the parameter field A has to satisfy the same 
boundary conditions as A^ and therefore can be KK expanded likewise 

X(x, y) = X k (x) — ; cos ( — -y | 

{ ' { V^d + M \ 2R ) 

This implies that the shift in A 5 generated by such a gauge transformation is 

I ™ k 1 ( k \ 
A^x, y) ^AHx,y) + -Y [ X k - 7 = sin {-y j 

which means that the field A 5 must obey Dirichlet boundary conditions and can be 
decomposed as 



A 5 (x,y) = J2A 5 k (x)-^ sin ^y^j 



Therefore, after the KK decomposition, we obtain an infinite tower of 4D gauge 
transformations which leave the theory invariant, the 4D gauge bosons A^ and the 
4D scalars A 5 transforming as 

A%(x) —>A%(x) + -dn k (x) k > 

1 k (3 - 9) 

A 5 k (x) — >A|(x) + -— X k (x) k > 1 

g IK 
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A closer look at 03.9D reveals at least two things of special interest. First, there is 
a gauge transformation generated by Ao which only acts on Aq and doesn't involve 
any of the 4D scalars At. Second, it is possible to perform gauge transformations 
which remove all of the A\ and completely absorb them into the At, fixing the 
gauge for all the 4D gauge fields with the exception of A 1 ^ . This suggests an elegant 
interpretation. 

As we chose Neumann boundary conditions for the A^, the fifth component of 
the 5-momentum fc 5 vanishes for Aq, leaving this mode as a massless gauge bo- 
son of an unbroken gauge symmetry generated by Ao- All the other modes A^. are 
massive and obtain their longitudinal modes from "eating" a A| field in a process 
very similar to the Goldstone mechanism in spontaneous symmetry breaking. Due 
to the mass terms, the gauge symmetry generated by the A&, k / is hidden in a 
simultaneous transformation of the A? and the "Goldstone" fields A\, which also 
is reminiscent of spontaneous symmetry breaking. Gauging away the A\ corre- 
sponds to unitarity gauge in a spontaneously broken gauge theory. 

Let's have a look at the Lagrangian to see whether this assertion is correct. The 
5D gauge kinetic Lagrangian is given by 

£gauge = —^F ab F a (3.10) 

with the 5D field strength tensor 

pab = Q a A b _ Qb A a 

The Lagrangian (13.10b can be split into 4D part C^d and a 5D part C^d 

W = —F^F^ + \ (d y A» + d»A 5 ) 2 



£-5D 

Performing the KK expansion and defining the KK masses as 

k 

mk = 2R 

we obtain the effective 4D Lagrangian 

^ oo oo 1 

w = — A E F k Vp ^ + E 9 (~ m k A k + d " A tf (3-n) 



2 



C-kA 



The £*;,4D are the kinetic Lagrangians for the 4D gauge fields, while the Ck,5D are 
easily identified as gauge invariant Stueckelberg Lagrangians. The Stueckelberg 
action is a well-known trick for giving mass to U(l) gauge bosons by spontaneous 
symmetry breaking, introducing a scalar, massless Goldstone boson. Therefore, 
the above interpretation is quite correct; in the 4D effective theory, the masses of 
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the A?, k > 1 fields arise from a mechanism equivalent to spontaneous symmetry 
breaking, yielding an infinite tower of massless Goldstone fields. 

In addition to the mass terms for A^, k > 1 and the kinetic terms for A\, the 
Lagrangian Ck.5D also contains a bilinear term which mixes A^ and d^A\. This 
can be removed by adding a suitable gauge fixing term to (|3.10|) similar to the 
gauges in spontaneous broken gauge theories 

C GF = 1 - : (dpA» - HdyA*) 2 



If we impose Dirichlet boundary conditions on A^, we find a result only slightly 
different. This boundary condition implies Neumann boundary conditions for the 
4D scalar A 5 . The massless Aq vanishes from the spectrum, and at the same time, 
a physical scalar Aq arises in which is not "eaten" by the gauge bosons. Impos- 
ing mixed boundary conditions would forbid the massless Aq mode as well as the 
physical scalar A^. 

The analysis of nonabelian 5D gauge theories proceeds along the same lines as 
the abelian case but is more complicated due to the more complex group struc- 
ture. This discussion can be found e.g. in [26 1 or [23 1. The most exciting difference 
when compared to the abelian case is the possibility of assigning different bound- 
ary conditions to different component gauge fields. As an example, consider a 5D 
SU(2) gauge field A a = A^r r . If we assign Dirichlet boundary conditions to the 
components belonging to the t\ and ti generators and Neumann ones to the T3 
component, then only the A\ tower will contain a massless mode mode in the KK 
decomposition. Thus, only the U(l) subgroup generated by T3 will make it into the 
4D effective theory as an unbroken gauge symmetry; all other gauge symmetries 
of the KK fields get spontaneously broken. 

Therefore, in a 5D version of the Standard Model, suitable boundary conditions 
could be used to break the SU(2)^ x U(l)y gauge symmetry down to the U(1)em 
subgroup without introducing a Higgs field. The price to pay would be towers of 
massless KK partners for all Standard Model particles. However, such a scenario 
is excluded by electro weak precision observables. Still, if we could cut off the KK 
towers after the first mode, this theory could be a starting point for the model sug- 
gested in chapter 1231 

Unfortunately, another feature of the nonabelian result is a more complicated 
transformation behavior of the KK modes under gauge transformation involving 
mixing between different modes [26|. In other words, we have to take all KK modes 
into account to retain gauge invariance. Simply cutting off the tower wouldn't do 
any good but instead destroy gauge invariance! However, there are loopholes: for 
example, dimensional deconstruction as discussed in section 13721 provides a gauge 
invariant way of cutting off the KK towers. 
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Fermions 

To describe 5D fermions, we have to find a five dimensional representation of the 
Clifford algebra, i.e. a set of five matrices 7° obeying the anticommutation relations 



{ 7 a ,7 fe } = 2 5 ab (3.12) 



Fortunately, this search does not lead far from well-known territory: if we take the 
usual set of 4D gamma matrices 7 M , then 7 5 has all the properties we need. Defining 



4 

7 as 



7 4 = i 7 5 



we obtain a set of matrices that obey (13. 12b . As a result, 5D fermions can be de- 
scribed by 4 component complex spinors just as the 4D ones. 

However, writing down the generators for 5D Lorentz transformations 



7°, 7' 



we note that the cj 5m don't commute with the 4D chirality projectors II-i- (c.f. dA.31) ). 
This means that 5D Lorentz transformations mix left- and right-handed fields, de- 
feating the notion of chirality: there is no such thing like a chiral fermion in five 
dimensions. Nevertheless, as 5D Lorentz invariance is broken down to the 4D sub- 
group anyway by the boundaries of the interval, it is possible to obtain a 4D ef- 
fective theory that makes a distinction between the chiralities by choosing suitable 
chiral boundary conditions. 

As an example consider a 5D fermion \f with left- and right-handed compo- 
nents ^ l/r = II-i-V]/. The massless 5D kinetic Lagrangian reads 

£^ = i^ a d a ^ = - ^7 5 ^ = 

i {^Lf^L + *) - ^idy^R + V R d y y L (3.13) 

If we impose Neumann boundary conditions on the left-handed field ^ l and Dirich- 
let ones on the right-handed field ^ r and insert the corresponding KK expansions 
into H3.13I) , we obtain the effective 4D Lagrangian as 

~— 00 /— — k — — 

= *0,L^0,L + Yl [ ^k,Lf^k,L + ^k,Rf^k,R - {^k,L^k,R + ^k,R^k,L) 
k=l ^ 

(3-14) 

From (|3.14[) it is then clear that the 4D effective theory contains a left-handed mass- 
less fermion ^o,L together with a tower of massive Dirac fermions. 

Other types of boundary conditions as well as an explicit 5D mass term change 
this picture and may result in more complicated situations. In particular, a mixing 
between different KK modes is possible requiring further diagonalization to obtain 
the mass eigenstates. However, as this is of no importance for this work, we refrain 
from delving deeper into these issues and stop the discussion at this point. Some 
further remarks on fermions in compactified 5D theories can be found e.g. in |27| . 
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Figure 3.1: Sketch of the deconstruction of a fifth dimension compactified on an 
interval [0; 2Rtt] to a lattice of N + 1 points. 



3.2 Dimensional Deconstruction 



Dimensional deconstruction is the reduction of one or more extra dimensions to a 
finite point lattice. This procedure is well-known from lattice gauge theory where 
Minkowski space is approximated by a four dimensional finite point lattice which 
allows for the numerical calculation of correlators by solving the path integral via 
Monte Carlo integration. 

In the context of model building, the deconstruction of an extra dimension 
opens up a new class of models which are akin to extra dimensional models but 
differ in some important aspects. In particular, deconstruction allows to cut off 
the KK towers in a gauge invariant fashion and provides a class of possible UV 
completions. 

Let us consider again the example of a 5D theory compactified on an interval 
which we discussed in the last section 

(x, y) G R A x [0; 2i?vr] 

and replace the fifth dimension by a discrete lattice of TV + 1 equidistantly spaced 
sites yo,...,yN 

, 2R7T 

y n = nd , d= 

This setup is sketched in fig. 13.11 



Gauge theory revisited 

In order to find out how a 5D gauge field A a is properly represented in the decon- 
structed theory, let us briefly recall some basics about gauge theories. Such a theory 
is invariant under a group of spacetime dependent transformations^ 

4>(x) — > u(x)4>(x) 

6 For the rest of this section, an unitary representation will be assumed together with the canonical 
trace normalization of the Lie algebra generators T t 
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of the fields <p. 

This invariance gives rise to the necessity of some prescription for how to com- 
pare fields at different points in space time in a gauge invariant way. Even in a local 
field theory, this is necessary because the definition of the partial derivative 

d^{x) = lim + 

e-K) £ 

(with the unit vector ra M pointing in direction fi) involves the field at two infinitesi- 
mally separated points and therefore is not gauge covariant 

{U(x)4>(x)) + U(x)d>*<p{x) 

Such a prescription is given by Wilson lines W(x\, X2), operators valued in the 
gauge group which transform as 

W(xi,x 2 ) — >■ U{x l )W{x l ,x 2 )U{x 2 ) ] (3.15) 

The gauge field is valued in the Lie algebra of the group and generates Wilson 
lines for infinitesimal displacements e M 

W(x, x + e) = I - ige^ix) + O (e 2 ) (3.16) 

with an arbitrary nonzero number g which ends up as the gauge coupling in the 
covariant derivative. Different finite Wilson lines W(x\,x 2 ) can be generated by 
iterating (|3.16|) along paths connecting x\ and x 2 (see e.g. |5J for details). (|3.16|) 
implies that the gauge field transforms as 

A"(x) — ► U(x) (V(x) + ±dA U(x)t 

Due to their transformation behavior, the Wilson lines can be used to connect 
different points on the manifold and define a covariant derivative: 

D , m = ]im W{x,x + en^{x + en»)-^x) 
t->o e 

Inserting H3.16I) into (13.17b we retrieve the well known expression for the covariant 
derivative 

D^(j)(x) = d^ftx) - igA^<j>(x) 

which by construction is gauge covariant as desired. 

The transformation law (|3.15b also implies that Wilson lines can be plugged 
together to form new Wilson lines 

W(x 1 ,x 2 )W(x 2 ,x 3 ) = W(x 1 ,x 3 ) 

In particular, it follows that the closed Wilson line defined as 

K{x) = W(x, x + e)W(x + e, x + e + rj)W(x + e + 7],x + i])W(x + 77, x) (3.18) 
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with two small displacements and i]^ transforms as 

K(x) — > U{x)K(x)U{x) ] (3.19) 
Plugging (|3.16|l and the Taylor expansion of A^ 

A»(x + e) = A"(x) + e^A^x) + O (e 2 ) 
into the definition of K (x) (13.181) , we obtain 

K(x) = I - ige^ (0M"(a;) - 9^) - ^ [A^(x), A"(x)]) +0 (r? 2 , /i 2 ) 

recovering the well known expression for the field strength tensor F^ LU . The trans- 
formation behavior \3.\9) must hold in every order in e M and rf , and it follows 
that 

F^(x) — »• £/(x)i^(x)£/(x) t 

Therefore, we can use to write down a gauge invariant Lagrangian that in- 
volves bilinears in the derivatives of A M and which therefore is suitable to give 
dynamics to the gauge field 

£gauge = ~2^~ r Fjav 

Summing up the discussion: the gauge field acts as connection that relates 
gauge transformations at adjacent points on the manifold to each other and there- 
fore allows to compare the value of fields at different points in a gauge covariant 
way. This is strongly reminiscent of Riemannian geometry where the Christoffel 
symbols provide such a connection w.r.t. coordinate transformations. The connec- 
tion can be integrated to provide a parallel transport along finite paths, much like 
the Wilson lines for gauge transformations. Indeed, there is a geometrical descrip- 
tion of gauge theory in terms of fiber bundles [28 J. The field strength tensor is 
built in a similar way as the curvature tensor by considering the transport along an 
infinitesimal parallelogram. 

Deconstructing gauge theory 

Let us now carry the above description of gauge invariance over to the decon- 
structed 5D theory. Going from the continuous coordinate y to discrete sites y n , 
we've got a copy of the 4D gauge group Q n at every lattice site, and the theory is 
invariant under transformations of the fields 

(j)(x,y n ) = (p n (x) ->■ U n (x)(f) n (x) (3.20) 

If we compare the values of one of the 4> n at two different 4D coordinates <j) n {x\) 
and <j> n {x2), then we can still use 4D gauge fields An(x) localized at the lattice sites 
to generate Wilson lines W n {x\,X2) which compensate for the change in gauge. 
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Go 1 Qi J // Gn-M Gn I 

E 2 (x) Sjv-i(x)^^r ^0*0 
<(*) ^^(x) A*(x) 

Figure 3.2: Deconstructing a gauge theory. At each lattice site, there is a 4D copy 
Gn of the gauge group G with the gauge field An, the sites being linked with Wilson 
line fields E n . 

However, as far as the fifth coordinate is concerned, there is no such thing as an 
infinitesimal displacement in y anymore after deconstruction. Therefore, the no- 
tion of a field A n (x) which connects points infinitesimally separated in y doesn't 
make any sense at all. The smallest displacement in y we can manage in the de- 
constructed theory is a hop between adjacent lattice sites. Since this is a finite dis- 
placement, we need finite Wilson lines T, n (x) which connect the sites and which 
transform as 

E n (x) — > U n -i(x)Z n (x)U n (xy (3.21) 

Therefore, in the deconstructed theory, we must replace the 5D gauge field 
A a (x, y) with N + 1 4D gauge fields A£(x), A" N {x) and N Wilson line fields 
Ex (a;), . . . , Ejy. This setup is shown in fig. 13.21 From these objects, we can build 
up Wilson lines that connect arbitrary points on the lattice (x\,y ni ) and (x2,2M 2 )- 
The E n are valued in the gauge group and can be parameterized using the group 
generators 

S n (x) = exp ( i^^n) (3.22) 



The scale v is introduced to make the component fields a n k dimensionful quantities 
and will later be fixed by matching the mass spectrum of the deconstructed theory 
to that of the continuous one. As the E n are valued in the gauge group, they must 
have a nonvanishing vacuum expectation value. If the component fields a n k don't 
develop nontrivial vacuum expectation values, then (|3.22[) tells us 

(E n ) = I (3.23) 

A kinetic term for the E n can be obtained by finding the analog to F fl5 (x, y) in 
the continuous case. According to our previous discussion, F^ 5 is obtained from 
the Wilson line along a parallelogram spanning an infinitesimal distance both in a 
4D direction and in y. As we can only go finite distances d in the discrete theory, a 
deconstructed analog is the quantity 

K n {x) = Z n (x)Wn(x,x + e)Z n (x + e)^W n -i(x + e,x) (3.24) 

which corresponds to an infinitesimal displacement in 4D direction and a finite 
hop from y n -\ to y n in the discrete extra dimension (c.f. fig. 13.311 . Expanding (|3.24l) 
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. £ n (x + e)t -V 

f^M< Vn J 



x o,i,2,3 y /W n {x,x + e) 

X 



x + e ,„."X S n (x + e) t : 7^ x + e 
W n -i(x + e,x) 

Mn-M ► ■ Vn) 

S n (x) 

Figure 3.3: Definition of the quantity K n as the Wilson line along a parallelogram 
spanning two adjacent lattice sites and an infinitesimal 4D distance. 

in a similar fashion as (|3.18[) , we obtain 

K n (x) = I - C/1 (0"E„Oc) - ig (A^x^nix) - Z n (x)A%(x))) S n (x)t +0 (e 2 ) 

* v ' 

By construction, F„ must transform as 

F£(x) -> C/ n (x)F^(x)[/ n (x) t 
and therefore, we can write down a gauge invariant kinetic term for the S n 

AT 
71=1 

£s is chosen such that the component fields a n k are canonically normalized. 

As the covariant derivative of the S n is fixed by their transformation behavior 

D^ n = d^ n - ig «_!S n - S n <) 

it follows that 



and finally, we can rewrite the complete Lagrangian for the deconstructed theory 
as 

N N 

£dec = " g E ^ F r^, + « 2 ^ Tr (C£ n ) ( J D M S n ) t (3.25) 

n=0 n=l 

If we expand the S n fields in their component fields cr n/ t according to 03.22D we 
find that the lowest order (aka the vacuum expectation value of the S n ) generates 
mass terms for the component gauge fields A^ fe : 

2 2 N y 
£mass = — Y~ ^2 (^n-l,k + ^n,k ~ = ~^^ nmJ Ki,k^m,k^ 

n=l ~~ 
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Figure 3.4: Deconstructing a gauge theory with Dirichlet boundary conditions. The 
condition A fl (x, y)| yg { 2R-k) = removes the gauge groups and gauge fields at the 
yo and y^ sites of the lattice (c.f. fig. 13.26 . 



with the real, symmetric mass matrix 

/I 



M 



2 2 

v g 



■1 











1 



0\ 





-1 

1 J 



The eigenvalues of this matrix (see e.g. |29| for the eigenvectors) give the gauge 
boson masses 



4v 2 g 2 sin z 



nir 



n 



0,. 



(3.26) 



2(iV + l) 

For every component gauge field, the spectrum contains a massless mode together 
with iV massive ones. As we have N sigma fields, we additionally get N scalar 
fields for every group generator which is just the amount of Goldstone bosons that 
are required to constitute longitudinal modes for the massive vector bosons. 
In the limit of large N we can approximate (I3.26D by 



vg 



nir 
N + l 



Therefore, if we fix the scale v as 



N + l ^ 1 
2gRir gd 



we recover the mass spectrum of the continuous case with Neumann boundary 
conditions from the deconstructed theory in the large iV limit. It seems that our 
way of deconstructing the theory has been the correct prescription for a 5D gauge 
field compactified with Neumann boundary conditions. 

What about Dirichlet or mixed boundary conditions on A^(x, y)l Setting one or 
all components to zero at both ends of the interval corresponds to setting the gauge 
fields at yo and / or ypj to zero and basically kills gauge invariance at these sites 
(c.f. fig. 13.4b . However, we have no way of transferring the appropriate boundary 
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conditions on A 5 (x, y) to the deconstructed version of the theory and can only hope 
that the theory takes care of this itself. 

This is indeed the case as can been seen by counting the scalar and vector fields 
involved: if we remove the Aq and fields, we obtain N — 1 massive gauge fields 
for every generator, while the number of scalars remains N. After all the Goldstone 
bosons have been eaten, a physical scalar remains for every generator, which we 
have argued in the last section also happens in the continuous case with Dirichlet 
boundary conditions. In the mixed case, we have N massive vectors and N scalars 
which all get eaten for every generator, also matching the continuous case nicely. 

Although deconstruction as presented here does a rather nice job at approximat- 
ing 5D gauge theories, the formalism had first been proposed in [30J as a tool for 
constructing renormalizable UV completions to 5D theories. At low scales, the de- 
constructed theory virtually cannot be distinguished from a true 5D theory. How- 
ever, while such a theory is inherently nonrenormalizable and must break down 
at some cutoff scale, the deconstructed version can be made renormalizable rather 
easily. 

For example, note that the S n are similar to the £ field introduced in chapter l2~3l 
to describe the spontaneous breaking of electroweak gauge symmetry. Therefore, 
the deconstructed theory can be UV completed just like the Standard Model by 
introducing Higgs fields H n for each of the S n fields. 

A different, technicolor-like approach to an UV completion is to describe the 
S n fields similar to the Pions in QCD as the excitation of a condensate of strongly 
interacting fermions. At energies above the confinement scale of this new strong 
interaction, the S n fields would be replaced by the "quarks" and "gluons" of the 
new strongly interacting sector, a theory which is known to be renormalizable. 

For our purposes, it is important to note that deconstructing the theory retains 
gauge invariance. Therefore, deconstruction presents a consistent gauge invariant 
way for cutting off the infinite KK towers that arise in the continuous case. 



Matter fields 

The correct prescription for the discretization of matter fields has already been 
given in l|3.20|l : the field <p(x, y) is replaced by a set of iV + 1 fields 4> n (x). The Wil- 
son line fields S n then allow for a straightforward discretization of the covariant 
derivative^ D 5 

dy<p(x, y n ) — > D n = 

The discretized covariant derivative can be used to build a discretized version of 
the Lagrangian for the matter fields, e.g. for a real scalar 

N N 

C * = 2 £ {D^nD^n - m 2 cg) DlD n (3.27) 

n=0 n=l 



7 In fact, different possible discretization prescriptions can be devised, each leading to a slightly 
different mass spectrum, but all having the same high energy limit. 
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Figure 3.5: The structure of the Three-Site Higgsless Model in moose notation. See 
the text for explanation. 

Expanding the second term in (I3.27D we obtain 




We find that the discretized covariant derivative terms yield contributions to the 
diagonal mass terms of the 4> n as well as Yukawa-type "hopping" terms that cou- 
ple fields at adjacent lattice sites and contribute off-diagonal elements to the mass 
matrix. Diagonalization of the mass matrix then leads to a finite tower of massive 
modes which matches the infinite KK tower in the N — > oo limit. 

If the discretization procedure is applied to fermions, a new issue arises. Plug- 
ging in the discretized derivative and calculating the mass spectrum, an unphysi- 
cal doubling of modes is observed: for every mode in the continuous theory two 
modes arise in the lattice approximation, one of which is unphysical. This problem 
of fermion doubling as well as the solution by adding a so-called Wilson term to the 
Lagrangian is well known from lattice gauge theory. A detailed discussion of the 
issue in the context of deconstructed QED can be found in (311 . In the Three-Site 
Model, however, this problem does not arise due to the extremely small size of the 
lattice, so we won't elaborate on this point. Instead, the interested reader is referred 
to above reference. 

3.3 Constructing the Lagrangian 

With the introduction of 5D gauge theories and dimensional deconstruction in the 
last two sections, we now have the tools at hand for the actual construction of a 
Lagrangian implementing the Three-Site Higgsless Model [4 J sketched in chapter 



CHAPTER 3. THE MODEL 



— TOP-DOWN APPROACH 



37 



In II23B it is shown that in 5D gauge theories the unitarity cutoff is pushed to 
the NDA cutoff scale by the exchange of KK gauge bosons. Therefore, a 5D gauge 
theory is a good starting point for a theory in which unitarity violation is delayed 
by heavy W' and Z' partners. Dimensional deconstruction can then be used to cut 
off the higher KK modes in gauge invariant way retaining only the Standard Model 
particles and one generation of KK partners. 

Fig. 13.51 shows the structure of the Three Site Model in moose notation [32|. 
The gauge group consists of two SU(2) factors and one U(l) factor which can be 
understood as a 5D SU(2) gauge theory broken by boundary conditions to U(l) 
on the right brane and deconstructed to three lattice sites. All three group factors 
are given their own gauge couplings g, g and g 1 which is required in order to allow 
the model to be consistent with experimental data and which corresponds to a y 
dependent gauge coupling g(y) in the continuous version of the theory. QCD is 
left unchanged in the Three-Site Model with only one set of gluon fields in the 
spectrum. The gauge kinetic Lagrangian is 

1 1 1 1 

£gauge = - 2 E ^ F n V Fn„u - F^ v - -Tr GTG^ (3.28) 

n=0 

It should be noted that this setup is identical to the BESS model ||33| . 

Deconstructing the gauge group gives us two nonlinear sigma type Wilson line 
fields Si and S 2 which connect adjacent lattice sites. Under gauge transformations 
generated by parameter fields Ao = AofcTfc, Ai = Xikfk and A 2/ the Wilson lines 
transform as 

Si — > e iAo S ie ^ Al , S 2 —> e ai S 2 e- iA2T3 (3.29) 

The transformation behavior (|3.29[) fixes the covariant derivative of the S fields 

i^Si = ^Si - igA^x + igViA* 
D"S 2 = d^ 2 - igA^ 2 + ig'A^ 2 T 3 

With the covariant derivatives fixed, we can write down the kinetic term for the 
S1/2 fields 

2 

£x = v 2 ^Tr {D^ n f (D M S n ) (3.31) 

71=1 

Expanding (|3.31|) in terms of the component fields of E]/ 2 gives mass terms for the 
gauge bosons. The mass scale is set by v and fixed by the W and Z masses. 

The fermion sector is inspired by deconstructing 5D fermions with chiral bound- 
ary conditions. The Three-Site Model comes with two copies of each Standard 
Model fermion, the chiral components of which are distributed according to fig. 13.51 
Each Standard Model isospin doublet gives rise to left-handed doublets under the 
respective SU(2) at sites and 1, a right-handed SU(2)i doublet at site 1 and two 
right-handed fermions at site 2. As shown in tab. 13.1 1 the fermions at site 2 inherit 
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I'OL / &1L / *lJi 


^2R 


Neutrinos 


1 

2 





Leptons 


1 

2 


-1 


Up type Quarks 


1 

6 


2 
3 


Down type Quarks 


1 

6 


1 

3 



Table 3.1: U(l)2 charge assignments for the fermions on the three site moose. 



their U(l)2 charge from the electromagnetic charge of the corresponding Standard 



Model fields, while those at the other two sites inherit the hyperchargqj. 

According to the above charge assignments, the covariant derivatives of the 
fermions (omitting the gluon contribution in the case of quarks) read 

D^ 0L = (d» - igA% - ig'A^Y) V 0L , D^ lL = (d* - i~gA» - ig'A^Y) * 1L 
D^ lR = [&> - igA$ - ig'A%Y) V 1R , D»*% = [d» - ig'A^Q) 

with the charge operator Q and the hypercharge operator Y. The fermion La- 
grangian then read^l 

1 2 
n=0 n=l 

+ v (x^ol^ir + X^i^ir + ^i L S 2 ^2R + h.c.J (3.32) 

In (|3.32|) , the Yukawa couplings between fermions at adjacent lattice sites are left 
as free parameters, an effect which can be achieved in the continuous theory by 
the introduction of brane kinetic terms. In particular, the Yukawa couplings X' u > d 

between vl/u and the ^f R ^ a U° w to accommodate the different masses of up- and 
down-type fermions observed in nature. The X' u / d are the only Yukawa couplings 
that have a nontrivial flavor structure, all other couplings are taken as flavor uni- 
versal. 

With the definitions 03.28D , (I3.31D and (I3.32D , the complete Three-Site Lagrangian 
reads 

^3-site = <£gauge + £s + £/ (3.33) 



8 While being a slight deviation from straightforward deconstruction as it corresponds to a nonlo- 
cal coupling of the fermion fields to the gauge field on the right brane in the continuous theory, this 
setup is necessary in order to reproduce the correct hypercharges of the Standard Model fermions 
without introducing additional neutral vectors. 

9 In all expressions containing fermions, an implicit sum over all fermion flavors present in the 
Standard Model is assumed. 

10 Flavor mixing is not included in the original model, but CKM-type mixing can be easily achieved 
through the \' u / d , retaining the flavor-universal nature of the other couplings 0). 
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How does this model relate to that sketched in chapter I2.5P In the limit of large g 
and large A we can integrate out the fields at the bulk site in order to obtain a "two 
site model" with one linear sigma field 

S = S1S2 

Comparison to chapter 12.31 identifies this two site model as the ordinary Standard 
Model with SU(2) playing the role of SU(2) L and U(l) 2 that of U(l)y! The bulk 
lattice site adds the additional SU(2) group factor and the set of fermion partners 
suggested in chapter 1231 and the desired mixings arise from the terms that connect 
the lattice sites via the fields, corresponding to the 5D kinetic terms in the 
deconstructed setting. 

In accordance with JH, we introduce the following definitions for later conve- 
nience 



9 

x = — 

a 

A 
A 



9' 

y = — 
9 

. A' 



t = tan 6 



ff_ 
9 



cos 9 



(3.34) 



A 



As argued above (and discussed in more detail in the next chapter), the Standard 
Model is recovered in the limit of large g and A, and therefore, we should expect 
that x, y, ei and e ^ are small quantities if the new structure is indeed a small per- 
turbation on top of the Standard Model. 
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Chapter 4 

Model Properties 



Who ordered that? 



(I. I. Rabi on the discovery of the muon on 1936) 



At the end of the last chapter, the Lagrangian of the Three-Site Model as presented 
in 01 was assembled. The present chapter is devoted to a detailed discussion of 
the model. In the first section, the mass spectrum is calculated. In the second sec- 
tion, the free parameters present in the model in addition to the Standard Model 
parameters are identified. A number of experimental constraints together with the 
resulting parameter space is discussed. With the third section follows a brief dis- 
cussion of the couplings between the Standard Model particles and their new heavy 
partners, and the last section is devoted to the widths of the new heavy particles. 
An autogenerated sample spectrum is printed in appendix [0 

4.1 Masses 

The first step in studying the phenomenology of the Three-Site Model is the calcu- 
lation of the mass spectrum. The mass terms for the gauge bosons can be readily 
extracted from (13.311) 




i 

M CC W + ^W~ + -M"B?B~ 



(4.1) 



with the charged gauge bosons 



W ±fl 



V2 



(n = 0,l) 



the neutral gauge bosons 



B£ = Ah (n = 0,l) 



B$ = A* 
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and the mass matrices 

M cc = v 2 ~g 2 (^ x , M NC = v 2 ~g 2 \-x 2 -tx ] (4.2) 

where we have used the definitions (13.34b . 

The fermion mass terms follow from the Yukawa couplings in (|3.32|l 

£mass,fermions = (HlMl) Mf ^ ^ ^ (4.3) 

(an implicit sum over all fermions present in the Standard Model is assumed) with 
the fermion mass matrices 

'e L 



M f = v\^ e ,j (4.4) 
Integrating out the bulk: the two-side model 

To get some intuition about the connection between the Standard Model and the 
Three-Site Model, let's pick up again the argument presented at the end of chapter 
13.31 and integrate out the bulk lattice site in the limit 

g — > oo , A — > oo (4.5) 

where we keep the other parameters (in particular A and A') finite in the limiting 
process such that the limit implies 

Q Q 1 A , A'r 

x = ^->0 , y = ^^0 , e L = T ^0 , e' f = -J- -»■ 
9 9 A 1 A 

for the quantities defined in (|3.34|) . 

The masses of the bulk gauge bosons and of the bulk fermions that are to be 
integrated out can be read off from (|4.1[) and (|4.3|) 

M Al = V2gv , M bulk = v\ (4.6) 

The resulting model lives on only two lattice sites, connected by one nonlinear 
sigma field 

X = S1S2 

To stress the connection to the Standard Model, we will use the same identifiers as 
in chapter |2~31 for the two-side model fields 

The covariant derivative of S is then 

D^T, = <9 M £ - igW^T, + ig'B^T 3 (4.7) 
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and the kinetic Lagrangian for £ is 

£e,2S = v 2 Tr (D M E)t (£)M S ) (4.8) 

As we want to reproduce the Standard Model in the limit (|4.15|) , we can express v 
through mw and g by virtue of (|2.13|) 

mw 
v = 

9 

To find out how v relates to v, let's perform the matching to the Three-Site Model 
and integrate out the bulk gauge field from the diagram 




The vertices can be read off (|4.1|) , and the resulting change in the Lagrangian is 

A£ m , w = H) {-iv 2 g~g) {-iv 2 g~g) W + »W~ = -l v 2 g 2 W + »W- 



modifying the mass term for W. 



± 
i-' 

C m ,w = v 2 g 2 W+>+W~ + AC rn , w = \ 2 g 2 W^W- 

Comparing this to (|4.7|) and (|4.8|) identifies v as 

v 

V = —= 
V2 



(4.10) 



The factor \/2 has an interesting effect: plugging v into the formula for the NDA 
cutoff (|2.18[) in order to estimate the UV cutoff of the Three-Site Model, we obtain 

A N da = 8ttv = 8V2ird = 8\/27r— « 4.4 TeV 

9 

Comparison to (12.1811 shows that the NDA cutoff scale is elevated by a factor \/2 in 
the Three-Site Model above that of the Standard Model. 

In fermion sector, Yukawa couplings between the left-handed fermions sitting 
at site and the right-handed ones at site 1 (formerly site 2) are obtained by inte- 
grating out the bulk fermions from diagrams of the type 



(4.11) 
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The vertex factors can be read off from (14.41) and inserted to obtain the effective 
Yukawa couplings 

£Yukawa,2S = v (W f L e L \'V f R + h.c.) = M bulk (¥ f L e L f? f * f R + h.c.) (4.12) 

with an implicit sum over all Standard Model fermions. We can read off the fermion 
masses from (14.12b 

rrif = M hu i k e L e' f 

After choosing the flavor-universal ex, the are to be chosen for each Standard 
Model fermion separately to yield the correct Yukawa couplings and masses. 

If we would drop the assumption of taking the limit A — > oo such that e* — > 
and instead keep the e'f finite (which is perfectly possible by letting ei go to zero 
fast enough), diagrams like 



(4.13) 




R *R 
would induce finite wave function renormalizationsL] 



— >• Jl + e' 2 M 



R 7 V / R 

In addition, diagrams of the type 




(4.14) 



W/B 



would induce finite couplings e.g. of two \E ' R to a W , leading to deviations from 
the Standard Model even in the limit g — > oo , X — > oo! Similar arguments apply if 

1 Due to the chirality of ty' R , the momentum expansion of the diagram <4.13t starts at the order f, 
leading to the change in the Lagrangian 

AC = ie' s 2 ¥ R $¥ R 
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€l is kept finite, and therefore, the correct limit for recovering the Standard Model 
is in fact 

x^O , M bulk ^oo , e L ^0 , e^O (4.15) 

while keeping g, g' and v finite! 

In the limit (|4.15l) , we fully recover the Standard Model and obtain the matching 
conditions (|4.10[) and (|4.12|) . By comparison to chapter 12.31 we can immediately 
write down the spectrum to lowest order in an expansion in x, el and e'j : 



The massless photon and gluons. 

The Standard Model and Z bosons with masses 



vg \g 2 + g' 2 

• The Standard Model fermions with masses given by 

mf = M hu i k e L e' f 

• Two W^' and one Z' which are degenerate with mass 

mw' = fnz 1 = V%gv 

• One partner fermion for each Standard Model fermion with mass 

M bulk = Xv (4.16) 

Exact spectrum 

To get the corrections for finite x, €l and e' f to the above picture, we must explicitly 
diagonalize the mass matrices (|4.2|) and (|4.4[) . This was done analytically with the 
results then being expanded as a series in x resp. in order to obtain the expres- 
sions for masses and wavefunctions presented in this section. 

As the gauge sector of the Three-Site Model corresponds to a deconstructed 
5D SU(2) gauge group with Neumann boundary conditions for the T3 component 
field and mixed ones for the other two, an unbroken U(l) gauge symmetry together 
with a massless gauge boson is to be expected. Indeed, setting 

Ao = Ai = A2T3 

in l|3.29|) , we find that the vacuum expectation value of the T>i/ 2 is not changed 
under the transformation and, therefore, this combination generates an unbroken 



2 Indeed, while being more intuitive considering the parameterization of the model we will be 
developing in the remainder of this chapter, (14.15ft together with the requirement of keeping g, g' and 
v finite is just a rephrasing of the limiting conditions stated already at the beginning of this section. 
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m 



fW 

Jo 



l-^-0(x 4 



/r = § +^-o(x 5 ) 



m^ = 2ffV(l + ^ + 0(x 4 



1 



l-^-0(x 4 



f z 
Jo 



x 2 c 3 (l+2t 2 -3t 4 ) 



+ 0(x i ) 



rrir 



ac(l-i 2 ) x 3 c 3 (l-t 2 )~ 



16 



+ £> (a; 5 ) 



2 SC 2 (3-2t 2 -t 4 ) 
8 



+ O (x 4 ) 



JO 



*(l-3t 2 ) 
16 



+ O (x 5 ) 



m' 



v = 2yV (l + |J + (x 4 



/f 



x 2 (l+i 2 ) 



+ O (x 4 ) 



/I 



xt 



+ 



16 



+ O (x 5 ) 



Table 4.1: Masses and wavefunctions of the gauge bosons as obtained from diago- 
nalizing the mass matrices M and M NC (c.f. (|3.34|) for the abbreviations). 



U(l) gauge symmetry. Looking at the fermion charges tab. 13.11 this is readily iden- 
tified as the electromagnetic U(l) em , and defining the electromagnetic gauge cou- 
pling e 

4 = 4 + 4 + "72 ( 4 - 17 ) 

e 2 ^2 ^2 ^2 

the corresponding photon field is 

A» = -Bft + -B* + -,B» (4.18) 
9 9 9 

The masses and wavefunction^l of the massive gauge bosons are given in tab. 14.11 
as an expansion in x. As expected, the masses coincide at lowest order with the 
above result obtained by integrating out the bulk site. The most noticeable effect 
of the corrections is to lift the degeneracy between the W and the Z'. Turning our 
attention to the wavefunctions, we find that the Standard Model W and Z bosons 
are localized at the brane lattice sites with f\ vanishing in the limit x — > 0, while 
the W ±r and Z' are localized at the bulk site with /o and /2 being of order O (x). 



' The wavefunctions are to be read as e.g. 
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M hu i k e L e' f 



rrij 



1 + e' 



\ 



2 ( e' t 2 + 1 



+ 0(4) 



/ 



f f 
JQL 



-1 + 



+ o (4) 



'1L 



1+.'/ 2 (e' 2 +l)<' 



+ o (e|) 



1+e' 



+ 



l 2\ $ 



(l+e/) 



0(4) 



■/2.R 



'1+e', 



+ 



e 2 e' 2 



(1+4 2 ) 



+ o (4) 



«V = ^4uik 

/ 

1 + 



1 + e' 



\ 



2 e'/ + 1 



+ 0(4) 



/ 




1R 



e' 2 e 2 



/1+e/ (l +e /) 5 



0(4) 



J2R 



4-1 + ^^ + 0(4) 

^ (1+4 2 )" 



Table 4.2: Fermion masses and wavefunctions expanded in el as obtained from 
diagonalizing the fermion mass matrix (c.f. (|3.34|) for the abbreviations). 



The masses and wavefunctions resulting from diagonalizing the fermion mass 
matrix (14.411 are given in tab. 14.21 as an expansion in e^. Again, in the limit (|4.15l) 
the result from integrating out the bulk site is recovered, with the corrections lift- 
ing the degeneracy between the heavy partner fermions. As in the case of the 
gauge bosons, the Standard Model fermions are predominantly localized at the 
brane sites, while their heavy partners reside mainly in the bulk with the modes 
becoming fully localized in the limit (|4.15|> . 

In order to illustrate the formulae in tab. l4.1l and tab. 14.21 fig. l4.1l shows the wave- 
functions of W, Z and massless left- and right-handed fermions as well as those of 
their heavy KK partners, the values being taken from app. |B] These profiles are 
representative for the whole allowed region of parameters space and clearly ex- 
hibit the (de)localization properties of the different modes as discussed above. The 
wavefunctions of t /if and b/b' differ slightly from those shown in fig. 14.11 due the 
comparatively large mass of the respective KK light modes, the most notable dif- 
ference being a slight derealization of the right-handed components which are 
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X 



W 



Figure 4.1: The wavefunctions / n of W, Z and massless left- and right-handed 
fermions as well as those of their heavy KK partners for mw = 500 GeV, = 
3.5 TeV and ei = 0.237 (ideal derealization, see section 1421 below) , the values are 
taken from app.|Bj The scale is linear, and the dotted lines correspond to f n = ±1. 

(nearly) completely localized for the (almost) massless fermions and their KK part- 
ners. 



4.2 Parameter space and constraints 

While it indeed manages to delay the unitarity cutoff without invoking any scalar 
fields, the Three-Site Higgsless Model is not an improvement over the Standard 
Model as far as the number of free parameters is concerned. To fix the model and 
make predictions, all parameters present in the Standard Model (except for the 
Higgs mass, obviously) must be specified in addition to a couple of additional ones. 

In the gauge sector we've got four free parameters: the three gauge couplings 
g, g' and g as well as the scalar v. After fixing the W and Z masses as well as the 
electric charge, only one degree of freedom remains in this sector, which can be 
chosen as x. From the expansions in tab. I4.1l it follows that to leading order in x 



m 



w 



m 



w 



x 

T 



(4.19) 



for 



Therefore, we can fix x from the mass ratio of W and W, and the limit x 
recovering the Standard Model is equivalent to mw — > oo. 

Going over to the fermion sector we have two flavor-universal parameters 
and A and one set of parameters e'* that are to be chosen separately for each fermion 



4 It is possible to choose two different scales vi and V2 for the link fields Ei/ 2 - The choice vi = 
V2 = v maximizes the delay of unitarity violation |34|, however, this relation is not honored by loop 
corrections 1351. 
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flavor. The e'^ are fixed by the masses of the Standard Model fermions, while A and 
€l remain as free parameters. We choose to replace A by the heavy fermion mass 
scale Mb u ik (|4.16|l which has a more intuitive physical meaning. 



Electroweak precision observables at tree level: ideal derealization 

In chapter 12.51 the introduction of heavy fermions was suggested as a loophole 
for evading the bounds coming from precision measurements at LEP and LEP-II. 
In these experiments, the neutral- and charged current correlators were measured 
very precisely, revealing an astonishing degree of agreement with the higher order 
predictions of the Standard Model and putting tight bounds on any contributions 
from new physics. 

Under the assumption that all contributions of new physics manifest them- 
selves as modification to the gauge boson self energies, the resulting constraints can 
be summarized as bounds on a small number of parameters. The LEP bounds are 
commonly parameterized by the aS, aT and all parameters |l9l or equivalently 
by the three ei/2/3 parameters [21] which correspond to bounds on the coefficients 
in an expansion of the self -energies to order p 2 . 

The data obtained at LEP-II allows for additional constraints on new physics 
contributions to the self-energies, corresponding to an expansion to order p A . To- 
gether with the LEP-I constraints, these bounds can be expressed by a set of 7 pa- 
rameters introduced in II36II which also cover the bounds expressed by aS, aT and 
aU. 

While these parameters were designed to parameterize deviations from the 
Standard Model that appear only as "oblique corrections" in the gauge boson self- 
energies, the authors of II37II show that non-oblique corrections appearing in Hig- 
gsless models can also be absorbed into the extended set of parameters defined in 
H36J. 

An generic analysis of the precision observables in a large class of deconstructed 
Higgsless models which also contains the Three-Site Model was published in [38 J. 
There, the general form of the charged and neutral current scattering amplitudes 
in the vicinity of the Z and W poles was derived and it was shown that a suitable 
tuning of the fermion wavefunctions dubbed "ideal derealization" can minimize 
the tree level corrections to the precision observables. In particular, the Ap param- 
eter which measures the low energy difference between the charged and neutral 
current coupling strengths vanishes, and the tree level contributions to aS and aT 
reduce to 



aS = 4s 2 w c 2 w m 2 z ( -5- " -5- J , 

n \ m Z n m W n J 



aT 




1 



m w n 



(4.20) 

where the sums run over all KK gauge bosons and cw and sw are sine and cosine 
of the Weinberg angle 

a m w 
cw = cos vw = 

m z 
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For the case of the Three-Site Model, plugging the expansions of the masses tab. 14.11 
into l|4.20|l and using (|4.i9|) as well as 



raw 
m z 



+ 0(x) 



(4.21) 



gives the naturally small results 



aS 



m 



+ 0(x* 



aT 



m 



Wj 



+ o 



rrtr 



m 



m 



The condition for ideal derealization is the vanishing of the couplings of the 
light Standard Model fermions to the heavy W partners. In the Three-Site Model 
there is only one W , and the corresponding coupling is given by the overlap of the 
wavefunctions 

9wff - -j= [9 Jo Jol Jol + 9h Til Hl j 

As e'f vanishes for massless fermions and the other two parameters of the fermion 
sector Mb u ik and are flavor-universal, all massless fermions have the same wave- 
functions f*. Due the mass matrix being symmetric, the wavefunctions f w and 
f w are orthogonal and, therefore, choosing 

h 



OL 



f 



(X 




or equivalently 




x 



n 



w 



w 



(4.22) 



is sufficient to implement ideal derealization. Examining the wavefunctions given 
in tab. l4.1l and tab. !4.2l reveals that (|4.22|) encodes a relationship between and x 



4 = ^ + 0(x 2 ) 



(4.23) 



which reduces the number of free parameters by one, allowing us to fix the Three- 
Site Model by specifying just mw and in addition to the Standard Model 
parameters, (in the case of ideal derealization). Also, (|4.23|) automatically implies 
€l — > in the limit x — > as required by (|4.15|) in order to recover the Standard 
Model. 



Electroweak precision observables at one loop 

Due to the extremely efficient suppression of the tree level contributions to the 
precision observables in the case of ideal derealization, the loop contributions can 



CHAPTER 4. MODEL PROPERTIES 



51 



10 r 



w 1 max 

aT mav = 5-10 



2.5-10"* 

3 



max 

Ideal delocalization 



300 400 500 600 700 800 

W mass m w > [GeV] 

Figure 4.2: Lower bound on from estimating the t' / b' one loop contribution 

to aT in the ideally delocalized scenario. 

be expected to be dominant. Because the model contains heavy copies of the top 
and the bottom quarks, new contributions to aT can arise from diagrams of the 
type 

t' 
b' 

As 6l honors isospin, the isospin violating effect of this diagram persists in the limit 
el — > 0. The leading contribution to aT in this limit is estimated in EH as 

From tab. 14.21 we can approximate e' t as 



(4.25) 



Plugging this into (14.241) with an upper bound aT max on aT and using the condition 
for ideal delocalization (|4.23[) , we obtain a lower bound on Mb u ik 

1 m t 1 m t m W> ,,r^^ 

M bulk > Hr = o 4.26 

" 4vrV2aT max ve\ 87rV2aT max vm 2 w 

This constraint is shown in fig. !4.2l for two different values of aT max (the exact value 
depends on the "reference Higgs boson mass" chosen for the determination of aT, 
see HI). 
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Figure 4.3: Region in the Mt, u ik - gw'ff plane consistent with aS and aT at one 
loop for mw' = 500 GeV and different values of the UV cutoff A; gw is the isospin 
gauge coupling. Plot taken from [35 J. 

A nice consequence of (14.261) is the fact that it also induces an upper bound on 
the e' f 

1 ml 

Together with (14.2311 and (14.1911 , this implies that any limit my/' —> oo, M^uik — > oo 
which respects the precision constraints leads to e'f — > and therefore recovers the 
Standard Model according to d4.15lF . 

A complete one loop calculation of aS and aT was started in Il39ll40l and com- 
pleted in 1 35 1 . The result is quite interesting as, although a sizable piece of param- 
eter space remains, the case of ideal derealization is excluded. As an example, 
fig. 14.31 shows the region in the - gw'ff plane which is consistent with aS and 

aT at 95% confidence level for mw' = 500 GeV. The resulting gw'ff coupling is 
constrained to be smaller than some 2 — 3% of the isospin gauge coupling, but still 
finite. 

As the Three-Site Model must be treated as an effective field theory, a loga- 
rithmic dependence on the UV cutoff scale A remains in the result which fig. 14.31 
demonstrates to be very moderate. Increasing the W mass moves the allowed area 
to higher values of Mbtjk in agreement with the bound (|4.26l) obtained from aT for 
ideal derealization. At some value of the W mass, the bound on Mt, u ik crosses the 
UV cutoff, indicating the breakdown of the theory that contains the bulk fermions 



5 Two comments are in order. First, if A/ bu i k exceeds the UV cutoff, then the above estimate of 
aT from loops with dynamical KK fermions is inconsistent. Instead, an analysis in the effective 
theory obtained by integrating out the heavy fermions from an UV-completed version of the model 
should be performed. This is carried out in 1 4 1 with a toy completion and the result is shown to be 
unchanged. 

Second, the case of ideal derealization is excluded by the one-loop result for aS (see below). How- 
ever, g W i f f has still to be very close to zero, so the above constraints on A4, u ik an d the e'f are still 
realistic. 
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as physical degrees of freedom. While this is of little consequence to tree level cal- 
culations, the effective field theory analysis performed in [35j breaks down at this 
point. 



Triple gauge boson couplings 

After the discovery of the W and Z, LEP-II was the first experiment capable of mea- 
suring the triple gauge boson couplings predicted by nonabelian gauge theory. The 
results put bounds on these couplings which differ between the Three-Site Model 
and the Standard Model. 

In the absence of CP- violation, al@ possible dimension 4 operators that generate 
couplings between two charged and one neutral gauge bosons can be parameter- 
ized in Hagiwara-Zeppenfeld notation |4"T| 



£ 3g = -ie (1 + Ak a ) W+W~A>" - ie (l + Ag£) (W + ^W~ - W~^W+) A v - 

^ (1 + Ak z ) W+W-Z^ - ie^ 
sw s w 



ie— (1 + Ak z ) W+W~ Z» v - ie— (l + Agf ) {W + ^W~ - W~^W+) Z v 



{4,27) 



In order to obtain the values of the couplings and the parameters Ak^/z and Ag l , 
the mass eigenstates W, Z and A must be expressed through the gauge eigenstates 
situated at the lattice sites. The normalization of the wavefunctions together with 
the photon wavefunction (|4.18|) implies Aka = Agf = by electromagnetic gauge 
invariance. The remaining two parameters must be equal and evaluate to 



(1 + Ak z ) = (l + Agf) = gf z (#) 2 + gf* (/f) 2 
sw sw 

Plugging in the wavefunctions and masses tab. 14.11 we obtain 



A ^f = A«z = ^ + O (x 3 ) 
Using (|4.19l) and (14.211) results in a lower bound on the W mass 



.2 



m w 



i = — - 



1 m% 



2 Ag? 

From [42 1, the upper bound on Agf is 0.028, and we obtain 



m W ' > 380 GeV (4.28) 



6 The Hagiwara-Zeppenfeld parameterization excludes operators containing pieces like 9 M zT M 
whose contributions are suppressed with the electron mass at LEP and which are therefore not ob- 
servable in the LEP / LEP-II data. 
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Note that, together with the above bounds on M^ik, this implies a lower bound on 
MbuiJ^ which we can read off from fig. l4.2l as 

M buJk > 1.8 TeV - 2 TeV (4.29) 

(depending on the exact bound on aT assumed, see above). 

4.3 Couplings 

The calculation of couplings in the Three-Site Model follows the pattern already 
shown for gyji f f and gwwz in the last section: the couplings are obtained from 
reexpressing the mass eigenstates (which are delocalized over the whole lattice) in 
terms of the gauge eigenstates localized at each lattice site, e.g. for couplings of 
arity three 

9XYZ = ^gnfn fn fn 
n 

In the language of deconstruction, these sums are the deconstructed version of the 
overlap integrals which give the couplings in the continuous case (c.f. (13.71) ) 

9X YZ = Jdy 9 (y)f x (y)f Y (y) fz (y) 

(allowing for a possible dependence of the 5D couplings on y). 

Gauge sector 

The only free parameter in the gauge sector of the Three-Site Model not present in 
the Standard Model is x and, therefore, all gauge sector couplings must be inde- 
pendent of Mb u ik (and cl in the the case of nonideal derealization). 

As shown at the beginning of the previous section, the Three-Site Model is iden- 
tical to the Standard Model without a Higgs in the limit d4.15|) . This implies that 
couplings between the Standard Model gauge bosons match those in the Standard 
Model up to corrections of order O (x 2 ) or higheid, where x is constrained to be 
smaller than « 0.42 by KM and (|4l9l . 

By the same reasoning, the couplings g and g' at the lattice sites and 2 are 
equal to the isospin and hypercharge gauge couplings up to corrections of order 
O (x 2 ) . The coupling g at the bulk site can be written as 

g g m w > 

9 = ~ ~ 7T 

x 2 raw 

7 An additional bound on Mt, u ik arises from potential dangerous contributions to b — > sj orig- 
inating from the small but finite right-handed Wtb coupling present in the model. However, this 
constraint is weaker than 1 14.291 , see |4 [ for details. 

8 As the sign of the gauge couplings is purely conventional and can be changed by redefining the 
gauge field, all couplings have to be either even or odd in x which is the reason why the corrections 
to couplings between Standard Model particles start at O (x 2 ). 
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Figure 4.4: The four-point coupling between different combinations of W and W' , 
double logarithmic plot to emphasize the leading order dependence on x. 



and exhibits a linear growth with m^y. We can't reasonably expect the delay of 
unitarity to work out of the W is heavier than O (1 TeV) and, therefore 



9 < g\ 



=1.5 TeV 



6.2 < 4vr 



This means that the gauge coupling of the SU(2) at site 1 is moderately strong with 
radiative corrections being roughly of order y|^- m 25% in the worst case. 

In the limit x — > 0, the Standard Model gauge bosons are localized on the branes 
with the fi being of order x, while the W and Z' are localized in the bulk with 
/ /2 of O (x) (c.f. tab. I4.1|) . Naively one would therefore expect that all couplings 
involving different KK modes are suppressed by powers of x. However, this is not 
true because the coupling at the bulk lattice site is proportional to x~ x . Taking this 
into account and recalling the definition of x (|3.34|) , it is not difficult to see that a 
coupling between gauge bosons involving n W' or Z' (for n > 0) must be of order 
O (j; 2_n ) relative to the corresponding Standard Model coupling. This behavior is 
demonstrated in fig. 14.41 

The only exception from this rule are couplings involving only W, W and pho- 
tons. In this case, orthonormality of f w and f w together with the photon wave- 
function (14.181) implies that all couplings apart from 

9ww~i = Qw'w't = e ' 9WW-TY = gw'w'Tf = e2 

vanish. 

As QCD is not touched by the additional structure in the Three-Site Model, the 
couplings between the gluons are unchanged with respect to the Standard Model. 



Fermion sector: ideal derealization 



The fermion sector is more complicated as it depends on x as well as on Mbuik- 
However, the Standard Model must be recovered in the limit (|4.15l) , so all couplings 
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Figure 4.5: Left: Right-handed Wtb coupling relative to the isospin gauge coupling 
g (forbidden in the Standard Model) in the Three-Site Model in as a function of the 
mw' and Mbuik- 

Right: Similar plot of the left-handed Wtb coupling which vanishes for massless 
fermions due to ideal derealization. 

involving only KK light particles must equal the respective Standard Model values 
up to corrections of order O (x 2 ) and O ( e / 2 ) (this also covers the dependence on 
el as it is proportional to x by ideal derealization). 

In the ideally delocalized scenario, x enters the wavefunctions of the fermions 
through their dependence on ei, while M^uik enters through e'j. For a massless fer- 
mion we have e'j = and therefore, the wavefunctions and couplings of massless 
fermions only depend on x. In addition, the right-handed zero modes of the mass- 
less fermions have to completely decouple from their left-handed counterparts and 
therefore are completely localized at site 2, while the KK-heavy modes are com- 
pletely localized in the bulk in this case. 

Consulting tab. 14.21 shows that the massless KK-light fermions are localized at 
the branes with // suppressed like O (cl) while their KK-heavy partners are lo- 
calized in the bulk with f^ 2 similarly suppressed. This implies that the counting 
rule derived for the gauge sector also holds for W coupling to massless fermions or 
their KK-partners. In the case of massive fermions, e'j goes to zero as M^ u \^ goes to 
infinity, recovering the same pattern in this limit and yielding corrections of order 
e'f for finite M^^. 

The left-handed couplings between a W' and two massless fermions are not 
only suppressed by O (x 2 ) but fully vanish as required by ideal derealization due 
to a cancellation between the contributions from different lattice sites. In the case 
of massive Standard Model fermions, the cancellation is not complete, but the re- 
sulting coupling is still very small as demonstrated in fig. 14.51 right by the example 
of Wtb. 

As the right-handed KK-light modes of the massive fermions and the W have 
a nonvanishing overlap in the bulk, a small right-handed coupling between the W 
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Figure 4.6: The left- and right-handed Z'uu couplings relative to their Zuu coun- 
terpart as a function of the W mass mw> ■ Double logarithmic plot. 

and the massive Standard Model fermions suppressed like O ( e / 2 ) is present in 
the Three-Site Model. This coupling is showcased in fig. l4.5l left for Wtb. 

The couplings to the Z are more complicated due to the nontrivial charge of the 
fermions at site and 1 under the U(l) at site 2. Taking this into account, the left- 
and right-handed couplings read 

9z hh ,L =T 3 (g /* ft L f z + g ff L ff L /f ) - Yg'f 2 z ( £ /& /& 

\n=0 

9z f2f2 ,R =T 3 ~g f{ R f( R /f + Qgf fi R f? R ff + Yg> ff f{ R f( R 

The first two terms of both couplings come from the wave function overlap and 
directly fall under the same counting rule as above in the e'* = case. The third 
term comes from the part of the coupling to B 2 which is nonlocal on the lattice and 
could potentially violate our counting rule. To see that this is not the case, note that 
the contributions of these terms are either of order O (x) or O (1), and the latter only 
if the Standard Model Z is involved. However, in both cases, the only combination 
of KK modes that could give aO(l) contribution is Zf'f for which the counting 
rule also gives 0(1) and therefore, the leading order in x as given by the rule is 
unchanged for all possible combinations of modes at the vertex. 

The submatrix of the neutral gauge boson mass matrix M NC (c.f. (|4.2|l ) which 
links the gauge fields at sites and 1 is identical to the charged gauge boson mass 
matrix M cc . Therefore, we can expect the wave function of the Z' at these lattice 
sites to be approximately proportional to that of the W (a conjecture which can be 
confirmed by a look at the actual wave functions in tab. I4.1|) . Thus, the cancellation 
leading to ideal derealization also suppresses the left-handed couplings between 
the Z' and the Standard Model fermions. This is confirmed by fig. 14.61 which shows 
the right- and left-handed couplings of Z'uu relative to their Zuu counterparts. 
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While the right-handed coupling is of order x as expected if there is no special rela- 
tion between the wave functions at different lattice sites, the left-handed coupling is 
suppressed by roughly another order of magnitude by the incomplete cancellation. 

Putting the pieces together, in the limit e'j — > (either by vanishing Standard 
Model mass or large Mb u ik), the couplings of the KK-light fermions to W and Z are 
equal to their Standard Model model values up to corrections of O (x 2 ), while cou- 
plings involving n (n > 1) KK-heavy particles go like O (x 2 ~")B As far as the left- 
handed couplings to the Standard Fermions are concerned, the W' is completely 
fermiophobic and the Z' at least close to it. 

Due to electromagnetic gauge invariance, the photon doesn't couple to fermions 
of different mass. The couplings of a photon to two Standard Model fermions (or 
both the KK partners) are equal to the Standard Model couplings. Similarly, the 
gluons either couple to two KK-light or two KK-heavy quarks with the couplings 
borrowed from the Standard Model. 

Nonideal derealization 

In |35| it is shown that, in contrast to the ideally delocalized scenario, a nonvanish- 
ing coupling between the W and the light fermions is necessary for compatibility 
with the precision observables at the one loop level. This coupling gy/> / / for which 
the bounds are derived is defined in the low-energy effective theory obtained by 
integrating out the heavy fermions and renormalized at the W mass scale. 

There are two operators contributing to gy/ 1 ff in the one loop analysis. The first 
one 

arises from integrating out the bulk fermions from diagrams of the type 




and corresponds to the Yukawa coupling el in the full model which includes the 
heavy fermions. The second operator arises from loop corrections 

2 = ^ol (d„ (Siif E{)) 7 m *ol 

Although this operator also encodes a coupling between the left-handed Standard 
Model fermions and the W, it has a nontrivial momentum structure. However, a 

9 As a slight exception to this rule, several couplings to right-handed fermions vanish completely 
due to the exact localization of the right-handed components of the massless fermions and their KK 
partners. 
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Figure 4.7: Examples of the change induced in the left-handed Zuu type couplings 
by tuning e l always from ideal derealization. 



nonlinear field redefinition in the spirit of on-shell effective field theory B3ll can 
be used to apply the equation of motion to A\ and convert O2 to the same form 
as 0\ at the price of introducing additional higher dimension operators which are 
suppressed by a power of the gauge couplings and two powers of v. After the 
redefinition, O2 also gives a contribution to gy/' / / which is accounted for by in 
the full theory. Therefore, at tree level, the finite gy/ 1 / / can be accommodated in the 
full theory by tuning away from the value required by ideal derealization. 

From the fermion and gauge boson wavef unctions tab. 14.21 and tab. 14.11 we can 
estimate the value of ei required to generate a given g\y> / / 

f „ 9W'ff , x<1 

£l ~ — ; h — 

9 2 

Using the condition for ideal derealization d4.23|) , we can solve for the relative 
change in 

f-, , 9Wff\ 
£L ~ £L,o \l+x ^ J J 

with cl,o being the value required for ideal derealization. According to II35II , gw'ff 
is constrained to be smaller than some 2 — 3% of g, and we find that the required 
change to is of the same order. 

At leading order, the wavefunctions of the right-handed fermions are indepen- 
dent of el (c.f. tab. 14.21) , and therefore, only the left-handed couplings are sensitive 
to the departure from ideal derealization. Naively, we would expect the resulting 
relative change in the left-handed couplings to be of the same order of magnitude 
as the relative change in e^. However, this is not true for the couplings of the Stan- 
dard Model fermions to the W and the Z' . 

For these couplings, the change of disrupts the cancellation among the con- 
tributions from different lattice sites, leading to a potentially much bigger change. 
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This is demonstrated in fig. 14.71 which shows the change induced in several of the 
left-handed Zuu type couplings. All couplings change within a few percent or even 
significantly less, with the exception of the Z'uu coupling which changes by nearly 
300% when e/, is changed by 10%! However, although the relative change is sizable, 
the coupling still remains nearly two orders of magnitude smaller than the isospin 
gauge coupling and therefore close to fermiophobic. As the corresponding right- 
handed couplings to the Z' are considerably larger and not affected by nonideal 
derealization, we shouldn't expect a big change in most physical observables. 

For the W couplings to Standard Model fermions, the situation is different: al- 
though the coupling is still virtually fermiophobic, the change induced by nonideal 
derealization opens up the possibility of producing this resonance from fermions 
which is impossible in the case of the ideally delocalized scenario (c.f. chapter[8]>. 

The gauge sector as well as the couplings to photons and gluons are not affected 
by nonideal derealization. 

Goldstone bosons 

From expanding the Wilson lines Sj / 2 in terms of the component fields and insert- 
ing the expansion into the Three-Site Lagrangian, we would obtain infinitely many 
couplings involving arbitrarily many Goldstone bosons. However, choosing uni- 
tarity gauge completely removes the Goldstone bosons from the spectrum, which 
is what we will be doing for the remainder of this work. 

4.4 Widths and decay channels 

In order to regularize the mass poles appearing in transition amplitudes, we have 
to calculate the widths of the new heavy resonances. In addition, some knowledge 
on the decay channels and branching ratios is required in order to identify suitable 
final states for the discovery of these particles. In this section we will only state the 
results of the calculation, the details of which can be found in appendix [Cl 

Heavy gauge bosons 

At leading order, the W and Z' can either decay into two Standard Model gauge 
bosons or two Standard Model fermions fermions 



All other potential 1—^2 decay channels are forbidden due to the near-degeneracy 





of W and Z' 



1 + O (x 2 ) 



m z > 
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Figure 4.8: Left: Relative widths of W and Z' as a function of mw 1 - 

Right: Branching ratio of the heavy gauge bosons going to fermions as a function 

of mw>- 



and the bound on the heavy fermion mass scale (I4.26D . The resulting relative widths 
are shown in fig. !4.8l left to be very similar for W' and Z' and of the order of 1 — 3%. 

Compared to the fermionic decay channel, the bosonic decay channel receives 
an enhancement factor 



l6my V m 2 z 

which is of order 100 and which comes from the decay into longitudinal modes. 
Therefore, although there are more final states and phasespace available for the 
fermionic channel, the bosonic final state is highly favored. This is demonstrated in 
fig. l4.8l right which shows the fermionic branching ratio of Z' and W as a function 
of mw' (in the case of the latter, cl is tuned away from ideal derealization by 5%). 
Evidently, the W decays to more than 95% into gauge bosons, while the fermionic 
branching ratio of Z' is even lower than 1%. 

Heavy fermions 

For the heavy fermions there are two two possible decay channels at leading order: 




The decay into the heavy isospin partners is forbidden due to the near-degeneracy 
of the heavy fermions. The O (a s ) QCD contributions to the decays of the heavy 
quarks can be expected to be of order of several percent and may therefore be 
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Figure 4.9: Right: Contour lines of constant relative e' width T e / re j = — ^- in the mw' 
- M bvSk plane. 

Left: Relative widths of e', v! , t' and b' (including one loop QCD corrections) as a 
function of A/bulk- 



worthwhile to include. These corrections can be written diagrammatically as 



A\M\ 2 = 




(4.31) 



and are included in the calculation of the heavy fermions widths presented here; 
the analytical result can be found in appendix ICl 

Fig. l4.9l left shows contour lines of constant relative width T e i K \ for the e' in the 
mw> ~ Afbulk plane. Evidently there is a huge range of possible values from about 
10% to more than 100%. However, direct detection as resonances at the LHC will 
only be possible for moderate values of Mt, u ik, say < 4 TeV anyway and in this part 
of parameter space, the relative width does not exceed 30 — 40%. 

The right-hand plot of fig. 14.91 shows the relative width of e', it', b' and t' as a 
function of M\, u \k for mw> = 500 GeV. As the tree-level widths of the heavy fermi- 
ons (with exception of the t' and b') turn out to be virtually identical, the difference 
between the e' and the it' width is solely due to the QCD loop corrections which 
are around 7.5% and show only very little variation over the parameter space. The 
t 1 is considerably broader than the partners of the light quarks, which is due to the 
wavefunctions of the t being fairly nonlocal to accommodate the high top mass, 
enhancing the couplings between t' and t. 
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The big surprise, however, is the b' which turns out to be nearly two times as 
broad as the partners of the light quarks. The reason for this unexpected effect is 
the right-handed W^b't coupling which is considerably enhanced by the nonlocal 
wavefunction of the right-handed top quark. This can be confirmed by looking at 
the sample spectrum in appendix[BJ 

Concerning the branching ratios, fig. I4.10l left shows the fraction of d' decaying 
to heavy gauge bosons as a function of my/ 1 for different values of M^^. While 
this varies a bit over parameter space, it's always around 50%. For d! which de- 
cay versus the light resp. heavy decay channels, fie. l4.10l right shows the fractions 
which goes into a W resp. W , demonstrating these channels to be preferred with 
branching ratios of around 65% over the Z / Z' channels. Both results are repre- 
sentative for the other heavy fermions with little variation over different values of 
Mb u ik, again excluding the t' and b' due to the high top mass. 

For these two, the respective results are shown in fig. 14.111 and fig. 14.121 Al- 
though the functional dependence on my/' and Mb u ik is different, the branching 
ratio of both t' and b' into heavy gauge bosons is around 50% similarly to the part- 
ners of the light fermions. However, the branching ratios into and W are very 
different and exhibit a strong dependence on mw'- Further investigation shows 
that, while still not drastic, the dependence on Mt, u ik is much more pronounced as 
for the other KK fermions. 
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Figure 4.10: Left: The branching ratio of the d! into heavy gauge bosons as a func- 
tion of mw 1 for different values of Mt, u ik- 

Right: The fraction of the d! taking the light / heavy decay channels going to W / 
W. 
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Figure 4.11: The same type of plots as fig. l4.10l for the t' . 
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Figure 4.12: The same type of plots as fig. l4.10l for the b' '. 



Chapter 5 

Tools for Phenomenology 



This admitted, we may propose to execute, by means of machinery, the mechan- 
ical branch of these labours, reserving for pure intellect that which depends on 
the reasoning faculties. 

(L. Menabrea, "Sketch of The Analytical Engine, invented by 
Charles Babbage", translated by Lady A. Lovelace) 

So far, we have dealt with the formal structure of a quantum field theory (the Three- 
Site Model in our case) which defines the theory and allows to calculate observables 
such as scattering probabilities. To test the theory, we must connect these observ- 
ables to the quantities that are measured in collider experiments or, more specific to 
this work, at the LHC. This is usually done with the help of computers via Monte 
Carlo eventgenerators. 

In this chapter, we first briefly review the general structure of an event simu- 
lation. After that, the WHIZARD / O'Mega package which has been used for the 
numerical simulations presented in this thesis is introduced together with a small 
review of the challenges which have to be tackled by Monte-Carlo eventgenerators. 
In the third section, the implementation of the Three-Site Model into WHIZARD / 
O'Mega done as part of this thesis is discussed. Finally, in the last section, a brief 
overview over an interface between FeynRules [44J and WHIZARD / O'Mega is 
given. 

5.1 The structure of an event simulation 

The simulation of collision events at particle colliders can be roughly split into three 
building blocks. The first is the calculation and integration of the differential cross 
section which consists of evaluating the squared matrix elements and performing 
the phase space integration. The second block is the actual generation of an ensem- 
ble of events that is distributed according to the differential cross section. The third 
block consists of simulating the characteristics of the detector used for the actual 
observation of the particles produced in the events. 
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5.1 THE STRUCTURE OF AN EVENT SIMULATION 



Obviously, the calculation and integration of the differential cross section is the 
necessary prerequisite for the generation of events. However, the detector simu- 
lation part completely factorizes from the other two steps. As it also is time con- 
suming and requires an intricate knowledge of the inner workings of the actual ma- 
chinery employed in the experiment, it is usually a good idea to perform studies on 
the feasibility of detecting new physics without a detailed detector simulation first. 
The result then is a machine-independent assessment of the chances for detection 
which, if it is positive, can be further improved by adding a detector simulation. 

The calculation of the cross section also is an intricate process. The complica- 
tions start with the initial state which, at a hadron collider like the LHC, consists 
of complicated bound states of quarks and gluons which we cannot describe from 
first principles. The final state usually consist of leptons, photons and mesons, and 
the latter again are bound states of QCD whose formation cannot be described with 
perturbation theory. Luckily, nature has arranged for asymptotic freedom and the 
factorization theorems of QCD |45| which allow us to split this complicated pro- 
cess of hadrons — > mesons + X into different regimes which decouple in good 
approximation. 

More specifically, we may describe the initial hadrons at high energies as en- 
sembles of non-interacting partons with the probability of finding a given parton 
carrying a given momentum fraction being described by parton distribution func- 
tions (PDFs) which are universal and independent of the process under considera- 
tion. In addition, the QCD interactions that lead to the fragmentation of the hard 
partons and the final hadronization to the observable bound states can also be ap- 
proximately factorized. 

In the end, as a good approximation, we can reduce the complicated scatter- 
ing process to a "hard" core process consisting of two ingoing partons scattering 
into a handful of hard outgoing partons, photons and leptons. This then is convo- 
luted with the PDFs to accommodate the hadronic initial state, and parton showers 
are added to describe the evolution of the small number of hard partons to a large 
number of soft ones. In a final step, empirical models can be applied to describe the 
formation of bound states from these soft partons in a process called "hadroniza- 



The only piece of this decomposition which we can calculate in perturbation 
theory is the hard matrix element describing the partonic core process (which for- 
tunately usually also is the only part sensitive to new physics). Algorithms for 
performing the parton shower can also be obtained from first principles. However, 
there is to this day no way of calculating the PDFs from first principles, and only 
empirical models exist for the hadronization process. 

In this work, we concentrate on the hard process and don't include parton 
showers, hadronization or a detector simulation. This is called a "parton-level" 
simulation which essentially relies upon the assumption that every parton pro- 
duced in the hard scattering event is distinctly observable as a jet in experiment. 
While this is certainly not correct, it is a good assumptions for embarking on a 
study of the phenomenology of the Three-Site Model: it is only in the case of signif- 



tion". 
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icant deviations from the Standard Model at this level that we can hope to observe 
something at the LHC. 

5.2 WHIZARD / O'Mega 

In this thesis, the software package WHIZARD / O'Mega [46J is used for the simu- 
lation of parton level events. This package actually consists of at least two separate 
programs, the optimizing matrix element compiler O'Mega and the core Monte 
Carlo and infrastructure code of WHIZARD. 

O'Mega: fast generation of tree level matrix elements 

Easy as it may seem from calculating the matrix elements for 2 — > 2 by straightfor- 
ward application of the Feynman rules, the evaluation of tree level matrix elements 
actually is very demanding on computer hardware and requires clever algorithms 
if it is to run in finite time. 

The reason for this is the enormous growth of the number of Feynman diagrams 
if we increase the number of external lines. Consider, for example, <f> 3 theory^. If we 
know all the tree level diagrams contributing to the n-point correlator of n particles 
with momenta pi, . . . ,p n , we can easily obtain those for the n + 1-point function 
as all possible ways of attaching the p n +i leg to one of the lines in the n-point 
diagrams, e.g. 



Pi 




Therefore, we can calculate the number N n+ \ of diagrams contributing to the n + 1- 
point function from the number N n of n-point diagrams as 

N n+1 = L n ■ N n 

with the number of lines in any n-point diagram L n . As attaching a leg increases 
L n by 2, we deduce 

L n = 2n — 3 

and therefore find 

71-1 

N n = (2n - 5)iV n _i = (2n - 5)(2n - 7) ■ . . . • 1 = JJ (2k - 3) 

k=2 

1 While this is a pathological example of a field theory due to the indefinite potential, it actually 
is the easiest example when it comes to the combinatorics of Feynman diagrams. The only difference 
to other, physical models with arity three couplings (like QED) is the "dressing" of the ^ 3 -topologies 
with particle flavors which doesn't affect the asymptotic combinatorics. 
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5.2 WHIZARD / O'MEGA 



By induction we have 

N n = (2n - 5)!! = > 2 n " 3 (n - 3)! 

v ' 2 n - 2 (n - 2)1 ~ v ; 

Thus, the number of diagrams grows like factorial , and while we only have three 
diagrams for n = 4, we already have 105 diagrams for n = 6 and 10395 diagrams 
for n = 8. 

Keeping in mind that the matrix element has to be evaluated at every phase 
space point under consideration it is obvious that this factorial growth of the num- 
ber of Feynman diagrams kills off any attempt to perform a brute-force, diagram- 
by-diagram calculation. Things get even worse if a straightforward analytic calcu- 
lation is attempted and therefore, traditional tools like CalcHEP [47J which perform 
very well for 2 — > 2 processes have severe problems handling more than four par- 
ticles in the final state. 

Luckily, this is not the end of the story. If we consider again the addition of 
an additional leg to a diagram, we note that recalculating the whole diagram is a 
waste of resources. Reconsider the insertion of a line 

Evidently we don't have to recalculate the blobs A and B, but just have to cut the 
line connecting the blobs and reattach them to a vertex together with the new leg. 
Also, if A and B are not only pieces of Feynman diagrams but instead complete 
correlation functions with all legs amputated except the leg extending to the mid- 
dle vertex in the above diagrams (so-called "one particle off-shell wavef unctions" 
lPOWs), we can calculate a whole class of diagrams in this way. This observation 
suggests an iterative approach to the calculation of Feynman amplitudes which is 
potentially much more efficient than a straightforward calculation of all Feynman 
diagrams. 

The potential savings in choosing a clever algorithm can be acknowledged by 
counting the number of possible lPOWs. In our example of 3 theory, these are 
uniquely enumerated by the distinct nonzero combinations of momenta which can 
be formed from external lines, the number of which 

^G) + Q + - + (n-l)^ 1 + 1) "- 2 -"- 2 

only grows exponentially. Algorithms like the one used in the optimizing matrix 
element generator O'Mega exploit this redundancy to reduce the factorial growth 
in complexity of the calculation to an exponential one 11481 . 

O'Mega represents the amplitude as a directed acyclical graph (DAG) of lPOWs 
from which it weeds out all redundancies that would arise from the recalculation 



CHAPTER 5. TOOLS FOR PHENOMENOLOGY 



69 



of lPOWs. Finally, this DAG is transferred to a backend which translates it into a 
representation suitable for performing numerical calculations. The only complete 
backend included in the O'Mega distribution is a FORTRAN 95 backend which 
emits a process library coded in FORTRAN 95 that can be called to evaluate the 
amplitude. The resulting code is clean, human readable, fast and can be automat- 
ically instrumented for numerical checks of gauge invariance. The matrix element 
is calculated in a completely numerical way by numerically multiplying spinors, 
four-vectors, Lorentz tensors etc. 

WHIZARD: phase space, Monte Carlo integration and more 

Once we have a way of calculating the matrix elements, we can go over to integrat- 
ing them over phasespace, a step which turns out to be challenging for a number 
of reasons. 

First, the dimensionality is high — subtracting four degrees of freedom due to 
momentum conservation and one due to rotational invariance around the beam 
axis, we have 3n — 5 phasespace dimensions for n final state particles plus two ad- 
ditional dimensions for the parton distributions. Therefore, in a 2 — > 4 process, we 
already have a 9-dimensional integral to evaluate. The complexity of deterministic 
numerical methods (e.g. Gauss-Kronrod) grows exponentially with the number of 
dimensions and therefore, such algorithms are not suitable for this type of integral. 

This problem can be overcome by Monte Carlo methods. A simple method for 
Monte Carlo integration is drawing N random points p^ in the integration domain 
Cl and then approximating the integral of / by 

r v N 

/ dxf(x)KV(f) = -Y,f(P)k 
Jn k=i 

with the volume V of Cl. 

If / is a probability density bounded from above by a number Fq, then we can 
use the set of pairs 

W = {(p 1 ,/(p 1 )),...,(p A r,/(Piv))} 

to obtain an ensemble of points distributed according to / by drawing N numbers 
Zk between and F and then discarding all points pk for which f(pk) exceeds z^. 
The resulting set of n surviving points 

U = {pifc 1 ,...,Pfe„} = [pk | f(pk) < z k } 

can be significantly smaller than W but will obey the desired distribution. The set 
W is usually called a set of "weighted events" with statistical weights given by the 
f(pk), while U is called a set of "unweighted events" with the processes leading 
from W to U termed as "unweighting". If we want to simulate events at a particle 
collider, then it is the unweighted events we are interested in as they give a realistic 
estimate of the fluctuations induced by the statistics. 
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It can be shown that this Monte Carlo integration comes with an error that scales 
with -7= for large N independently of the number of dimensions of the integral. 
The coefficient a is determined from the variance functional A 



A[/] = f dx f(xf -( [ dx f(x)\ 



Evidently, this is a small number only if the function / is well-behaved over SI (aka 
doesn't fluctuate much). 

However, the differential cross section we want to integrate is far from well- 
behaved: by the virtue of the propagators which have mass singularities, it fluctu- 
ates wildly, being close to zero in large regions of phasespace and nearly singular 
in others. The obvious problems arising from this large variation can be alleviated 
by choosing the ensemble of points p n not uniformly distributed in Q but instead 
with respect to a nontrivial probability density g(x). Provided we perform the re- 
placement 

the Monte-Carlo method presented remains correct, but with the error coefficient a 
now determined by A ^ . 

Choosing a suitable g allows to reduce the Monte Carlo error drastically at the 
expense of more calculational time spent on the generation of the pj. Therefore, g 
should be chosen such that the generation of the points is cheap enough to retain 
the speed gained by reducing A. As manual tuning of g usually is too cumber- 
some, a suitable density is usually determined in an adaptive process designed to 
minimize A over a restricted function spaced. 

Monte-Carlo eventgenerators like WHIZARD usually start with an adaption 
phase during which g is optimized, then proceed with an integration phase where 
the integral is obtained and finally an event generation phase where (unweighted) 
events are generated. Ideally, new points are generated in all three phases to avoid 
pollution of the error estimates by statistical correlations. 

In order to choose a suitable class of functions over which the density g is op- 
timized, WHIZARD has some limited knowledge of diagrammatics and needs in- 
formation on the vertices of the model. Prior to the start of the adaption phase 
WHIZARD then classifies all singular regions in phasespace that might contribute 
to the integral ("integration channels"). This information is passed to the Monte- 
Carlo core VAMP (491 (a multichannel modification of the VEGAS algorithm 115010 
which chooses a suitable phasespace map <pk for each channel. In these coordinates, 



2 It is an interesting fact that, for a restricted function space, the qa which minimizes the variance 
functional A usually does not maximize the acceptance in the unweighting step (the percentage of 
points in W that is kept). It is also possible to optimize the acceptance instead of the variation in 
the adaption phase, increasing the performance of the unweighting step at the price of reducing the 
accuracy of the estimate for the integral. 
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a multidimensional grid function together with a weight a n is assigned to each 
channel, and the full parameterization of g reads 



with the sum running over all channels. 

In the adaption process, g is optimized by adjusting the grids g n as well as the 
weights a n . After the adaption, the best g obtained in the adaption phase is used 
to calculate the Monte Carlo integral of the cross section. Finally, WHIZARD can 
use the grid and the integral to generate an ensemble of unweighted events which 
correspond to a given integrated luminosity. 

The emphasis of the WHIZARD package lies on the fast generation of parton 
level events for BSM physics. Although its primary source of matrix elements is 
O'Mega, it can also interface other matrix element code^l like MadGraph |5T| or 
CalcHep. It is capable of convoluting the partonic cross section with PDFs (via 
PDFLIB 1 52] or LHAPDF (51) and can interface PYTHIA (5j| for parton shower- 
ing. The high numerical speed achieved by the combination WHIZARD / O'Mega 
allows for the simulation of processes with 6 — 8 final state particles which is very 
hard to achieve with other tools. Other features include an integrated facility for 
event analysis which allows to create histograms on-the-fly from the generated 
event data. 

The package is written in FORTRAN 95 with some perl glue being used to as- 
semble the source and the matrix element code prior to compilation and evalua- 
tior@. For the exact treatment of color, the color flow decomposition [55| of Feyn- 
man amplitudes is used. 

5.3 Implementation 

In order to study the phenomenology of the Three-Site model with WHIZARD / 
O'Mega we have implemented the model into this package. In its current state, 
the code implements the Three-Site Model in unitarity gauge without flavor mix- 
ing and supports both ideal and non-ideal derealization. It has been validated by 
numerical checks of the Ward identities in the limit v — > as well as by direct com- 
parison to the CalcHEP implementation used by |34] and (lately) to the FeynRules 
version done by N. Christensen et al. (see chapter I5T4)). 

The model implementation will be part of a future version of WHIZARD. In 
the meantime, the coupling library as well as modified versions of the O'Mega and 
WHIZARD packages which include the model can be downloaded from 

http: //theorie. physik .uni-wuerzburg.de/~ enspeckn/ 

3 This functionality is present in the 1.9x branch of WHIZARD but will be dropped from version 
2.0 onwards. 

4 This is only true for the 1.9x branch. WHIZARD 2.0 and higher compile the matrix elements at 
run time and link them dynamically without relying the perl component which is dropped. 
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5.3 IMPLEMENTATION 



O'Mega 

In order to perform its duties as a matrix element generator, O'Mega needs some 
information about the model under consideration. This comes under the guise of a 
O'Caml module which contains the necessary definitions and functions needed by 
O'Mega to operate. We have written such a module implementing the Three-Site 
Model. 

The implementation offers several options that can be (de)activated at the users 
leisure to change the set of particles vertices included. More specifically it is pos- 
sible to allow for flavor off-diagonal couplings (although these are set to zero at 
the moment, so this would currently just be a waste of resources), to leave out the 
heavy fermions and to discard the couplings of the W to the leptons and to the 
first two quark generations. In addition, it is possible to replace the unitarity gauge 
propagators with Feynman gauge ones (this is necessary for using the massless 
limit). 

Upon compilation of O'Mega, binaries are generated from the model modules 
which can then be called to generate code for the desired matrix elements. For the 
Three-Site Model, different binaries are generated for the different vertex sets de- 
scribed above both in colored and colorless versions (the colorless ones completely 
exclude QCD). A command line switch can be used to replace the unitarity gauge 
propagators. 

More details on the implementation together with some snippets of the actual 
code can be found in app. ID.2I 

WHIZARD 

In addition to the model implementation in O'Mega, WHIZARD requires a model 
file which defines the particles and vertices present in the model together with all 
parameters that can be changed at run time (and optionally also a list of secondary 
parameters calculated from these). 

Writing the model file is a straightforward task. The only part which requires 
some trickery is the vertex list which WHIZARD requires to be loosely sorted with 
respect to the mass of the particles meeting at the vertex. Even if flavor violating 
couplings are excluded, the Three-Site Model contains 220 electroweak three-point 
couplings. To avoid the cumbersome and error-prone task of collecting and sorting 
all these couplings per hand, the list in the model file is generated directly from the 
O'Mega module by an O'Caml script. 

The resulting vertex list may be too long for processes with 6+ particle final 
states, causing the phasespace mapping step to take forever or fail altogether. In 
this case, removing all vertices containing KK fermions provides a workaround 
without affecting the quality of the result overmuch. 

The parameters^ defined in the model file allow to set all Standard Model masses, 
the t, W and Z widths, the electromagnetic and strong coupling constants, the W 



5 The actual names of the parameters are documented in the model file. 
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mass and M\, u \^. This is enough to fix the model in the case of ideal derealization 
(c.f. chapter l4!2|) . Nonideal derealization can be activated by toggling a flag^t mak- 
ing it necessary to supply a value ei to complete the set of model parameters in this 
case. The widths of the KK partners and the couplings are automatically calculated 
from the input parameters (see below). The inclusion of the QCD corrections into 
the heavy quark widths can be toggled by a flag, and another flag can be used to 
dump the spectrum and couplings (similar appendix|B]) to a file. 

Unfortunately, the architecture of WHIZARD in its current form (version 1.9x) 
requires some patches to WHIZARD itself to make the model work. However, in 
a future version of WHIZARD, it will be possible to use the model without any 
changes to the actual WHIZARD code, allowing for the integration of the Three- 
Site Model into the package. More details on the required changes can be found in 
app.|D3j 

Coupling library 

To make the implementation work, a helper is required which calculates all masses, 
couplings and widths from the input parameters (c.f. chapterH]). This task is hand- 
led by a dedicated FORTRAN 90 library. 

This library takes the input parameters and uses them to first calculate the 
masses and wavefunctions of the particles. The calculation is performed using ex- 
act analytic formulae in contrast to the expansions in x and €l given in chapter l4~2l 
The wavefunctions are then used to calculate the couplings as shown in chapter [43] 
(with the exception of the couplings to the photon which are determined by elec- 
tromagnetic gauge invariance). After the couplings are calculated, the library then 
calculates the widths of the KK particles by looping over all possible two particle 
final states, using the analytical formulae given in app. IC.ll 

If activated, the O (a s ) corrections to the heavy are also included using the an- 
alytic expressions of app. IC.2l In this case, the LoopTools library [56] is utilized for 
the numerical evaluation of the one loop integrals appearing in the virtual correc- 
tions and for the dilogarithm appearing in the three-particle phase space integrals. 

In order to validate the O'Mega implementation of the model by numerical 
checks of the Ward identities, the library supports a setup in which the full SU(2)o x 
SU(2)i x U(l)2 gauge symmetry is left unbroken. In this mode, the model is ini- 
tialized from x, e, t and an arbitrary angle <j). The scale v and the Yukawa couplings 
A, A, A' are set to zero, resulting in vanishing mass matrices. The gauge couplings 
are initialized as 



As the mass matrices vanish, we can choose the wavefunctions freely (provided 
the wavefunction for KK partners are orthonormal). The W / W and the fermion 

6 As WHIZARD only supports real numbers as parameters, a negative number is read as "false" 
and a positive number as "true". 



g = ey 1 + x 2 + t 2 




g = gt 
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wavefunctions are parameterized by the angle 

f w = f f L = 

f f R = 



COS ( 

sin< 



sm< 

COS ( 



— sine; 

COS (ft 

— COS ( 

sin <b 



The wavefunction of the photon is left unchanged to avoid a change in the cou- 
plings to the photon (which are hardcoded), and the Z and Z' wavefunctions are 
chosen to form an orthonormal system together with p 
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with the normalization factors being suppressed for the sake of readability. This 
setup tries to avoid accidental cancellations in the couplings which would impair 
the effectiveness of the desired consistency check. 

The library is not tied to WHIZARD but can also be used in a standalone fash- 
ion. Several tools are included which generate spectra (like e.g. app. [B), density 
plots (e.g. fig. I4.5|) etc. together with a wrapper that allows to access the masses, 
widths and couplings from Mathematica (in fact, many of the plots shown in chap- 
ter H] were generated from Mathematica this way). 

More information about implementation and usage of the library can be found 
in app. lD.ll together with some snippets of code; more documentation can be found 
in the README file and in the woven source contained in the distribution which 
can be downloaded from the URL quoted at the beginning of this chapter. 



5.4 FeynRules — > WHIZARD interface driver 

The implementation of a new model into an eventgenerator usually is a time- 
consuming and error-prone process, and careful testing is necessary to ensure that 
everything works as expected. The situation isn't improved by the fact that every 
eventgenerator has its very own format for specifying new models, follows its own 
conventions and has its own deficiencies and weak points which need to be worked 
around. Therefore, a generic format for specifying new models models which is as 
close as possible to the formal representation by a set of fields and a Lagrangian 
would be desirable. 

The FeynRules package |44| is a step in this direction. FeynRules is a Mathe- 
matica package which can be used for the automatic generation of model files for 
different Monte Carlo eventgenerators. To this end, the particles in the model are 
specified in a format which is derived from FeynArts [57]. FeynRules then takes the 
Lagrangian (which can be entered in a way which is very close to the usual formal 
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notation) and calculates the Feynman rules from it. The generated rules can then 
be passed to an interface driver which generates a model file suitable for use with 
the event generator of choice. 

This way, the cumbersome task of implementing the model has to be done just 
once and can be performed at a fairly abstract level with all the particularities of 
the different generators being take care of by the interface drivers. Having the 
model available in different eventgenerators is a very important step as all of them 
have their own strengths and weaknesses and cross-checking a result over different 
generators greatly improves the reliability of a prediction. 

So far, the list of eventgenerators officially supported by FeynRules consists of 
CompHEP / CalcHEP, Sherpa |58| and MadGraph. During the preparation of this 
thesis, a new interface driver has been developed in corporation with N. Chris- 
tensen which generates output for WHIZARD / O'Mega which will be included in 
a future distribution of FeynRules. 

The driver is capable of producing the model module for O'Mega together with 
the model file and FORTRAN glue required by WHIZARD. It supports nearly all 
the features offered by FeynRules, including the free choice of any gauge. The 
generated code has been successfully validated for the Standard Model and for the 
Three-Site Model, validation for the MSSM is in progress. 
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Chapter 6 

W' Strahlung 



After giving an overview over the model and the tools used for studying it in the 
last chapters, we now move on to presenting the results of the actual Monte Carlo 
simulations performed in this thesis. In this chapter, we start off with a discussion 
of the production of W which are radiated from a light gauge boson in a process 
similar to Higgsstrahlung 1159 II . The process is presented in the first section, the 
details of the simulation are discussed in the second section and the result of the 
simulation is presented in the last part. 

6.1 The Process 

As the Three-Site Model is designed to mimic the Standard Model up to small cor- 
rections as far as the Standard Model particles are concerned, the most prominent 
signatures of the model at the LHC should be the resonances of the new heavy 
particles. While the LHC runs at an energy of 14 TeV in the proton center-of-mass 
system, the energy typically available for the partonic hard process is much smaller 
due to the parton distributions and only in the range of several TeV. Therefore, the 
heavy fermions will be hard to detect as resonances due their large mass > 1.8 TeV 
(c.f. chapter©, and the W'/Z' are thus the preferred candidates for direct detection. 

The couplings between the heavy gauge bosons and the Standard Model fer- 
mions are suppressed due to their fermiophobic nature, and diagrams containing 
KK fermions are suppressed with the bulk mass at the W'/Z' mass scale. 

Therefore, the processes that first spring into mind as promising candidates for 
producing the heavy gauge bosons on-shell involve the coupling of the KK particle 
to a W or Z line. Of this type there are basically two kinds of processes which are 
shown in fig. 16.11 namely the radiation from SM gauge bosons or the fusion of two 
SM gauge bosons. 

As shown in chapter 14.41 the heavy gauge bosons decay to nearly 100% into 
Standard Model gauge bosons which then subsequently decay into fermions. Ex- 
cluding all final states with more than two jets due to the large QCD backgrounds 
and all final states with more than one neutrino due to the missing momentum in- 
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Figure 6.1: Diagrams contribution to the signal for the production of KK gauge 
bosons in strahlung type (left) and fusion type (right) processes. 



formation, the following final states suitable for W'/Z' production remain: 



strahlung 


fusion 


W 


Z' 


W 


Z' 


jj + 4/ 

jj + 31 + V 

51 + u 




jj + 3/ + v 





Thus, W strahlung seems to be the most promising candidate as it can be stud- 
ied without involving any neutrinos, the missing momentum information of which 
would have to be made up for with tricks like using the transverse mass 11601 instead 
of the invariant one or reconstructing the neutrino momentum from kinematic re- 
lations (see chapter 0. 

The only source of notable contributions from the additional structure in the 
Three-Site Model to this process is the gauge sector, and therefore, the W mass mw> 
is the only free parameter which has an influence on the result of the simulation. In 
particular, we are free to perform the simulations in the ideally delocalized scenario 
and can infer the dependence of the result on parameter space solely from varying 

In the context of this thesis, a simulation of the W strahlung process via the 
jj + 4/ final state was performed using WHIZARD / O'Mega and the implemen- 
tation of the Three-Site Model presented in chapter |5l Unfortunately, after the sim- 
ulation was started, an already finished study of W production via strahlung and 
fusion was published in [34 J. However, the simulation performed in the context 
of this thesis still is of some value as it simulates the full six particle partonic final 
state using unweighted events, while the simulation performed in above reference 
uses CalcHEP and only simulates the four particle intermediate state jjZZ using 
weighted events. 



1 Tuning e l away from ideal derealization in principle allows for an additional contribution from 
a W' in the s channel, the coupling of which to the radiated W' is larger than that of the W. However, 
the coupling of the W' to the Standard Model fermions is still so small in this that the resulting change 
to the cross section in the peak region is only about 20% at best. 
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6.2 Simulation Details 

In order to study the visibility of the W resonance at the LHC, we have performed 
a full parton level simulation of the process pp — > jj + 41 with 

pe{g,u,d} , j £ {g,u,d,c,s,b} , le{e~,n~} (6.1) 

(and any of the respective antiparticles) for a total energy of 14 TeV in the proton 
CMS. Unweighted events have been generated for an integrated luminosity of fC = 
lOOfb -1 . The parton distributions were taken from the CTEQ6M series [61 J and 
evaluated at the parton CMS energy in a running scale scheme. 
All diagrams contributing to the signal are of the form 



Q 3 




Q r 



(the three-point couplings of the photon don't mix KK modes) and therefore, in- 
variant mass cuts can be used to reduce the non-resonant contributions to the back- 
ground. In particular, we have applied a cut 

60 GeV < rrijj < 100 GeV 

on the invariant mass of the jet pair and a cut 

71 GeV < m u < 111 GeV 

on the invariant mass of the lepton pairs, trying both possible combinations in the 
case of four leptons of the same generation and discarding events where no unique 
identification of the Z bosons could be achieved this way. 
To further suppress the background, a px cut of 

p T > 20 GeV 

was imposed on the final state particles. Another cut 

-0.99 < cos6> < 0.99 

was applied to the polar and intermediate angles of all final state particles, and in 
order to avoid infrared singularities, we imposed the condition 



E > 10 GeV 
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Figure 6.2: Left: W' resonance in the invariant mass of the Z/W pair for the W 

strahlung process; different values of my/i and fC = 100 fb -1 . 

Right: The W resonance for m w > = 600 GeV and fC = 100 fb" 1 , 400 fb" 1 . 



on the energy of the incoming partons which corresponds to a small x cut 

x > 1.4- 10~ 3 

(x being the fraction of the proton momentum which is carried by the parton). 

As there is no way of determining which Z originated from the W', we have 
counted both possible combinations into the histograms similar to the analysis in 
l34l . Simulations have been performed in the ideally delocalized scenario at three 
points in parameter space 

{m w ,,M hulk ) € {(380 GeV, 3.5 TeV), (500 GeV, 3.5 TeV), (600 GeV, 4.3 TeV)} 

(6.2) 

where Mt, u ik has been varied only for consistency in order to satisfy the precision 
constraints as discussed in chapter l42l 



6.3 Results and Conclusions 

Fig. l6.2l left shows the W resonance peaks in the invariant mass distribution of the 
W/Z pair obtained from the simulation as described above. For mw = 380 GeV 
and mw = 500 GeV, the peaks are clearly visible above the background, while 
only a very small bump is present in the data for mw = 600 GeV. This bump 
is reproduced in fig. 16.21 right together with the corresponding distribution for a 
higher integrated luminosity of 400 fb -1 . With this increased luminosity, the mw = 
600 GeV resonance can also be clearly discerned. 

In order to get a quantitative estimate of the significance of the resonances, let's 
define the signal N s as the number iV of events in the ±20 GeV region around 
the peak in the Three-Site Model minus the the same quantity iVj, in the Standard 
Model. 

N S = N-N b 
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Figure 6.3: The statistical significance of the W resonance peaks in W' strahlung as 
a function of the integrated luminosity together with 3a and 5a discovery thresh- 
olds. 

We can then define the statistical significance of the signal as the deviation from the 
Standard Model relative to the standard deviation 

s= N * _ N-N b 
&N b \fN~b 

which scales like \J~f£~ 

For the determination of the background Nf,, we have performed a Standard 
Model simulation of the process for an integrated luminosity of fC = 4000 fb _1 
and downscaled the result by a factor of 40 in order to reduce the error on the 
significance induced by uncertainties in the background estimate. Fig. I6.3I shows 
the resulting significance as a function of the integrated luminosity together with 
the 3a and 5a discovery thresholds. The result for mw> = 380 GeV, 500 GeV looks 
encouraging with a 5a discovery being possible within the first 7 fb _1 resp. 25 fb -1 . 
For mw' = 600 GeV however, things are different, and nearly 300 fb -1 would be 
required for a 5a discovery in this channel, which is (hopefully) at least still within 
LHC running time. 

These results should be compared to the results obtained from simulating the 
ZZjj intermediate state in 11341. To this end, the W resonance for m\y' = 500 GeV 
published by these authors is reprinted in fig. 16.41 left. The number of events con- 
tained in the peak is about 30 which is slightly smaller than the same number in our 
simulation (50 events, c.f. fig. l6.2l left). However, as the CalcHEP implementation of 
the model used in 11341 and our WHIZARD / O'Mega implementation completely 
agree on the 2—^2 process pp — > W'Z and have also been validated against each 
other in numerous other processes, this discrepance is not overly disturbing and 
can be attributed to differences in the cuts, parameters and to the simulation of an 
intermediate state in |34|. 
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200 400 600 800 1000 1200 400 500 600 700 800 900 1000 



M(Z i/) GeV M W] (GeV) 

Figure 6.4: Left: The W resonance peak in the W strahlung process for my/i = 
500 GeV as simulated in [34 1, plot taken from the same reference. 
Right: The integrated luminosity required for a 3a resp 5a discovery as a function 
of the W mass as obtained in 1341 , plot taken from the same reference. WqZqZq 
refers to the W strahlung process, while WqZqjj refers to the fusion process (see 
chapter [67TT> . 



Reprinted in fig. !6.4l right is the integrated luminosity required for a 3a or 5a dis- 
covery as a function of the W mass as estimated in 11341 , the WqZqZq line referring 
to the W strahlung process under consideration. The 5a results are 10 fb , 20 fb -1 
for my/' = 380 fb -1 , 500 fb -1 and roughly consistent with our results quoted above. 
However, according to the authors of I34|, a 5a discovery for my/ 1 = 600 GeV 
should be possible within 100 fb -1 which is much more optimistic than our esti- 
mate of 3001b- 1 . 

Unfortunately, |34| gives no details about how fig. 16.41 right was obtained and 
whether the higher mass region of the curve reflects the results of actual simula- 
tions or is an extrapolation. However, the strong decay of the signal observed in 
our simulations when going from my/' = 380 GeV to my/' = 600 GeV which leads 
to the high luminosity required for discovery is completely consistent with the re- 
sult from performing a simulation of pp —7 W'Z. 

Concluding the discussion of our W' strahlung simulation, the discovery of the 
W in this process at the LHC would be possible for the whole segment of param- 
eter space shown to be compatible with the electroweak precision observables in 
l35l . However, depending on the W mass, a huge amount of data taking might 
be necessary before any conclusive result is within reach. Another disadvantage of 
this process is that it only probes the electroweak structure of the Three-Site Model 
which coincides with other models, e.g. BESS II33II . In order to probe the fermion 
sector of the Three-Site Model which distinguishes it from other scenarios, other 
processes like the s channel production discussed in the next chapter must be ex- 
plored. 



Chapter 7 

KK Gauge Bosons in the s Channel 



In the previous chapter, a study of W strahlung was presented. The major draw- 
back of this process is that it only probes the gauge sector and doesn't reveal any 
information about the fermion sector. In sharp contrast, the production of the W 
and Z' in the s channel is highly sensitive to the structure of the fermion sector due 
to the fermiophobic nature of the couplings between the KK gauge bosons and the 
light fermions. In this chapter, we discuss the prospects of exploiting this kind of 
process at the LHC. 

After presenting an overview over the different channels of this type that might 
be accessible at the LHC in l7.1[ we move on to the details of our simulations in !7.21 
In 17.31 and 17.41 respectively, the results for Z' and W production are discussed. It 
turns out that there are contributions from both W' and Z' to jjlu type final states, 
and a possible strategy do disentangle these is presented in 17.51 The chapter is 
concluded in !7.6t and 17.71 contains a small collection of additional plots which have 
been omitted from the other sections for the sake of clarity. 

Most of the work presented in this chapter is also covered in 162]. The results 
presented here at the parton level were supplemented by a detector simulation in 
the master's thesis of F. Bach II63II . 

7.1 The Processes 

The main motivation behind the Three-Site Model is the delay of unitarity violation 
without introducing a Higgs while still satisfying the precision constraints imposed 
by the LEP / LEP-II data. As discussed in chapters |2] and |3l this implies strong 
constraints on the couplings between the KK gauge bosons (which are responsible 
for delaying unitarity violation) and the Standard Model fermions in this scenario. 
Therefore, if the LHC were to discover new resonances compatible with a W and a 
Z' , further tests of the fermion sector would be necessary in order to probe for such 
a type of new physics. 

While sufficient for the bare discovery of the W resonance, the W strahlung 
process discussed in chapter [6] does not probe the structure of the fermion sector 
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Figure 7.1: g W 'ff, 9z'uu,l and g z >dd,L normalized to the site gauge coupling g as 
a function of the derealization parameter sl- The gray rectangle marks the range 
for gw'ff allowed by the precision observables as derived by the authors of |35|. 



and needs to be supplemented by other channels more sensitive to this part of the 
Three-Site Model. At a hadron collider, the simplest process involving the coupling 
of the W'/Z' to the light fermions is the production of these particles in the s chan- 
nel. As the KK gauge bosons decay to nearly 100 percent into light gauge bosons 
(c.f. chapter |4~4)| , the signal in this kind of process is provided by diagrams of the 
type 




which involve one order of the heavily constrained couplings between light fer- 
mions and KK gauge bosons. Depending on the magnitude of those couplings, 
the total invariant mass of such final states should exhibit a W'/Z' resonance, the 
magnitude of which depends on the fermion (de)localization. 

As discussed in chapter 14.31 the relevant couplings are restricted to be non- 
vanishing at the order of several percent of the isospin gauge coupling. Fig. 17.11 
shows the coupling of the left-handed^ light quarks to the W and to the Z' for 
mw' = 500 GeV as a function of the fermion derealization parameter (c.f. chap- 
ter [42]), the corresponding plots for my/' = 380 GeV, 600 GeV can be found in fig. 
17.111 Remarkably, the left-handed coupling of the Z' to the up quark gz>uu,L and 
that to the down quark gz'dd,L change in exactly opposite ways under a variation of 



1 Recall that in the case of (approximately) massless fermions, the right-handed couplings to the 
W' vanish and those to the Z' are independent of et (c.f. chapter l4~3l . 
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the derealization parameter. This is no accident and can be understood by using 
the normalization of the fermion wavefunction to recast the expression for the Z'ff 
type couplings (14.30b as 

9z> f f,L = ±\ (gtf (fU) 2 + 9fi Z ' (f( L ) 2 ) + Yg'ff (7.1) 

where Y is the hypercharge and the sign depends on the isospin of the fermion 
(recall that the only source of a dependence of d7.ll) on is the wavefunction of the 
left-handed fermion f[). 

As the proton contains both up- and down quarks (albeit in the ratio 2 : 3) and 
9z'uu,l an d gz'dd,L both start with a positive value at the point of ideal derealiza- 
tion, the effect of tuning e l away from nonideal derealization partially cancels out 
from the Z' production cross section. In addition, there are also right-handed cou- 
plings of the same order of magnitude which don't depend on at all and there- 
fore, we shouldn't expect Z 1 production in the s channel to be very sensitive to ex,. 
W production, however, is completely forbidden in the case of ideal derealization 
and therefore highly sensitive to the value of ex- 

The partonic final state of the processes is reached through a cascade of the KK 
gauge boson going into two Standard Model bosons which then in turn decay into 
two fermion pairs. Excluding final states with four jets due to the QCD background 
and final states with more than one neutrino due to the missing momentum infor- 
mation, the following final states remain: 



w 


Z' 


jjlv 

jjll 

lllv 


jjlv 



Evidently, the only such state suitable for the detection of the Z' resonance consists 
of two jets, a lepton and a neutrino, while we can choose from three different final 
states in the case of the W, jjlv featuring the biggest branching ratio and lllv the 
smallest one. 

7.2 Simulation Details 

If we want to have any chance of observing the Z' in the s channel or to take advan- 
tage of all three final states suitable for the observation of the W' , we have to cope 
with final states containing one neutrino. Because the neutrino escapes the detector 
without leaving any traces, we have no direct information on its momentum which 
we require in order to measure the total invariant mass of the final state. Therefore, 
some trickery is necessary either to reconstruct the missing information or to live 
without it. 
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Neutrino mass shell condition —¥ 

W mass shell condition \ /■ 

Longitudinal neutrino momentum 

Figure 7.2: The two curves generated by the mass shell conditions for W and neu- 
trino in the case of a W decaying to lu. The two points of intersection give the two 
possible solutions for the neutrino momentum. 



The problem of missing momentum information is well known e.g. from SUSY 
decay chains where the transverse mass |60j is usually used instead of the invariant 
one in order to observe a resonance without requiring the full momentum informa- 
tion. However, our attempt to apply the transverse mass to the Z'/W' production 
processes under consideration was rather fruitless, the resonance in this observable 
being almost completely washed out. 

Luckily, this is not the only trick available for dealing with the missing neutrino 
momentum. Even though we cannot observe the particle itself, we still can infer 
the projection^ of its momentum on the transverse plane from transverse momen- 
tum conservation (assuming that the whole missing transverse momentum pT,miss 
originates from the neutrino). The mass shell conditions of neutrino and W bo- 
son (assuming that the neutrino originates from a W approximately on-shell) then 
give two additional conditions for the energy po and the longitudinal momentum 
component p i 

P 2 0-Pl-\P±\ 2 = (7.2a) 

2 

PoQo - PlQl - P±q± = (7.2b) 

(with the momentum q of the corresponding lepton and the W mass mw)- 

(|7.2a|l describes a hyperbola in the pl-Po plane, while (|7.2b|) describes a straight 
line as depicted in fig. 17.21 These curves can have up to two points of intersection, 
one of which corresponds to the correct momentum of the neutrino. The two solu- 
tions can be obtained analytically as 

gp { m w + 2 P±Q±) ± QlA q L (m w + 2p ± q ± ) ± A 

P0 ~ 2q (ql-ql) ' PL ~ 2 (gg - q\) 



CHAPTER 7. KK GAUGE BOSONS IN THE S CHANNEL 



87 



with the abbreviation 



A = q ^ (m 2 w + 2p ± q ± ) 2 + Ap\ (q 2 L - q 2 ) 



It is easy to see from (|7.2b|) that the modulus of the slope of the straight line is 
always smaller than 1 



Therefore, there will be always (excluding the rare case of straight line being a 
tangent of the hyperbola) either two points of intersection or none at all, and this 
method of reconstructing the neutrino momentum will always yield two solutions, 
only one of which will usually be close to the correct value. 

To cope with this, one can either try to find a kinematical criterion for discrimi- 
nating between the two solutions (at the price of losing part of the signal) or count 
them both into the histograms. Of these two approaches, we have chosen the lat- 
ter as it does not waste any precious signal events. The price to pay, however, is a 
doubling of the background. 

Performing the actual simulation reveals that for roughly 10% (a number quite 
sensitive to the W mass used in the reconstruction) of the events, the straight line 
doesn't intersect the hyperbola at all but instead lies slightly below the minimum, 
missing it by a small margin. Although we have chosen to simply remove these 
events from the histograms, more elaborate ways could be devised to deal with 
this, for example choosing the point on the hyperbola closest to the straight line 
or keeping mw as a free parameter which is determined for each event such that 
(|7.2a|) - (|7.2b|) have one degenerate solution and using it as a cut variable. A more 
in-depth discussion of this issue can be found in [63|. 

In order to determine the discovery potential of the LHC in these processes 
at parton level, we have performed a Monte-Carlo simulation similar to that pre- 
sented in the last chapter using our implementation of the model in WHIZARD 
/ O'Mega with y/s = 14TeV and an integrated luminosity of fC = 100 fb -1 . We 
have performed a complete simulation of the partonic processes pp — >■ X with X 
being one of jjlu, jjll and lllv (using the assignments d6.ll) and the same parton 
distributions as in chapter [6). 

Fig. 17.31 shows the general structure of the diagrams contribution to the reso- 
nance signal. As the final state fermions are decay products of Standard Model 
gauge bosons, we can leverage cuts on the invariant mass rriff of the visible fi- 
nal state fermion pairs in order to reduce the background^]. More specifically, we 

2 Note that after the momentum reconstruction, the invariant mass mi„ of the lepton-neutrino 
pairs is equal to mw per definition, rendering an invariant mass cut on mj„ pointless. 



Ql 



< 1 




while the asymptotes of the hyperbola (|7.2a|) have the slopes ±1 
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Figure 7.3: Left: General structure of the signal diagrams for Z' production in the s 
channel. Right: Dito, but W production. 

demanded rriff to lie within the ±5 GeV region around the invariant mass of the 
gauge boson which the fermions are to reconstruct: 



pp -> Z' -> jjlv 


75 GeV < mjj < 85 GeV 


pp — > W' — > jjlv 


86 GeV < mjj < 96 GeV 


pp — >• W — > jjll 


75 GeV < mjj < 85 GeV 
86 GeV < mu < 96 GeV 


pp^-W'^r lllv 


86 GeV < mu < 96 GeV 



Note that the narrow cut window of ±5 GeV is necessary in order to discrimi- 
nate between the W and Z' contributions to jjlv (as the masses of Z' and W are 
quasidegenerate, the resonance peaks fall together). However, it is quite unclear 
if ATLAS or CMS will be able to achieve enough precision in the measurement of 
the jet momenta for the separation to work out this way. This issue is discussed at 
length in section [731 

In the case of three leptons of the same generation in the final state, it is un- 
known which of the leptons originates from a decaying W. Therefore, we tried 
both possible combinations in this case and discarded the event if no unique iden- 
tification of Z and W could be achieved this way. 

In addition to the invariant mass cuts, we also applied px cuts of 

p T > 50 GeV 

on the momenta of all visible particles and also on pr.miss in order to suppress the 
background even more. Another cut of 

-0.95 < cos6> < 0.95 

was applied to the polar and intermediate angles of all visible particles, and as in 
the last chapter, we demanded 

E > 10 GeV 

for the incoming partons in order to avoid infrared divergences. 

For each final state, we performed simulations for the three parameter space 
points (16.211 , choosing three different values for e l from the interval allowed by 11351 
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Figure 7.4: Left: Invariant mass distribution in pp — > jjlv for different values of 
my/ 1 and M^uUs (Z' production). 

Right: The my/ 1 = 500 GeV peak obtained from the reconstructed neutrino mo- 
menta vs. the corresponding distribution obtained directly from Monte Carlo data. 

in addition to the point of ideal delocalization (c.f. fig. 17. II and 17. 11 1) for each value 
of my/ 1 . As in the case of W strahlung, the influence of (which was varied 

only for consistency) on the processes is much too small to be resolved, with the 
results depending virtually only on my/' arid the fermion delocalization e^. 

7.3 Z' Production 

Fig. 17.41 left shows the invariant mass distributions for pp — » Z' — » jjlv for my/i = 
380 GeV, 500 GeV and 600 GeV (again, recall that my/> ~ mz>) in the case of ideal 
delocalization. For all three values of the mass, the resonance peaks are clearly 
visible, their overall size declining with growing W mass as expected from the 
behavior of the parton distribution of the antiquark. 

In order to check the reliability of the neutrino momentum reconstruction, fig. 
17.41 right compares the my/i = 500 GeV resonance peak obtained by reconstructing 
the neutrino momenta to the distributions created from the exact momenta taken 
directly from the eventgenerator. As expected, the reconstruction roughly doubles 
the number of background events, while the size of the peak remains more or less 
the same. However, the peak in the reconstructed data exhibits broad "sidebands" 
which closer investigation reveals to originate from the second solutions to the neu- 
trino momentum in events within the peak. 

The histograms in fig. 17.41 have been generated for the case of ideal delocaliza- 
tion in spite of the claim made in [35 J (and discussed in chapter |4~2)| that this case 
is in fact excluded by electroweak precision data. Instead, e l must be tuned away 
from ideal delocalization in order to allow for a small but nonvanishing W'ff cou- 
pling. The resulting variation of the my/' = 500 GeV resonance is shown in fig. 17.51 
left for three different cl chosen from the allowed interval (c.f. fig. l7.1|) . 

As argued in 17.11 this is rather small and causes a small increase of the event 
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Figure 7.5: Left: The effect of tuning sl away from ideal derealization (cf. fig. I7.ll 
ITTA on the Z' peak. 

Right: The significance as defined in the text as a function of the integrated lumi- 
nosity. The dotted lines mark the 3a resp. 5a discovery thresholds. 

count in the peak as compared to the case of ideal derealization. Therefore, we 
will confine ourselves to presenting the results for this scenario, with the effect 
caused by varying cl being small and pushing in the direction of slightly larger 
significance. 

To get a quantitative handle on the significance of the signal and to estimate the 
minimal luminosity necessary for discovering the Z', let's proceed similarly to the 
analysis in chapter I7.3I and define the raw signal N to be the number of events in 
the ±20 GeV region around the peak. In order to estimate the background, we have 
generated SM events for an integrated luminosity of JC = 400 fb -1 , analyzed this 
data the in the same way as the Monte-Carlo data for the Three-Site Model and then 
downscaled the resulting distributions by a factor of 4 to reduce the error coming 
from fluctuations in the background. We denote the number of background events 
in the ±20 GeV region around the peak obtained this way by N b . 

We then again define the signal N s as 

N S = N-N b 

The determination of the standard deviation of the background is complicated by 
the fact that the pairs of entries in the histograms for the different solutions to the 
neutrino momentum lead to a nontrivial statistical correlation in the histograms. 
For simplicity, let's assume that the background is simply doubled by the momen- 
tum reconstruction 

N b = 2N' b 

The standard deviation of a^ b of N b must scale accordingly 

resulting in the statistical significance defined as in chapter 1631 being reduced by a 



CHAPTER 7. KK GAUGE BOSONS IN THE S CHANNEL 



91 



factor of y/2 

s _ Ns_ _ N-N b 

which nevertheless scales like \J~Jc~ 

The significance of the signal in the ideally delocalized scenario thus calculated 
is shown in figure [731 right as a function of the integrated luminosity together with 
the 5a and 3d discovery thresholds. The 5a thresholds are approximately lfb -1 , 
2 fbr 1 , 5 fb- 1 for m w > = 380 GeV, 500 GeV, 600 GeV, respectively. Considering the 
fact that tuning ex into the region allowed by the precision observables does not 
significantly change the signal, the Three-Site Z' may be discovered in this process 
as early as in the first 1 — 2fb _1 and even in the worst case can be expected to 
manifest itself in the first 10 — 20 fb _1 of data. 



7.4 W Production 

In the case of W production we have the three different final states jjlu, jjll and 
lllu which we can use to look for the resonance peak. The left column of fig. 17.61 
shows the resulting invariant mass distributions for these final states, again for all 
three my/' values \6.2\ under consideration with ex chosen to yield roughly the 
maximum allowed value of gw'ff ( c -f- fig- 17. II and fig. 17.11)1 . The peaks are clearly 
visible in all three final states, with the total number of events contained in them 
dropping in the order (biggest to smallest) 

jjlu — > jjll — > lllv 

due to the declining branching ratio (as can be easily seen by simply counting the 
number of final states). 

However, going from jjlu to jjll, the need for the reconstruction of the neutrino 
momentum with the associated doubling of the background events also goes away, 
and the possibility of cutting on the invariant mass of the W enhances the ratio 
of signal to background. Furthermore, the hadronic background is completely re- 
moved when going to the fully leptonic final state lllu, leading to the much cleaner 
signal visible in fig. I7.6l in spite of the even smaller branching ratio. 

Contrary to the Z 1 case, the signal in the W production process must be pro- 
portional to the square of the gy/' / / coupling and is therefore highly sensitive to 
the fermion derealization parameter ex- The corresponding dramatic change in 
the my/' = 500 GeV resonance induced by changing ex within the allowed part of 
parameter space is shown in the right column of fig. l7.6[ the change in the peaks for 
the other W' masses in fig. 17.101 From these distributions, it is evident that the peak 
can be nearly completely removed within the allowed region of parameter space 
by choosing a suitable value for ex- 

As would be expected, the resonance completely vanishes in the (albeit forbid- 
den) case of ideal derealization. The histograms show that this is also true for the 
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Figure 7.6: Left column: The W resonance peaks for the three different final states 
under consideration. Different values of mw', near-maximal gwff- 
Right column: Close-up of the my/ 1 = 500 GeV resonance peak for different values 
of £l from the interval allowed by the precision observables. 
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Figure 7.7: The integrated luminosities required for a 3a resp. 5a discovery of the 
W' in the s channel as a function of el- Different final states, my/' = 500 GeV. 



jjlu final state(with the exception of my/' = 380 GeV in which case a small peak 
from misidentified Z' events remains), demonstrating that the ±5 GeV cut around 
the W mass is sufficient to separate W and Z' contributions to this final state at 
least at parton level if the smearing cause by measurement errors is not take into 
account. 

Defining the significance of the signal in exactly the same wa >| as in [73] for Z' 
production (without the additional factor v2 for jjll), we can exploit the quadratic 
dependence of the signal on gy/i //in order to obtain an estimate of the integrated 
luminosity required to yield a signal with a given significance. The ha resp. 3a 
result is shown in fig. 17.71 as a function of for my/> = 500 GeV; the corresponding 
plots for the other two masses are shown in fig. 17.121 

From these curves it is clear that, while the actual result varies a bit over param- 
eter space, the performance of the different final states is fairly comparable at the 
parton level. However, once detector effects are included, it is likely that the fully 
leptonic final state is preferred due to the more precise measurement of the lep- 
ton momentum and the absence of complicated QCD backgrounds — an assertion 
which is confirmed by the detector simulation performed in [63 J. 

As far as the prospects of discovering the W' in this process are concerned, fig. 
17.71 and 17.121 clearly demonstrate that there are regions in parameter space easily 
accessible in the first 10 — 20 fb _1 , while it would take more than 500 fb _1 for a 5a 
discovery in other regions which is most likely out of question at the LHC. 



3 As the total number of events is much lower in the jjll resp. lllv final states when compared to 
jjlu, we generated the Standard Model backround for these by downscaling simulations for JC — 
1000 fb" 1 resp. f£ = 5000 fb _1 . 
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7.5 Disentangling the jjlv final states 

Since flavor tagging is impossible for light quark flavors, we have to rely on invari- 
ant mass cuts on the jet pairs in order to be able to separate the case of the two 
jets in jjlv coming from the decay of a W in Z' production from that of the jets 
being produced by a decaying Z in W production. However, it may very well be 
impossible to obtain a resolution of order ±5 GeV in the jet invariant mass from 
experimental data. In the following, we model the effect of the measurement error 
on the W'/Z' separation with a gaussian smearing of the invariant mass of the jets. 

In the ideal case of exact rrijj measurement, events coming from the decay of a 
intermediary W / Z are distributed according to a Breit-Wigner distribution 

Tih(?n r}~ 

Pb(x, to, T) dx = ^ dx 

(x 2 — TO 2 ) + T 2 TO 2 

with the normalization factor 



,2 



( r \ 3 -if 1 r 

n b (m, T) = - — ^ 1 H ^ sin - atan — 

4to j y to J \2 TO 

Emulating the measurement error in the jet mass by convoluting pt>w with a gaus- 
sian of standard deviation a 

1 _ z 2 

vAx, a) dx = , e 5^ dx 

8 V2^a 

we obtain the smeared distribution 

p sm (x,m,T,a) dx = / dy p b (y,m,T)p g (x - y,a) 
Jo 



Figure 17.81 shows the effect of this smearing on the Breit-Wigner peaks of the 
Z and the W. Turning on the smearing and increasing a causes the sharp Breit- 
Wigner peaks to decay rapidly, and for a = 10 GeV, only two very broad bumps 
are left. The consequence is that, if a cross section has one contribution which stems 
from the decays of a virtual Z and one coming from a virtual W, any attempt to iso- 
late the Z contribution by cutting on the resonance will inevitably also select events 
coming from the W decay contaminating the sample (and vice versa). Therefore, 
our analysis of the jjlv final state will show a W peak even in the case of ideal 
derealization which is caused by jet pairs from a decaying W misidentified as a Z. 

If we try to isolate the W peak with a cut on the invariant mass m,jj 

Lw < m jj <Uw 

and the Z peak with a cut 

L z < mjj < U z 
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Figure 7.8: The effect of a gaussian smearing on the Breit-Wigner shape of the W 
and Z resonances for various widths a of the gaussian. 



then the resulting event counts Nw, N z can be calculated from the true event counts 
Nw, Nz coming from a decaying W or Z via a matrix T as 



Nw \ _ ( Tww T wz \ ( N w 
N z j \T ZW T ZZ )\N Z 



with entries 



Tij = dmp sm {m,mj,Vj,a) 

JLi 

Inverting T we can calculate the event counts Nw and N z 

■ T- 1 



N w 
N z 



Nw 
Nz. 



(7.3) 



The entries of T give the probability of (mis)identifying an event and can be readily 
calculated numerically; for example, choosing cuts 



L w = 60 GeV 



U w = 85 GeV 



yields 



0.64 
0.29 



0.27 
0.62 



L z = 86 GeV 



1.9 

-0.89 



U z = 111 GeV 



-0.85 
2.0 



This way we can in principle use T to disentangle the contributions from the W and 
Z resonances to the signal in the presence of a measurement error which causes the 
Breit-Wigner peaks to lose their shape. However, in order to apply this to actual 
data, it is vital to separate the signal from both the reducible and the irreducible 
backgrounds as these don't obey a Breit-Wigner distribution. 
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Figure 7.9: Left: Signal in the W detection channel for the case of ideal derealiza- 
tion smeared with a gaussian error. 
Right: The same for the case of nonzero gw' / / 



In order to estimate the significance of a signal obtained this way, we calculate 
the standard deviation of N{ according to 

Calculating the fluctuations of Nj is a bit tricky as the neutrino momentum recon- 
struction adds statistically correlated pairs of events to the sample. According to 
the analysis performed in the last two sections 17.31 and [7!4l we obtain the number of 
events N in the peak by subtracting the background (aka Standard Model expecta- 
tion) N b from the total number of events iV tot 

N = N tot - N b 

The Standard Model value N b comes from a theoretical calculation and is by 
definition free of statistical fluctuations^. The total number of events AT tot contains 
N events originating from the resonance which we idealize to be free of double- 
counting and N b background events which contain the correlated double counts. 
Putting the pieces together, we have 

°N = a N tot = \/N + 2N b = y/N toi + N b 

which we can now plug into (17.4b in order to calculate the statistical fluctuation we 
have to expect for the reconstructed event counts and finally arrive at 

**i = . / E C^ 1 ) 2 {Nt >i + N »J (7 - 5) 

V itw,z 

4 Of course, this is a simplification which is not quite true in our analysis where we obtain N b from 
a Monte-Carlo simulation. However, due to our generating and downscaling a surplus of events, the 
error on N b is smaller than that on N io t, and in a more serious analysis of experimental data, <Tjv b 
could be made arbitrarily small by investing enough time into its determination. 
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ideal derealization 




Ni 


Si 




Nj 


i = W 


3193 


17 


5126 


13 


i = Z 


1371 


7.5 


-96.10 


0.24 



e L = 0.254 




Ni 


Si 


Ni 


Ni 

ONi 


i = W 


3767 


21 


5628 


14 


i = Z 


2083 


11 


811.6 


2.0 



Table 7.1: Comparison of the signals N w /z obtained with an gaussian smearing of 
the invariant mass of the jets with a = 10 GeV to the "true" signals N w / Z calculated 
from the measured ones via the transfer matrix T -1 . 



For a simulation of the effect of the measurement error our analysis we have 
randomly distributed the invariant mass of the jet pairs within a gaussian with 
width a = 10 GeV centered around the correct value calculated from Monte Carlo 
data. We then performed the same analysis as in the last two sections 17.31 and 17.41 
with mw' = 500 GeV and m,b u ik = 3.5 TeV both for = 0.254 and for the ideally 
delocalized scenario. The only difference to the previous analysis are the cuts on 
rrijj which we enlarged to 

60 GeV < rrijj < 85 GeV resp. 86 GeV < m 3j < 111 GeV 

Fig. 17.81 shows the resulting effect on the W peak for the cases of ideal dereal- 
ization (left) and for €l = 0.254 (right). In both cases a peak is clearly visible, which 
in the ideally delocalized scenario is only composed of events with jets coming 
from a decaying W misidentified as a Z. 

The number of signal events N w / Z after smearing, the significance s w / z of these 
calculated similarly to the last sections, N w / Z obtained from applying the trans- 
fer matrix T _1 A7.3D and the resulting significance Ni/a^. obtained from A7.5D are 
shown in table 17.11 All peaks are significant with s > 5er; however, after applying 
the transfer matrix, the W peak vanishes within one standard deviation for ideal 
delocalizatiord, while in the case of ei = 0.254 a residue as big as 2a remains. The 
Z' peak remains significant after applying the transfer matrix although, however, 
the significance is reduced because the transfer matrix enlarges the error. 

7.6 Conclusions 

What is the message to be taken from the simulations presented in this chapter? 
At the beginning, we have argued that W'/Z' production in the s channel is a very 
sensitive probe of the fermion sector in the Three-Site Model and in particular of 
the delicate interplay of fermion and KK gauge boson wavefunction which allows 
the model to evade the precision constraints. 

In the simulations presented here we have shown that the heavy gauge bosons, 
albeit fermiophobic, can indeed lead to observable resonances in the s channel 

5 The negative event count in the W' peak is an artifact of the multiplication with T _1 . 
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at the LHC. In the case of the Z' , we have found that the corresponding reso- 
nance shows only little dependence on the fermion derealization parameter e l and 
should be expected to be visible within the first 10 — 20 fb _1 in the jjlv channel for 
any point in parameter space. At the LHC, such a signal together with the absence 
(or near-absence) of a corresponding signal in the dilepton and Drell-Yan channels 
would be a strong sign for a fermiophobic, heavy neutral vector resonance like the 
Three-Site Z' . 

As far as the W' is concerned, the simulations have also shown that the res- 
onance might be accessible at the LHC. However, the magnitude of the signal is 
highly dependent on cl, and there are regions in parameter space in which no sig- 
nal would be observed at the LHC. Even in this case though, the existence of the W 
could be established via the strahlung process presented in the last chapter (which 
is essentially independent of cl), and the absence of a corresponding peak in the s 
channel processes would be a clear sign of its fermiophobic nature. 

We have found three different possible discovery channels for the W, all of 
which show a comparable discovery potential at the parton level. However, the 
jjlv type final states contain contributions from both W' and Z' , and if ATLAS and 
CMS turn out to be unable to resolve the jet momenta with the accuracy required 
for separating the resonances via cuts on the invariant jet mass, some trick like the 
one presented in 17.51 needs to be applied if this final state is to be exploited for W' 
production. Otherwise, only the combined resonance of Z' and W can be observed 
in this channel which, luckily, at least wouldn't hurt the Z' detection much as the 
signal is stronger than the W one anyways, and the amount of W events contained 
in the peak could be inferred from the other two W final states available. 

All simulations presented in this thesis have been performed at the parton level 
only with no detector effects taken into account. In the case of those presented in 
this chapter, an ATLAS detector simulation was added in the master's thesis of F. 
Bach [63 J. This also includes additional backgrounds which were purposefully not 
taken into account in the parton level simulations as they depend on details of the 
detector and of the jet identification algorithm, e.g. the background coming from tt 
production with two jets escaping down beam pipe. 

The result shows that the Z' signature is rather robust and persists with a 5<r 
effect still being visible in the first 5 — 30 fb^ 1 even if detector effects are taken into 
account. The W resonance turns out to be more fragile, with the separation trick 
for jjlv not working well after the detector simulation and the signal in jjll be- 
ing diminished by the measurement error on the jet momenta and additional back- 
grounds. In light of this result, lllv becomes the preferred final state for this process 
at the LHC, and the detector simulation suggests that the parton level results carry 
over to the actual experiment in this case. For more details on these issues see [63J. 
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7.7 Additional plots 
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Figure 7.10: The W resonance peaks not covered by fig. 17.61 for the different final 
states under consideration. Different values of mw> and ei- 
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7.7 ADDITIONAL PLOTS 





Figure 7.12: The same as fig.EZlfor m w > = 380 GeV, 600 GeV. 



Chapter 8 

Heavy Fermion Production 



Oliphaunt am I, 
Biggest of all. 

(Excerpt from the "Oliphaunt song" in J.R.R. Tolkien's "Lord of 

the Rings") 

After discussing at length the prospects of discovering the KK partners of the Z 
and W as well as the possibility of getting a handle on their fermiophobic nature, 
the only aspect of the model which we have not addressed in our simulations so far 
are the heavy fermion partners. With masses above 1.8 TeV and widths of several 
100 GeV (c.f. chapter HJ, these particles certainly are near the limit of what may be 
accessible as a resonance at the LHC. Nevertheless, we will see in this chapter that 
this may be feasible at least in some portions of parameter space. 

As in chapters [6] and we will first discuss the different processes that could 
be exploited for probing for this kind of new physics in 18.11 Subsequently the 
details of the simulations we performed will be discussed in 18.21 while the results 
are presented and discussed in 18.31 

8.1 Processes 

Besides the fermiophobic W and Z' , the main feature of the Three-Site Model are 
the new heavy fermion resonances. As we have seen in chapter [21 their introduc- 
tion was required in the Three-Site setup in order to allow for the tuning of the 
couplings between the Standard Model fermions and the KK gauge bosons which 
in turn is necessary for evading the electro weak precision constraints. Therefore, 
regardless of their large mass scale Mb u ik > 1.8 TeV, these particles are an inte- 
gral part of the phenomenology of the Three-Site Model, and if other experimental 
results would hint at such a kind of new physics, a way to access these particles 
would be desirable. 
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Figure 8.1: Total cross section for pp — > jq' as a function of M^^. 

The production mechanism which might first come to mind are simple Drell- 
Yan like processes of the type 




However, such processes always involve a quark-antiquark pair in the initial state 
which is very unlikely to deliver the large amount of energy necessary for the cre- 
ation of a heavy fermion at a pp collider like LHC, and simulating pp — > e + e~' with 
WHIZARD / O'Mega indeed reveals a cross section of the order of 0.5 ab at the 
LHC. Clearly, even if we multiply this with with an optimistic factor of 42 in order 
to take all Standard Model fermions visible to the detector into account the result- 
ing cross section is still pathetic, and detecting the heavy fermions this way is out 
of question. 

However, there is another kind of potential production processes which are me- 
diated by t channel interactions 




This kind of process can contain two valence quarks in the initial state with much 
more energy being typically available in the partonic CMS than in the above case 
of a valence and a sea quark. The resulting cross section for pp — > jq' (q' being any 

1 (6 quarks x 3 colors + 3 leptons) x 2 (fermions + antifermions) = 42 
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Figure 8.2: Signal diagrams for t channel induced production of heavy fermions 
with the two final states as obtained from the different decay channels of the heavy 
fermions. 

of vf, u', df, d!) is shown in fig. l8.ll as a function of Mt, u ik, exploiting the fact that the 
heavy partners of the light Standard Model fermions are degenerate (c.f. chapter 
I4.1|) by summing over fermion flavors in order to enhance the signal. Evidently, the 
available energy is much higher, and the resulting cross section can be as big as r* 
66 fb for the smallest Afbulk = 1-8 TeV, dropping only moderately with nearly 10 fb 
being still available for M^uxk = 3 TeV. These numbers are much more encouraging 
than those in the Drell-Yan like case, and we can indeed hope that it might be 
possible to dig out a heavy quark resonance this way if is not too high. 

As discussed in chapter 14.41 the heavy fermions decay into a Standard Model 
fermion and a gauge boson, the ratio between decays into light and heavy gauge 
bosons being roughly 1:1. The Standard Model W and Z then decay into two fermi- 
ons, leading to four particle final states like fig. l8.1l left, while the W and Z' cascade 
into four fermions by means of two intermediary gauge bosons, resulting in a six 
particle final state similar to fig. l8.1l right. 

The six particle final states come with the advantage that a cut on the invariant 
mass of the heavy gauge boson can be leveraged in order to reduce the number of 
background events (given that the W'/Z' mass has been previously determined). 
However, a preliminary study of the six particle final state containing four jets 
turned out to be quite involved due to the large number of background processes 
already present at the parton level, and the resulting event count in the signal chan- 
nels was not big enough to leave much hope for a significant discovery potential 
in this process. The investigation of six particle final states of the form jjlllv is still 
an ongoing project, but as processes leading to this final state are suppressed rela- 
tive to those with four particles in the final state by another gauge boson branching 
factor, it is unlikely that they will offer any advantage over these. 

Therefore, we only present results for the four fermion processes in this work. 
Excluding the case of four jets due to the large QCD backgrounds, we remain with 
jjlv and jjll as potential final states. However, as the branching ratio into Z — > 
11 is significantly smaller than that of W — > lv, we should expect the jjlv final 
state to perform much better than j jll, provided that we can cope with the missing 
neutrino momentum. 

Apart from the mass and couplings of the heavy fermions, the prospects for 
discovering these particles as a resonance also depend on their total width. As the 
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Figure 8.3: The total u' width as a function of for the minimal value of the 

W mass mw' = 380 GeV and for the highest mw' allowed by the constraint aT < 
0.005. 

amount of phasespace available for the decay of a particle depends on the mass 
of the decay products, we should expect that the total width exhibits a strong de- 
pendence on mw'- Fig. I8.3I shows the total v! width as a function of Mbuik both 
for m\v' = 380 GeV (which is the smallest value allowed by the LEP constraint on 
the triple gauge boson couplings) and for the highest W mass which is still in ac- 
cordance with the condition aT < 0.005 (c.f. chapter |4~2)| . This plot demonstrates 
that the range of potential values for the width grows significantly with increasing 
Mmlk/ an d for Mb u ik = 3 TeV, the difference is more than 200 GeV. Clearly this 
variation of the width with m-w is big enough to have an influence on the signifi- 
cance of the resonance. 

An additional dependence on m\y comes from the Wf'f / Zf'f couplings 
which scale with x (c.f. chapter I4.3)> . Looking at (|4.19|) , it is obvious that, while 
the resonance becomes more narrow and easier to detect with growing W mass, 
the couplings responsible for the production of the heavy quark in fact go down, 
and it is not clear a priori which effect dominates the significance of the signal in 
the end. 

In contrast to the s channel production of heavy gauge bosons discussed in the 
last chapter, we have no reason to expect a large influence of the derealization 
parameter on this process. The only possible source of such a dependence ore 
processes involving a W in the t channel, the relative contribution to the amplitude 
of which would be 

9W'qq',L 9Wqq,L 
9Wqq',L 9Wqq,L 

(the right-handed components of massless fermions sit completely at the right brane 
with their partners residing in the bulk and therefore, the right-handed couplings 
gw'q'q.R vanish). This turns out to be less than 10% over the whole parameter space 
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and therefore, we are free to perform the simulation in the ideally delocalized sce- 
nario with the effect of nonideal derealization on the result being small and negli- 
gible when compared to other uncertainties. 



8.2 Simulation Details 



In order to examine whether it might indeed be possible to discover the heavy fer- 
mions at the LHC, we have performed full parton level simulations of the processes 
PP — > jjll an d pp — > jjlv for an integrated luminosity of JC = 400 fb -1 and a total 
energy of 14 TeV, again with the definitions d6.ll) and the same parton distributions 
as in chapters [6] and In the case of the final state containing a neutrino, we have 
applied the reconstruction method already used in chapter [7| and presented in !7.2[ 
again double counting the event if both solutions pass the cuts. To get a sufficiently 
smooth estimate of the background, additional simulations in the Standard Model 
were performed for an integrated luminosity of fC = 1600 fb _1 and the result then 
downscaled by a factor of 4. 

In order to suppress the backgrounds and identify the jet which originates from 
the heavy fermion, we concentrate on events in which the total energy available in 
the parton CMS is just barely enough to create the KK particle which is therefore 
only weakly boosted in this frame. In the CMS, this implies that the jet coming from 
the decaying heavy quark is emitted isotropically in all directions and is approxi- 
mately back-to-back to the W/Z coming from the decay, while the jet originating 
from the proton will be aligned with the beam axis. Note that these kinematics are 
very different from those exhibited by gauge boson fusion processes like 




which feature two high rapidity jets aligned with the beam axis and which can be 
expected to constitute the largest part of the irreducible background. 

To facilitate the selection of events whose kinematics match these criteria we 
applied the following cuts in the parton CMS 
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polar angle of jet j f originating from 

J. V V Vj LI 1 IX 


< |cos0j/| < 0.7 


polar angle of jet j originating from pro- 
ton 


0.95 < (cos^l < 1 


angle between j' and W/Z 


-1 < cosOjw < -0.8 


Pt of heavy quark 


< Pr,q> < 200 GeV 


Pt of f and W 


500 GeV < p TJ i , Pt,w 



We also applied a pr cut on all observable momenta (and also on pT,miss i n the case 
of pp -> jjlu) 

p T > 50 GeV 



as well as a cut on the total invariant mass 

mtot > 1 TeV 

for additional background suppression. In order to avoid infrared divergences and 
for simulating the acceptance region of the detector, we demanded 

-0.99 < cos 6* < 0.99 

for the polar and intermediate angles of all observable particles in the hadron CMS, 
and the same small x cut already used in the last two chapters [7|-[8] 

x > 1.4- 10~ 4 

was applied on the incoming partons (also in the hadron CMS). In the case of the 
jjll final state, we also applied an identification cut on the invariant mass of the 
dilepton system 

81 GeV < mu < 101 GeV 

In order to identify the part of parameter space in which the heavy quarks might 
be visible at the LHC, we have conducted simulations of the jjlu final state at ten 
different points in parameter space as shown in tab. 18.11 For each value of M\^, a 
simulation at mw' = 380 GeV (which is the lower limit given by the LEP / LEP-II 
constraint on the triple gauge boson couplings, c.f. chapter I4.2[) was performed. 
In addition, for M^ n \^ > 2 TeV, another simulation near the upper limit on mw 
(see fig. I8.3|) was conducted in order to probe the dependence of the result on the 
W mass. As the simulations of pp — > jjll at the points I and II already show that 
jjlu yields a much better signal, no simulations of jjll were performed at Illa/b - 
Vla/b. 
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parameters 


heavy quark properties 


point label 


M bulk [TeV] 


mw> [GeV] 


m [TeV] 


r [GeV] 




1.8 


380 


1.90 


150 


I 


2.0 


380 


2.10 


205 


II 


2.2 


380 


2.31 


273 


Ilia 


430 


2.29 


204 


im? 


2.5 


380 


2.63 


401 


IVa 


460 


2.58 


257 


IVb 


2.7 


380 


2.84 


506 


Va 


480 


2.78 


294 


Vb 


2.9 


380 


3.05 


628 


Via 


500 


2.98 


333 


Vlb 



Table 8.1: The parameter space points at which the simulations of heavy quark 
production were performed. The first two columns define the parameter space 
point (we always assume ideal derealization, c.f. section I87TT) , the second two give 
mass and width of the KK partners to the light Standard Model quarks, and the 
labels used to identify the points in the analysis are given in the last column. 




Figure 8.4: The heavy quark resonance in the jjlu final state for the parameter 
space points I (left ) and II (right), visible as an excess over the Standard Model back- 
ground. 
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8.3 SIMULATION RESULTS 



S o 



final state: iilv 



parameters: I 
parameters: II 



2 3 
j' + W invariant mass [TeV] 




parameters: I 
parameters: II 



2 3 
j' + W invariant mass [TeV] 



Figure 8.5: Left: The excess in the j' + W invariant mass with the heavy quark 
resonance. Right: The per-bin significance of the excess calculated as described in 
the text. 

8.3 Simulation Results 



Let's start our discussion of the simulation results with the jjlv final state. Fig 
shows the invariant mass of the combination of the jet j' potentially originating 
from a heavy quark (as identified by our cuts, see the preceding section) and the 
W for the parameter space points I and II (see tab. 18.1b . While the background is 
considerable, an excess in the histogram at the invariant mass of the heavy quark 
over the Standard Model background is clearly visible (the peak in the distribution 
around 1.2 TeV is not related to the new physics, but a kinematical artifact induced 
by the cuts). 

Calculating the difference between the event counts in the Three-Site Model 
and those in the Standard Model as shown in fig. 18.51 left for points I and II, we 
find that the excess is quite big, containing about 1000 events, and indeed has the 
shape of a broad resonance peak. In order to calculate its significance, we follow the 
same reasoning as in chapter 17.31 and assume that the double counting performed 
in the neutrino momentum reconstruction step simply doubles the event count and 
define the significance of the excess by 



Si 



(8.1) 



where the N 3 i are the event counts in the Three-Site Model and the N b < are the 
counts in the Standard Model. The resulting significance is shown on a per-bin 
basis in fig. 18.51 right, demonstrating the excess in the bins lying within the peak 
to be significant with s > 10. In addition, while the "naked" difference shown in 
fig. 18.51 left also exhibits a noticeable difference in the leftmost bins around 1 TeV, 
fig. l8.5l right reveals this to be not significant and in fact just a statistical fluctuation, 
the conspicuous magnitude of which is solely due to the number of background 
events in this energy range. 

In order to get a quantitative estimate of the significance of the whole resonance, 
we define the signal N s to be the number of events within the ±r range around 
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total integrated luminosity 

Figure 8.6: The significance of the heavy quark resonance as defined in the text at 
the points I and II as a function of the integrated luminosity together with the 3a 
and 5a discovery thresholds. 

the invariant mass^l of the heavy quarks. From this, we calculate the significance 
according to (|8.1|) just as we did for the individual bins. The result for the points I 
and II is shown in fig. I8.6I as a function of the integrated luminosity together with 
the 3a and 5a discovery thresholds. At both points, the signal exhibits nearly the 
same significance of « 15a for JC = 400 fb^ 1 and intersects the 5a threshold already 
around 45 fb _1 . 

The upper two rows of fig. l8.7l show the difference in the j' + W invariant mass 
between the Three-Site Model and the Standard Model for the parameter space 
points III/IV/VIa and III/IV/VIb (Va/b has been excluded in order to make the 
figures more legible) with the left column giving the raw event count (similar to 
fig. 18.51 left) and the right column again the significance of the excess (like fig. 18.56 . 
For both the points with minimal and maximal mw, the resonance remains visible 
while broadening as it moves to higher invariant mass values. 

Similar to fig. 18.61 the bottom row of fig. 18.71 shows the significance of the res- 
onance peaks as a function of fC for the points Illa/b - Vla/b. While the signifi- 
cance crosses the 5a threshold below 400 fb _1 for Illa/b - Va/b, the significance is 
just barely 5a for Via and even below that for VIb at 400 fb _1 . For nearly all values 
of ^4,^, the significance is worse for the maximum value of mw even though the 
resonances are considerably broader for mw = 380 GeV (c.f. tab. I8.1|) . The only 
exceptions are Va/b (Mbuii- = 2700) for which Vb gives a slightly better result than 
Va. However, as the determination of the significance comes with an uncertainty 
of about \a at fC = 400 fb _1 , this may well be a statistical fluctuation. 

Let us now move on to the jjll final state. For this process, fig. 18.41 shows the 

2 Choosing a wider interval gives a slightly higher significance of the signal at parameter space 
points with low r, but actually turns out degrade the signal at points with high width like Va and 
Via. 
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8.3 SIMULATION RESULTS 



linal state: jjlv 
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parameters: IVa 
parameters: Via 
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parameters: IVb 
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Figure 8.7: Top row: Similar to fig. 18.51 left, but for the points III/IV/VIa and 
III/IV/VIb. Middle row: Similar to fig. 1831 right, but for III/IV/VIa/b. Bottom row: 
Like fig. EE but for Illa/b - Vla/b. 
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f + W invariant mass in the Three-Site Model versus the Standard Model expecta- 
tion in a fashion similar to fig. l8.8l for the parameter space points I and II. Again, the 
resonance is clearly visible as an excess over the Standard Model background, but 
the total event count is much (more than an order of magnitude) lower than in the 
jjlv case. The corresponding difference of the Three-Site Model and the Standard 
Model is shown in fig. l8.9l left, making the resonance peak manifest, but the fluctu- 
ations are much higher than in the jjlv case due to the lower total event count. 

This is emphasized by the significance of the excess showcased in fig. 18.91 right 
(with the significance again being calculated via d8.ll) , but now without the factor 
y/2 in the denominator due to the absence of double counting). Already in this 
histogram, the significance of the signal turns out to be much lower than that ob- 
served in fig. l8.5l right for the jjlv case. The total significance of the signal is shown 
in fig. 18.101 as a function of fC. For 400 fb -1 , it is only about « 6.5a at point I and 
even below 5a for II, which is clearly much worse than what we observed for jjlv 
in fig. 18.61 As evidently jjlv is much better suited for the detection of the heavy 
quarks in this process, with a 5a discovery already being difficult for M^ u ^ = 2 TeV 
in the jjll final state, we didn't perform any more simulations for pp — > jjll. 
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8.3 SIMULATION RESULTS 
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Standard Model 



parameters: I 
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parameters: II, 14 bins. 
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Figure 8.8: Like fig. l8.4l but for the jjll final state. 
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Figure 8.9: Like fig. I8.51 but for the jjll final state. 
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Figure 8.10: Like fig. 18.61 but for the jjll final state. 
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8.4 Conclusions 

While the detection of the heavy fermions at the LHC is challenging due to their 
large mass and considerable width, our parton level results show that this should 
indeed be possible in the jjlv (and to a lesser extend in jjll) final states at least 
in the domain of the parameter space with M^ u \^ smaller than 2.9 TeV. While the 
signal significance depends heavily on Mb u ik, the dependence on mw> turns out 
to be less than might be expected form the associated change in the width of the 
resonance which indeed turns out to be (over)compensated by the decrease of the 
relevant couplings. In addition, as argued in the beginning of this chapter, the 
effect of changing the derealization parameter cl should be small and not more 
than 20% in the worst case. 

However, even in the accessible part of parameter space, several caveats remain. 
First, this kind of processes is restricted to the heavy quarks, and no information 
on the heavy leptons or neutrinos can be obtained this way. Second, due to the 
near-degeneracy of the heavy fermions masses, only the heavy quarks as a whole 
(and, if small contributions from flavor mixing are ignored, of these only the u' 
and d!) are accessible. For a resolution of individual flavors, tagging would be 
required which is impossible for the light flavors. Third, the background is rather 
large, reducing the broad resonances to an excess in the histograms, from which 
the Standard Model background must be subtracted in order to be able to resolve 
the resonance, and it is not clear a priori whether a sufficiently good background 
estimate will be available. 

On the pro side, we can expect that the inclusion of hadronization and detector 
effects into our parton-level results would not harm the outcome overmuch. The 
main effect of these effects would supposedly be a smearing of the invariant mass 
distribution due to the finite jet resolution, but as the resonances are very broad 
anyway, no big harm should come from that. More dangerous is the potential in- 
crease in the reducible background coming from other final states which we did 
not take into account in our parton level analysis; this should be studied carefully 
if such an analysis of the experimental date should ever be performed. 

So, can the heavy fermions be seen at the LHC? In the light of the results of 
our simulations, the answer to this question is a definitive "maybe", with the final 
verdict depending on the actual value of Mb u ik, the accuracy of our understand- 
ing of the Standard Model background, the amount of data the LHC will be able 
to collect and, also important, the actual energy the collider will reach — our sim- 
ulations were performed for 14 TeV, and less energy would diminish the result 
considerably. 



114 8.4 CONCLUSIONS 



Chapter 9 

Conclusions 



The physics is theoretical, but the fun is real. 

(Tagline of Sheldon's "Research Lab" game from "The Big 
Bang Theory — The Guitarist Amplification") 

While it certainly does not solve most of the issues the Standard Model is usually 
criticized for (apart from the hierarchy problem which obviously is void in Higgs- 
less model), we have argued in chapter [2] that a scenario like the Three-Site Model 
does arise in a rather natural fashion once we try to replace the Higgs with a set of 
SU(2) vector bosons as agents for the preservation of perturbative unitarity at the 
TeV scale. Discussing the properties of the new structure in the model and the con- 
sequences of the experimental constraints in chapter HI we have found that these 
constraints imply that the new physics is well hidden from discovery by forcing 
the W'/Z' to be essentially fermiophobic and the KK partners of the fermions to be 
very heavy. 

The goal of this thesis was the implementation of the model into the WHIZARD 
/ O'Mega framework and the application of this implementation to studying the 
phenomenology exhibited by the model at the LHC. The implementation as pre- 
sented in chapter [5] is complete and has been cross-checked against two indepen- 
dent implementations in CalcHep and FeynRules. In addition, the FeynRules — > 
WHIZARD driver also developed in the context of this thesis and presented in 
chapter I5T41 will greatly simplify the future implementation of other models of new 
physics into WHIZARD / O'Mega. 

Due to the constraints on the masses and couplings of the new particles, the 
LHC phenomenology of the Three-Site Model turns out to be quite interesting even 
though the new structure cannot compete with other extension of the Standard 
Model (e.g. SUSY) in terms of complexity. Using our WHIZARD / O'Mega imple- 
mentation, we have performed a study of the discovery prospects of the new heavy 
particles in three different types of processes, the results of which we presented in 
chapters |6]-|8l 

For the heavy gauge bosons, essentially two classes of production mechanisms 
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exist with the W'/Z' coupling either to a Standard Model gauge boson or to a Stan- 
dard Model fermion. As a representative of the first class, we have showcased 
a simulation of the W strahlung process in chapter [6l We find that this process 
indeed has the potential of uncovering the W at the LHC. However, while a 5<r 
discovery of the W should be possible in this process with considerably less than 
100 fb" 1 for W masses of between 380 GeV and 500 GeV (the first value being the 
lower bound on mw from existing experimental data), the integrated luminosity 
necessary for a 5a discovery rises to 400 fb -1 for m\yi = 600 GeV (which is the W 
mass where the heavy fermions move above the UV cutoff). Considering that no 
detector effects are taken into account in these numbers, mw' = 600 GeV should 
present the upper limit on the reach of the LHC in this process. 

As the W' strahlung process is essentially mediated by the gauge sector and 
does not contain any information on the fermiophobic couplings of the heavy gauge 
bosons to the Standard Model fermions, we have performed another set of simula- 
tions presented in chapter [7] where we studied the production of the W'/Z' in the 
s channel. For the Z', the only viable discovery channel leads to the final state jjlv 
which our simulations show to be suitable for the discovery of the Z' with the first 
10 — 20 fb^ 1 over the whole parameter space. For the W, three potential discovery 
channels exists with the final state jjlv, jjll and lllv. However, the separation of 
the W resonance from the Z' peak in jjlv turns out to be difficult due to the fi- 
nite experimental resolution of the jet momenta, leaving jjll and lllv as the most 
promising channels. We find that the discovery potential of these ranges from "pos- 
sible within the first 100 fb -1 " to "utterly impossible" for all W masses depending 
on the exact value of the derealization parameter e^. As the previously discussed 
W' strahlung process allows for a discovery of the W' independent of the dereal- 
ization parameter, the combination of these two processes might not only facilitate 
the discovery of the W but even give some hints about its fermiophobic nature. 

Apart from the heavy gauge bosons, the second major component of the new 
structure in the Three-Site Model are the heavy partners to the Standard Model 
fermions. In chapter [8]. a simulation of the t channel induced production of heavy 
quarks leading to either a jjll oder a jjlv final state was discussed. The result 
suggests that, while the mass and width of the heavy fermions makes them hard 
to be detected as a resonance, such a measurement might indeed be possible at 
least for values of Mbuik below 2.9 TeV. However, the success of this in the "real 
world" would be highly dependent on the exact value of Mbuik and the accuracy 
of our understanding of the expected Standard Model background which must be 
subtracted in order to establish the resonance. Still, the situation is not hopeless, 
and if physics reminiscent of the Three-Site Model W'/Z' were to be found at the 
LHC, the heavy fermions could in principle be probed for. 

Combining the results of all our simulations, it is safe to say that the LHC would 
be capable of probing for all major components of the new physics in the Three-Site 
Model. However, there are areas in parameter space which would remain partially 
dark at the LHC, and while the W' and Z' should eventually turn up, the heavy 
fermions might remain completely invisible. 
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Of course, the Three-Site Model is just one example of Higgsless new physics, 
and there is no reason why nature should choose this particular scenario (or why 
nature should be Higgsless anyway). Still, even in the likely case of it not being 
a valid description of nature, the the Three-Site Model remains an interesting ex- 
ample of a model of new physics which is quite successful at (partially) hiding 
its new structure from quick discovery. In fact, as all Higgsless models in which 
heavy vectors are responsible for the delay of unitarity violation must manage to 
pass the precision constraints, phenomena like the fermiophobic couplings of the 
heavy gauge bosons or additional heavy fermions occur in other such models as 
well, with similar difficulties concerning the detection at the LHC. 
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Appendix A 

Notation and Conventions 



All covariant derivatives appearing in this work are defined with a minus sign 

= — igA^ (A.l) 

In several places, the generators of the fundamental representation of SU(2) are 
referred to which are defined as 

n = f (A.2) 

with the Pauli matrices <Tj. Left- and right-handed fermions are defined via the 
chirality projectors U± as 

*l/r = n±* = (A.3) 

In the introduction to 5D gauge theories given in chapter|3l a convention is used 
of Greek Lorentz indices going from to 3 (thus labeling the 4D part of 5D vectors), 
while Roman Lorenz indices go from to 4. The 5th coordinate is labeled either as 
x 5 , y or (implicately in sums) x 4 . For the metric, the mostly minus convention is 
used 

«T = diag (1,-1,-1, -1)^ , 5 ab = diag(l,-l,-l,-l,-ir fe (A.4) 



It is also necessary to establish a convention on the nomenclature regarding the 
Standard Model particles and their heavy partners. If not explicitly stated differ- 
ently in the text, we refer to the Standard Model particles as "light" or "KK light" 
and to their partners as "heavy", "KK heavy" or "KK particles". 
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Appendix B 

A sample spectrum 



The following is a sample spectrum generated by running the spekt rum program 
described in appendix lD.il It is calculated for ideal derealization raw = 500 GeV 
and Mb u ik = 3.5 TeV and gives all masses, wavefunctions, widths and couplings 
(excluding QED and QCD). The heavy quark widths include the one loop QCD 
correction (c.f. chapter l4~4l and appendixICl 

Note that the scale v given in the first few lines must be divided by v2 to obtain 
the symmetry breaking scale v defined in chapter l33l due to a difference in conven- 
tion. For similar reasons, the parameter lambda must be multiplied by 2 to get the 
bulk Yukawa coupling A. 

gO = 0.6915048 gl = 2.0914001 g2 = 0.3565523 

-> x = 0.3306421 t = 0.5156179 e = 0.3133290 

-> eps_L = . 2369703 
v = 235.7870519 lambda = 10.4962241 



electron neutrino : light lepton with isospin +1/2 belonging to generation 
eps_r: 0.0000000 

wavefunction L : -0.9730523 0.2305845 

wavefunction R : 0.0000000 1.0000000 

mass: 0.0000000 
width: 0.0000000 



heavy electron neutrino 
eps_r: 0.0000000 
wavefunction L : 
wavefunction R : 1 
mass : 3596 . 9289707 
width: 628.3857980 



heavy lepton with isospin +1/2 belonging to generation 



2305845 
0000000 



. 9730523 
. 0000000 



up quark : light quark with isospin +1/2 belonging to generation 
eps_r: 0.0000000 

wavefunction L : -0.9730523 0.2305845 

wavefunction R : 0.0000000 1.0000000 

mass: 0.0000000 
width: 0.0000000 



heavy up quark : heavy quark with isospin +1/2 belonging to generation 
eps_r: 0.0000000 

wavefunction L : 0.2305845 0.9730523 

wavefunction R : 1.0000000 0.0000000 

mass : 3596 . 9289707 
width: 585.6279797 



electron : light lepton with isospin -1/2 belonging to generation 
eps_r: 0.0000000 

wavefunction L : -0.9730523 0.2305845 

wavefunction R : 0.0000000 1.0000000 

mass: 0.0000000 
width: 0.0000000 



121 



122 



heavy electron : heavy lepton with isospin -1/2 belonging to generation 
eps_r: 0.0000000 

wavefunction L : 0.2305845 0.9730523 

wavefunction R : 1.0000000 0.0000000 

mass: 3596.9289707 
width: 628.3857980 

down quark : light quark with isospin -1/2 belonging to generation 
eps_r: 0.0000000 

wavefunction L : -0.9730523 0.2305845 

wavefunction R : 0.0000000 1.0000000 

mass: 0.0000000 
width: 0.0000000 

heavy down quark : heavy quark with isospin -1/2 belonging to generation 
eps_r: 0.0000000 

wavefunction L : 0.2305845 0.9730523 

wavefunction R : 1.0000000 0.0000000 

mass: 3596.9289707 
width: 585.6279797 

muon neutrino : light lepton with isospin +1/2 belonging to generation 1 
eps_r: 0.0000000 

wavefunction L : -0.9730523 0.2305845 

wavefunction R : 0.0000000 1.0000000 

mass: 0.0000000 
width: 0.0000000 

heavy muon neutrino : heavy lepton with isospin +1/2 belonging to generation 1 
eps_r: 0.0000000 

wavefunction L : 0.2305845 0.9730523 

wavefunction R : 1.0000000 0.0000000 

mass: 3596.9289707 
width: 628.3859271 

charm quark : light quark with isospin +1/2 belonging to generation 1 
eps_r: 0.0015489 

wavefunction L : -0.9730524 0.2305840 

wavefunction R : -0.0014665 0.9999989 

mass: 1.2500000 
width: 0.0000000 

heavy charm quark : heavy quark with isospin +1/2 belonging to generation 1 
eps_r: 0.0015489 

wavefunction L : 0.2305840 0.9730524 

wavefunction R : 0.9999989 0.0014665 

mass: 3596.9328385 
width: 585.6454626 

muon : light lepton with isospin -1/2 belonging to generation 1 
eps_r: 0.0001313 

wavefunction L : -0.9730523 0.2305845 

wavefunction R : -0.0001244 1.0000000 

mass: 0.1060000 
width: 0.0000000 

heavy muon : heavy lepton with isospin -1/2 belonging to generation 1 
eps_r: 0.0001313 

wavefunction L : 0.2305845 0.9730523 

wavefunction R : 1.0000000 0.0001244 

mass: 3596.9289985 
width: 628.3858577 

strange quark : light quark with isospin -1/2 belonging to generation 1 
eps_r: 0.0011771 

wavefunction L : -0.9730524 0.2305842 

wavefunction R : -0.0011145 0.9999994 

mass: 0.9500000 
width: 0.0000000 

heavy strange quark : heavy quark with isospin -1/2 belonging to generation 1 
eps_r: 0.0011771 

wavefunction L : 0.2305842 0.9730524 

wavefunction R : 0.9999994 0.0011145 

mass: 3596.9312047 
width: 585.6492771 

tauon neutrino : light lepton with isospin +1/2 belonging to generation 2 
eps_r: 0.0000000 

wavefunction L : -0.9730523 0.2305845 

wavefunction R : 0.0000000 1.0000000 

mass: 0.0000000 
width: 0.0000000 
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heavy tauon neutrino : heavy lepton with isospin +1/2 belonging to generation 2 
eps_r: 0.0000000 

wavefunction L : 0.2305845 0.9730523 

wavefunction R : 1.0000000 0.0000000 

mass: 3596.9289707 
width: 628.4222032 

top quark : light quark with isospin +1/2 belonging to generation 2 
eps_r: 0.2202501 

wavefunction L : -0.9752869 0.2209422 

wavefunction R : -0.2046066 0.9788443 

mass: 174.0000000 
width: 1.5230000 

heavy top quark : heavy quark with isospin +1/2 belonging to generation 2 
eps_r: 0.2202501 

wavefunction L : 0.2209422 0.9752869 

wavefunction R : 0.9788443 0.2046066 

mass: 3674.4892058 
width: 739.9590965 



tauon : light lepton with isospin 
eps_r: 0.0022056 
wavefunction L : -0.9730525 
wavefunction R : -0.0020883 
mass: 1.7800000 
width: 0.0000000 



■1/2 belonging to generation 2 



0.2305835 
0.9999978 



heavy tauon : heavy lepton with isospin 
eps_r: 0.0022056 
wavefunction L : 0.2305835 
wavefunction R : 0.9999978 
mass: 3596.9368138 
width: 628.4026301 



-1/2 belonging to generation 2 



. 9730525 
0.0020883 



bottom quark : light quark with isospin -1/2 belonging to generation 2 
eps_r: 0.0052042 

wavefunction L : -0.9730536 0.2305789 

wavefunction R : -0.0049275 0.9999879 

mass: 4.2000000 
width: 0.0000000 

heavy bottom quark : heavy quark with isospin -1/2 belonging to generation 2 
eps_r: 0.0052042 

wavefunction L : 0.2305789 0.9730536 

wavefunction R : 0.9999879 0.0049275 

mass: 3596.9726381 
width: 910.9435197 



W boson : light W boson 
wavefuntion: 0.9858825 
mass: 80.4030000 
width: 2.0480000 



0.1674384 



heavy w boson : heavy w boson 
wavefuntion: -0.1674384 
mass: 500.0000000 
width: 5.0007454 



0.9858825 



Z boson : light Z boson 
wavefuntion: -0.8759133 
mass: 91.1880000 
width: 2.4430000 



heavy z boson : heavy z boson 
wavefuntion: 0.1657276 
mass: 501.6804843 
width: 4.7681804 



-0 . 9827488 



0.0820918 



w boson: light 
isospin +1/2 lepton: light 
isospin -1/2 lepton: light 

left-handed coupling matrix (rows <-> iso-up, columns <-> iso-down) : 



0.6641137 
0.0000000 
0.0000000 



0.0000000 
0.6641137 
0.0000000 



0.0000000 
0.0000000 
0.6641137 



right-handed coupling matrix (rows <-> iso-up, columns <-> iso-down) : 



0.0000000 
0.0000000 



0.0000000 
0.0000000 



0.0000000 
0.0000000 
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w boson: 

isospin + 1/2 quark: 
isospin -1/2 quark: 



light 
light 
light 



left-handed coupling matrix (rows <-> iso-up, columns <-> iso-down) : 



.6641137 
0.0000000 
0.0000000 



0.0000000 
0.6641137 
0.0000000 



0.0000000 
0.0000000 

0.6648179 



right-handed coupling matrix (rows <-> iso-up, columns <-> iso-down) : 



0.0000000 
0.0000000 
0.0000000 



0.0000000 
0.0000006 
0.0000000 



0.0000000 
0.0000000 
0.0003530 



w boson: light 
isospin +1/2 lepton: light 
isospin -1/2 lepton: heavy 

left-handed coupling matrix (rows <-> iso-up, columns <-> iso-down) : 



-0.0743928 
0.0000000 
0.0000000 



0.0000000 
-0.0743928 
0.0000000 



0.0000000 
0.0000000 
-0.0743921 



right-handed coupling matrix (rows <-> iso-up, columns <-> iso-down) : 



0.0000000 
0.0000000 
0.0000000 



0.0000000 
0.0000000 
0.0000000 



0.0000000 
0.0000000 
0.0000000 



w boson: 

isospin +1/2 quark: 
isospin -1/2 quark: 



light 
light 
heavy 



left-handed coupling matrix (rows <-> iso-up, columns <-> iso-down) : 



-0 . 0743928 
0.0000000 
0.0000000 



0.0000000 
-0 . 0743928 
0.0000000 



0.0000000 
0.0000000 
-0.0780258 



right-handed coupling matrix (rows <-> iso-up, columns <-> iso-down) : 



0.0000000 
0.0000000 
0.0000000 



0.0000000 
-0.0005135 
0.0000000 



0.0000000 
0.0000000 
-0 . 0716484 



W boson: light 
isospin +1/2 lepton: heavy 
isospin -1/2 lepton: light 

left-handed coupling matrix (rows <-> iso-up, columns <-> iso-down) : 



-0 . 0743928 
0.0000000 
0.0000000 



0.0000000 
-0 . 0743928 
0.0000000 



0.0000000 
0.0000000 
-0.0743932 



right-handed coupling matrix (rows <-> iso-up, columns <-> iso-down) : 



0.0000000 
0.0000000 
0.0000000 



0.0000000 
-0.0000435 
0.0000000 



0.0000000 
0.0000000 
-0.0007313 
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w boson: 
isospin 



1/2 quark: 
1/2 quark: 



light 
heavy 
light 



left-handed coupling matrix (rows <-> iso-up, columns <-> iso-down) : 



-0.0743928 
0.0000000 
0.0000000 



0.0000000 
-0.0743926 
0.0000000 



0.0000000 
0.0000000 
-0.0678180 



right-handed coupling matrix (rows <-> iso-up, columns <-> iso-down) : 



0.0000000 
0.0000000 
0.0000000 



0.0000000 
-0.0003903 
0.0000000 



0.0000000 
0.0000000 
-0.0016890 



W boson: light 
isospin +1/2 lepton: heavy 
isospin -1/2 lepton: heavy 

left-handed coupling matrix (rows <-> iso-up, columns <-> iso-down) : 



0.3678096 
0.0000000 
0.0000000 



0.0000000 
0.3678096 
0.0000000 



0.0000000 
0.0000000 

0.3678095 



right-handed coupling matrix (rows <-> iso-up, columns <-> iso-down) : 



. 3501807 
0.0000000 
0.0000000 



o . ooooooo 

0.3501807 
0.0000000 



o . ooooooo 

0.0000000 
0.3501799 



w boson: 

isospin +1/2 quark: 
isospin -1/2 quark: 



light 
heavy 
heavy 



left-handed coupling matrix (rows <-> iso-up, columns <-> iso-down) : 



0.3678096 
0.0000000 
0.0000000 



0.0000000 
0.3678095 
0.0000000 



0.0000000 
0.0000000 
0.3670548 



right-handed coupling matrix (rows <-> iso-up, columns <-> iso-down) : 



. 3501807 

. ooooooo 

0.0000000 



0.0000000 
0.3501801 
0.0000000 



0.0000000 
0.0000000 

0.3427682 



w boson: heavy 
isospin +1/2 lepton: light 
isospin -1/2 lepton: light 

left-handed coupling matrix (rows <-> iso-up, columns <-> iso-down) : 



. ooooooo 
0.0000000 
0.0000000 



o . ooooooo 

0.0000000 
0.0000000 



o . ooooooo 
0.0000000 
-0.0000005 



right-handed coupling matrix (rows <-> iso-up, columns <-> iso-down) : 



0.0000000 
0.0000000 
0.0000000 



0.0000000 
0.0000000 
0.0000000 



0.0000000 
0.0000000 
0.0000000 



w boson: 

isospin +1/2 quark: 



heavy 
light 
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isospin -1/2 quark: light 

left-handed coupling matrix (rows <-> iso-up, 



0.0000000 
0.0000000 
0.0000000 



0.0000000 
-0.0000004 
0.0000000 



columns <-> iso-down) : 

0.0000000 
0.0000000 
-0.0048388 



right-handed coupling matrix (rows <-> iso-up, columns <-> iso-down) : 



0.0000000 
0.0000000 
0.0000000 



0.0000000 
0.0000034 
0.0000000 



0.0000000 
0.0000000 
0.0020788 



w boson: heavy 
isospin +1/2 lepton: light 
isospin -1/2 lepton: heavy 

left-handed coupling matrix (rows <-> iso-up, columns <-> iso-down) : 



.4886032 
0.0000000 
0.0000000 



0.0000000 
0.4886032 
0.0000000 



. 0000000 
0.0000000 

0.4886032 



right-handed coupling matrix (rows <-> iso-up, columns <-> iso-down) : 



0.0000000 
0.0000000 
0.0000000 



0.0000000 
0.0000000 
0.0000000 



0.0000000 
0.0000000 
0.0000000 



w boson: 

isospin +1/2 quark: 
isospin -1/2 quark: 



heavy 
light 
heavy 



left-handed coupling matrix (rows <-> iso-up, columns <-> iso-down) : 



0.4886032 
0.0000000 
0.0000000 



0.0000000 
0.4886022 
0.0000000 



0.0000000 
0.0000000 
0.4693173 



right-handed coupling matrix (rows <-> iso-up, columns <-> iso-down) : 



0.0000000 
0.0000000 
0.0000000 



0.0000000 
-0.0030238 
0.0000000 



0.0000000 
0.0000000 
-0.4218681 



w boson: heavy 
isospin +1/2 lepton: heavy 
isospin -1/2 lepton: light 

left-handed coupling matrix (rows <-> iso-up, columns <-> iso-down) : 



0.4886032 
. 0000000 
0.0000000 



0.0000000 
.4886032 
0.0000000 



. 0000000 
. 0000000 
.4886012 



right-handed coupling matrix (rows <-> iso-up, columns <-> iso-down) : 



0.0000000 
0.0000000 
0.0000000 



0.0000000 
-0.0002564 
0.0000000 



0.0000000 
0.0000000 
-0.0043058 



w boson: heavy 
isospin +1/2 quark: heavy 
isospin -1/2 quark: light 



left-handed coupling matrix (rows <-> iso-up, columns <-> iso-down) : 
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0.4886032 
0.0000000 
0.0000000 

right-handed coupling 

. 0000000 
0.0000000 
0.0000000 



0.0000000 
0.4886026 
0.0000000 

matrix (rows <-> iso-up, 

. 0000000 
-0 . 0022981 
0.0000000 



0.0000000 
0.0000000 
.4885680 

columns <-> iso-down) : 



0.0000000 
-0.0099449 



w boson: heavy 
isospin +1/2 lepton: heavy 
isospin -1/2 lepton: heavy 

left-handed coupling matrix (rows <-> iso-up, columns 



1 . 9460904 
0.0000000 
0.0000000 



0.0000000 
1 . 9460904 
0.0000000 



0.0000000 
0.0000000 

1.9460909 



right-handed coupling matrix (rows <-> iso-up, columns <-> iso-down) : 



2.0618749 
0.0000000 
0.0000000 



0.0000000 
2.0618749 
0.0000000 



0.0000000 
0.0000000 
2 .0618704 



w boson: 

isospin +1/2 quark: 
isospin 



quark: 



heavy 
heavy 
heavy 



left-handed coupling matrix (rows <-> iso-up, columns <-> iso-down) : 



1.9460904 
0.0000000 
0.0000000 



0.0000000 
1.9460908 
0.0000000 



0.0000000 
0.0000000 
1.9508340 



right-handed coupling matrix (rows <-> iso-up, columns <-> iso-down) : 



2 .0618749 
0.0000000 
0.0000000 



0.0000000 
2 .0618714 
0.0000000 



0.0000000 
. 0000000 
2.0182299 



couplings to the light Z 



2x electron neutrino: 
L: -0.3765880 



0.0000000 



2x up quark: 
L: -0.2648388 



0.1117492 



2x heavy electron neutrino: 
L: -0.2072812 R: 



2x heavy up quark: 
L: -0.0955319 



-0.0854589 



electron neutrino , heavy electron 



neutrino : 
0000000 



up quark , heavy up quark: 
L: 0.0425077 R: 



0.0000000 



2x electron: 
L: 0.2089642 



-0.1676238 



2x down quark: 
L: 0.3207134 



2x heavy electron: 
L: 0.0396573 



0.0295843 



2x heavy down quark: 
L: 0.1514065 
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electron , heavy electron: 
L: -0.0425077 R: 



down quark , heavy down quark: 
L: -0.0425077 R: 



2x muon neutrino: 
L: -0.3765880 



0.0000000 



2x charm quark: 
L: -0.2648388 



2x heavy muon neutrino: 
L: -0.2072812 



-0 .1972081 



2x heavy charm quark: 
L: -0.0955319 



muon neutrino , heavy muon neutrino: 
L: 0.0425077 R: 0.0000000 



charm quark , heavy charm quark: 
L: 0.0425076 R: 



0.0002892 



2x strange quark: 
L: 0.3207134 



-0.0558744 



2x heavy muon: 
L: 0.0396573 



2x heavy strange quark: 
L: 0.1514065 



muon , heavy muon: 
L: -0.0425077 



-0.0000245 



strange quark , heavy strange quark: 
L: -0.0425077 R: -0.0 



2x tauon neutrino: 
L: -0.3765880 



0.0000000 



2x top quark: 
L: -0.2656636 



2x heavy tauon neutrino: 
L: -0.2072812 



-0 .1972081 



2x heavy top quark: 
L: -0.0947071 



-0 . 0772030 



tauon neutrino , heavy tauon neutrino: 



. 0425077 



0.0000000 



top quark , heavy top quark: 
L: 0.0408237 R: 



0.0394964 



2x bottom quark: 
L: 0.3207138 



-0.0558698 



2x heavy tauon: 
L: 0.0396573 



2x heavy bottom quark: 
L: 0.1514061 



. 1413287 



tauon , heavy tauon: 
L: -0.0425075 



-0.0004118 



bottom quark , heavy bottom quark: 
L: -0.0425067 R: -0.0009717 



couplings to the heavy Z 



2x electron neutrino: 
L: -0.0150208 R: 



0.0000000 
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2x up quark: 
L: 0.0044925 



0.0195133 



2x heavy electron neutrino: 
L: -0.9846090 R: 



-1 . 0422955 



2x heavy up quark: 
L: -0.9650956 



-1 . 0227822 



electron neutrino , heavy electron neutrino: 



-0.2434336 



0.0000000 



up quark , heavy up quark: 
L: -0.2434336 R: 



2x electron: 
L: -0.0142492 



-0.0292700 



2x down quark: 
L: 0.0052642 



2x heavy electron: 
L: 0.9553390 



2x heavy down quark: 
L: 0.9748523 



1.0325389 



electron , heavy electron: 
L: 0.2434336 R: 



down quark , heavy down quark: 
L: 0.2434336 R: 



0.0000000 



2x muon neutrino: 
L: -0.0150208 



2x charm quark: 
L: 0.0044928 



2x heavy muon neutrino: 
L: -0.9846090 



-1 . 0422955 



2x heavy charm quark: 
L: -0.9650959 



muon neutrino , heavy muon neutrino: 
L: -0.2434336 R: 0.0000000 



charm quark , heavy charm quark: 



-0.0142492 



-0.0292700 



2x strange quark: 
L: 0.0052640 



-0.0097554 



2x heavy muon: 
L: 0.9553390 



2x heavy strange quark: 
L: 0.9748525 



1.0325376 



muon , heavy muon : 
L: 0.2434336 



strange quark , heavy strange quark: 



0.2434333 



-0.0011617 



2x tauon neutrino: 
L: -0.0150208 



2x top quark: 
L: 0.0092162 



2x heavy tauon neutrino: 
L: -0.9846090 



-1.0422955 



2x heavy top quark: 
L: -0.9698193 



tauon neutrino , heavy tauon neutrino: 
L: -0.2434336 R: 0.0000000 
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top quark , heavy top quark: 
L : -0.2337897 R: 0.2087489 



2x tauon: 

L: -0.0142497 R: -0. 02 92 655 



2x bottom quark: 
L: 0.0052614 



2x heavy tauon : 
L: 0.9553395 R: 1. 0130210 



2x heavy bottom quark: 
L: 0.9748551 



tauon , heavy tauon: 
L: 0.2434326 R: -0.0021766 



bottom quark , heavy bottom quark: 
L: 0.2434281 R: -0.0051358 



couplings between gauge bosons 

2xW boson , Z boson : -0.5950754 

2xW boson , heavy Z boson : 0.0537663 

2xheavy W boson , Z boson : -0.2374152 

2xheavy W boson , heavy Z boson : -1.9944859 

W boson , heavy W boson , Z boson : 0.0625477 

W boson , heavy W boson , heavy Z boson : -0.3581991 

4x W boson , Ox heavy W boson : 0.4551806 

3x W boson , lx heavy W boson : -0.0564797 

2x W boson , 2x heavy W boson : 0.1322188 

lx W boson , 3x heavy W boson : 0.6995732 

Ox W boson , 4x heavy W boson : 4.1325152 

2x W boson , Z boson , photon : -0.1864544 

2x W boson , heavy Z boson , photon : 0.0168466 

2x heavy W boson , Z boson , photon : -0.0743891 

2x heavy W boson , heavy Z boson , photon : -0.6249303 

W boson , heavy W boson , Z boson , photon : 0.0195980 

W boson , heavy W boson , heavy Z boson , photon : -0.1122342 

2x W boson , 2x Z boson : 0.3580269 

2x W boson , 2x heavy Z boson : 0.1311974 

2x W boson , Z boson , heavy Z boson : -0.0543996 

2x heavy W boson , 2x Z boson : 0.0602782 

2x heavy W boson , 2x heavy Z boson : 4.1062808 

2x heavy W boson , Z boson , heavy Z boson : 0.4511168 

W boson , heavy W boson , 2x Z boson : -0.0520704 

W boson , heavy W boson , 2x heavy Z boson : 0.6951640 

W boson , heavy W boson , Z boson , heavy Z boson : 0.0884049 



Appendix C 

Two-body decays — Analytical 
Results 



In this chapter the analytical results for the decay widths presented in chapter 14.41 
are given. 

C.l Tree level 

Consider a heavy particle with momentum p\ and mass mi decaying into two par- 
ticles with momenta P2 I p% and masses 771,2 / tti^. 

P2 
P3 

In the center-of-mass system, we can parameterize the momenta as 

# = (mi, 0,0,0)" , p% = (E 2 ,p)» , p» = (E 3 ,-p) 

Furthermore, momentum conservation and the mass shell conditions can be lever- 
aged to express the energies and the modulus of the three-momentum as 

m? + mn — mi „ m? + m 2 . — mk 



-r ,,l 2 - ,,o 3 ,,u x 
— , 



2mi " 2m\ 

\p\ 2 = -^-2 ( m i ~ ( m 2 + m 3 ) 2 ^ (m 2 - (m 2 - m 3 ) 2 ) (C.l) 

The total width of the decaying particle can be expressed as the product of a 
phasespace factor and the square of the transition matrix element \M.\ 2 (see e.g. 

El) 



167T ml 




l 

dcosO \M\ 2 
1 
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C.l TREE LEVEL 



with the product Y[f^ m f running over all external fermions. After performing 
the spin sum (including an averaging factor -h in front of T), Lorentz invariance 
requires \M. | 2 to be independent of 8. Inserting (|C.1|> and adopting the convention 
of absorbing the fermion factor into the squared matrix element, we arrive at the 
result 

1 1 



L 



N 16ir m 



^3 (m\ - (m 2 + m 3 ) 2 ) {m\ - (m 2 - m 3 )' 



\M\' 



which is completely expressed in terms of the particle masses, the squared matrix 
element | M. \ 2 and the averaging factor i which is determined from the spin of the 
initial state particl^. 

The remaining quantity that has to be determined is the squared matrix element 
\M\ 2 which is given in the following for the different spin combinations that are 
relevant for calculating the widths of the KK particles in the Three-Site Model. 

Spin 1 — > Spin 1 + Spin 1 

• Coupling: standard triple gauge boson coupling 

• Mass assignment: 



m 2 



Amplitude: 




m 3 



spin 



mi 



ran 



A i A A i A A i A 

771-2 + m 3 m l + m 3 m l + m 2 



TJln 



i 2 m^ in-; in 

o I 



2 2 2 2 



in-, 



rrin 



+ 



m 



+ 



Uln, 



rrin 



2 2 



Spin 1 — > Spin \ + Spin \ 

• Coupling: g v - g A j 5 

• Mass assignment: 



mi 



m 2 




m 3 



1 An additional symmetry factor of \ could arise if there were identical particles in the final state. 
However, this does not occur for the decays relevant to the calculation presented here. 
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Amplitude: 

\M\ 2 = (g 2 v + g\) ( 2 (2m? - m\ - m 2 ) - 2^±^ + 4^) 



spins 



+ 12 ig\ - g 2 A) mim 3 



Spin \ — > Spin ^ + Spin 1 

• Coupling: g v - g A j 5 

• Mass assignment: 



m 2 




mi 

m 3 

Amplitude: 

\M\ 2 = (g 2 v + g\) (2 (mj + m 2 3 - 2m 2 2 ) + 2 

spins 



4 1 4 2 2 

mi + rrij m\m^ 
2 2 

771.2 m 2 



12(^- 3 i)mim 3 (C.2) 



C.2 (9 (a s ) corrections to the heavy quark widths 

According to (|4.31|) , the calculation of the O (a s ) QCD corrections to the heavy 
quark widths can be split into the calculation of the virtual correction and that of 
the real ones. In order to subtract the divergences arising in the loop integrals, we 
have to specify a renormalization scheme. For this calculation, we use a modified 
MS scheme where the heavy particle masses are renormalized on-shell such that 
the masses obtained by diagonalizing the Lagrangian are treated as pole masses. 

The consistency of the result has been checked by confirming infrared finite- 
ness. In addition, it has been cross-checked against the calculation of the O (a s ) 
correction to the top decay in 1651 . 



Virtual corrections 

For computing the corrections which arise from the exchange of virtual gluons, we 
decompose the squared and spin summed two body decay matrix element as 

Y / \M\ 2 = M + a s (Mi,ipi + M MZ ) + O {a 2 ) (C.3) 

spin 
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level matrix element 

2 

given by (IC.2D . The second term is the interference between the tree level matrix 
element and the QCD vertex correction 




The third term M± ; sz arises from the wave function renormalization 




In order to better understand the origin of this term and how to calculate it, con- 
sider the full two point function E(p 2 ) of a scalaJl 



The first term in the decomposition is just the tree 




In the vicinity of the mass pole, S has a Laurent expansion 

Zi 



np 1 



p2 _ m 2 



+ 0(1 



with the pole mass m and residual Z. At tree level, the residual is Z = 1. However, 
loop corrections give rise to a nontrivial value of Z which depends on the renor- 
malization scheme chosen. If we apply LSZ reduction (see e.g. HI) and amputate 
legs from a Green's function to obtain a transition matrix element, we have to keep 
track of the potentially nontrivial residual 



M 



—= lim (p 2 

V Z p 2 ->m 2 





(C.4) 



(|C.4|) shows that we have to multiply the amputated Greens function with the 
proper residual \[Z for each external leg in order to correctly obtain the transition 



2 Choosing a scalar is convenient to make the argument more transparent. Of course, the same 
reasoning applies to fields of arbitrary spin. 



CHAPTER C. TWO-BODY DECAYS — ANALYTICAL RESULTS 



135 



matrix element. This is precisely the origin of Mxsz- H we work in perturbation 
theory, then the residual can be expanded as 

Z = l + a s AZ x + O (a 2 s ) (C.5) 

and M\fiz is correctly calculated as 

M 1>sz = M {AZ ljf , + AZ XJ ) 

with the residuals of the two fermions / and /'. 

For the actual calculation of the residual, decompose the one particle irreducible 
(1PI) two point function as 

iU(p 2 ) = » ^ » = i(ftA{p 2 )+mZ(p 2 )) (C.6) 

with scalar functions A and S. The full two point function is then obtained by 
Dyson resummation of the 1PI function 

i , . i 1 

+ z ( iU ) z + • - • 



ft — m ft — m ft — m ft—ml + j^^ ft — m + U 

Using the decomposition (|C.6|) , this can be more properly written as 



1 + A 2 2/"l-S N2 
P 171 ( T+A 



As we have chosen to renormalize the mass on- shel§ the self energy S must have 
a pole at p 2 = m 2 , and we can perform Laurent expansions around m 

1 ~ ^ p2 ^ =i + o(v 2 - m 2 ) 
1 + A(p2) + U ^ m ^ 



1 i iZ 



S(p2) 1 + A p2 _ m2 ^V P 2 -m 2 



+ O (1) (C.7) 



(|C.7|) implies that the self energy can be expanded around the pole as 

a iZ 



m 



+ 0(1) 



3 If the calculation is performed in a different renormalization scheme, then we can always per- 
form a finite renormalization of the mass 

m — > m — mSm , E — > S + 8m + O (ce s ) 

such that the on-shell condition is fulfilled. As the final O (a s ) result dC8t only depends on the 
derivative £' which is independent on the finite renormalization 5m to order O (a s ), it is in fact valid 
for any renormalization scheme (provided that we set m to the pole mass in <C.8t ). 
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with Z being the residual we are looking for and which therefore can be calculated 

as 

J = T* B W 2 2 

Z dp* p 2 =m 2 

Plugging in l|C.7|l and expanding in a s , we finally obtain 

AZ = -A(m 2 ) - 2m 2 (A'(m 2 ) + S'(m 2 )) (C.8) 



The actual calculation of M\,xyi ar *d was carried out using FeynArts J57 
and FormCalc II66II . The result for M\ ipi is 



M 



i.ipi — 
+ 



37rm 2 



((mi - m 3 ) 2 - mfj 



+ B (mi 2 + 2m 3 mi + 2m 2 2 + m 3 2 ) 



+ Cq (mi 4 + 2m3?n,i 3 + (m 2 2 + 2m3 2 ) mi 2 + (2m^ 3 — 2m2 2 m 3 ) m\ 



- 2m 2 4 + m 3 4 + m 2 2 m 3 2 ) 



+ Ci (mi + 2m 3 m\ + (m 2 + 2m3 ) mi + m3 (m 2 + 27773 J mi 
— 2m 2 4 + m.3 4 + 4m 2 2 m 3 2 ) 

(A Q / O 0\ O / O 0\ 

mi + 2m3mi + 2 (2m 2 + 1113 ) mi + m.3 (m 2 + 2m3 J mi 



-2m 2 4 + m 3 4 + m 2 2 m 3 2 ) 



- Coo (2 (mi 2 + 2m 3 mi + 2m 2 2 + m 3 2 )) 

- Cn (m 3 (mi + m.3) (mi 2 + 2m 3 mi - m 2 2 + m 3 2 )) 

- C12 ((mi + m 3 ) 2 (mi 2 + 2m 3 m\ - m 2 2 + m 3 2 )) 

- C 22 (mi (mi + 771,3) (mi 2 + 2m 3 m\ - m 2 2 + m 3 2 )) 

+ i^((mi + m 3 ) 2 -m 2 
String \ 

+ So (mi 2 - 2m 3 mi + 2m 2 2 + m 3 2 ) 

+ Co (mi 4 — 2m,3777 i 3 + (m 2 2 + 2m3 2 ) mi 2 + 27773 (m 2 2 — m3 2 ) mi 

-2m 2 4 + m 3 4 + m 2 2 777 3 2 ) 
+ Ci (mi — 277737771 + (m 2 + 2m3 ) mi — 7773 (7772 + 27773 ) mi 

— 2m 2 4 + WI3 4 + 4m 2 2 m3 2 ) 
+ C 2 (mi 4 — 2m3mi 3 + 2 (2m 2 2 + m3 2 ) mi 2 — 777.3 (m 2 2 + 2m3 2 ) mi 

-2m 2 4 + m 3 4 + m 2 2 m 3 2 ) 
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Coo (2 (mi 2 - 2m 3 mi + 2m 2 2 + m 3 2 )) 

Cn ((771.3 - mi)m 3 (mi 2 - 2m 3 mi - m 2 2 + m 3 2 )) 

C 12 ((mi - m 3 ) 2 (mi 2 - 2m 3 mi - m 2 2 + m 3 2 )) 

C22 (mi(mi - m 3 ) (mi 2 - 2m 3 mi - m 2 2 + m 3 2 )) 



(C.9) 



where the same assignments as in (|C.2|) are used for the masses of the particles. 
Bi, d and appearing in dC9l) are the scalar and tensor integrals defined in II67B 
with 

Bi = Bi (m 2 2 ,ml,m\) 
Ci = Ci (ml, m 2 ,, mf, 0, m§, m 2 ) 
= Cy (m|, m\,m\, 0, m§, mf) 

For the O (a s ) correction to the the residual dC5l) , we obtain the result 



AZ=-^ + ^\ B (m 2 ,0,m 2 )+B 1 (m 2 ,0,m 2 ) - 

2m 2 ^ (B (p 2 ,0,m 2 ) - B 1 (p 2 ,0,m 2 )) 



(C.IO) 



with the fermion mass m. 

Both (|C.9|) and l|C.10|) contain a infrared divergence which originates from the 
vanishing gluon mass and which is understood to be regularized by a small gluon 
mass to. The divergence in uj drops out once the contributions of gluon radiation 
are added, allowing for a finite limit ui — > 0. By the virtue of gauge invariance, 
the dependence on the renormalization scale \i drops out of the NLO result dC.3l) , 
leaving the running of a s as the only source of the scale uncertainty. 



Real corrections 

The remaining piece to calculate is the decay width of the heavy quark going into 
a quark, a heavy gauge boson and an on-shell gluon 



Treal OC 




(IC21 



For the calculation, we choose a mass and momentum assignment analogous to 



Pi, mi 




P2,m 2 



138 



C.2 O {a s ) CORRECTIONS TO THE HEAVY QUARK WIDTHS 



and label the phasespace integrals appearing in the calculation in accordance with 
1571 as 

pi-im _ 1 f d 3 P2 d 3 P3 d 3 q , (±2q Pjl ) . . . (±2qp jm ) 

h ^ ~~ 2 i MM^ HP1+P2+P3 ~ q) (±2^)... (±2^) (C11) 

with the minus signs belonging to p\ and the plus signs to P2/3- Analytical expres- 
sions for the integrals (|C.1H) can be found in |67| . 

The calculation has been carried out with FORM [68| and we choose to decom- 
pose the results as 

(y. 

Treal = ~ Q o 1 {9V Ay + 9aA A ) (C.12) 

Ay is defined in terms of the integrals dC.llb as 

Ay = (m 3 - mi) 2 1+ 

({m?, — mi) 2 — w|J {{jnz + nil) 2 + 2n\^\ x 

(h + h + m\l n + (mj -m 2 2 + mf) I X3 + m|/ 33 ) + 

mi + ^ (™3 - «i) 2 j (4 + if) 

and can be obtained from Ay by swapping the sign of mi (just in the polyno- 
mials, not in the integrals dC.llll ). 



Appendix D 

Code 



This chapter is intended to give an overview over some of the code written in the 
course of this thesis. However, it does not intended to cover the whole source 
which can be downloaded from 

http: //theorie. physik .uni-wuerzburg. de/ "cnspeckn/ 



D.l Coupling Library 

The coupling library is split into four core modules and a number of support mod- 
ules and programs. The core consists of: 

• module tdefs: Global definitions. 

• module threeshl: Actual implementation of the model; calculates masses, 
couplings and widths. 

• module nlowidth: Contains the numerical code for the calculation of the 
O (a s ) QCD corrections, c.f. app. IC.2l 

• module t g 1 u e : Definitions and functions for access from the O'Mega-generated 
amplitudes. 

The support components are 

• module t script: Allows to access to masses, widths, couplings etc. by text 
identifiers like e.g. lhcoupling (T b w) for the left-handed t'bW coupling. 

• program spektrum: Dumps spectrum and couplings for a specific point in 
parameter space (c.f. app.|B|. 

• program threeshleval: Simple command-line tool which evaluates a quan- 
tity as specified by a t script function. 

• program paraplot: Generates a density plot of a tscript function over pa- 
rameter space, taking a list of constraints into account. 
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• scan: Bash script for creating a "catalogue" of density plots. 

• plotter . pi: Perl frontend to threeshleval for plotting quantities over mw' 
or M bulk . 

• threeshl_mlink: Mathematica wrapper around the libraries. Allows to ac- 
cess quantities via t script functions. 

In this section, some pieces of code will be presented to illustrate organization, 
implementation and usage of the package. The interested reader should consult the 
woven documentation generated from the noweb source which can be downloaded 
from the URL quoted at the beginning of this chapter which is more detailed and 
also gives the analytic formulae employed by the module. 

Definitions: module tdef s 

The module tdef s defines some basic constants 

integer, parameter : : double=selected_real_kind & 

(precision ( 1 .) +1 , range (l.)+l) 
integer, parameter : : slength = 256 
real (kind=double) , parameter : : & 

pi=3. 1415 92 65358 97 932385_double 

slength defines a standard length for character variables, and the other two con- 
stants are self-explanatory. In addition, the module defines the output_unit and 
error_onit either via the iso_f ortran_env (a FORTRAN 2003 feature) or as an 
ordinary constant in case this module is not supported by the compiler. 

The actual calculations: module threeshl 

The first duty of this module is to provide the necessary infrastructure to store the 
masses, widths, wavefunctions and couplings. To this end, each particle is assigned 
a number between and 64 whose bits encode different quantum numbers 



bits 


meaning 





0: Light state, 1: Heavy state 


1 


0: Fermion, 1: Boson 


2 


0: Lepton, 1: Quark (0: W, 1: Z) 


3 


0: Isospin \, 1: Isospin — \ 


4-5 


Fermion generation 



integer, parameter, public : : e_bcd=B'001 000', nue_bcd=B'000000', & 
mu_bcd=B'011000', numu_bcd=B'01 0000', tau_bcd=B'1 01 000', & 
nutau_bcd=B'1 00000', he_bcd=B'001 001 ', hnue_bcd=B'000001 ', & 
hmu_bcd=B'01 1001 ', hnumu_bcd=B'01 0001 ', htau_bcd=B'1 01 001 ', & 
hnutau_bcd=B'1 00001', u_bcd=B'0001 00', d_bcd=B'001 1 00', & 
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c_bcd=B'010100', s_bcd=B'011100', t_bcd=B'1 001 00', & 
b_bcd=B'101100', hu_bcd=B'000101', hd_bcd=B'001 1 01 ', & 
hc_bcd=B'010101', hs_bcd=B'011101', ht_bcd=B'1 001 01 ', & 
hb_bcd=B'101101', w_bcd=B'010', hw_bcd=B'01 1 z_bcd=B'1 1 0', & 
hz_bcd=B'1 1 1 ', a_bcd=63 

(although it doesn't fit into the quantum number scheme, the photon snaps nicely 
into its place). For the indexing of the coupling arrays, index ranges are defined 
which are nontrivial and disjunct in order to allow for better compile-time and 
run-time consistency checks 



specifier 


meaning 


l_chir, r_chir 


chirality 


l_mode, h_mode, lh_mode 


KK modes / mode combinations 


iso_up, iso_down 


isospin 


lat_0, lat_l, lat_2 


lattice sites 


gen_0, gen_l, gen_3 


fermion generations 


ftype_l, ftype_q 


lepton / quark 


ptype_b, ptype_f 


boson / fermion 


btype_w, btype_z, btype_a 


W 1 Z 



integer, parameter, public : : l_chir=100, r_chir=101, & 
l_mode=110, h_mode=lll, lh_mode=112 , iso_up=120, & 
iso_down=121, lat_0=130, lat_l=131, lat_2=132, & 
gen_0=140, gen_l=141, gen_2=142, ftype_l=150, & 
ftype_q=151, ptype_b=160, ptype_f=161, btype_w=170, & 
btype_z=171, btype_a=172 

The masses and widths of the particles are then stored in arrays indexed by the 
integer assigned to the particle 

real (kind=double) , public, target : : & 

mass_array ( : 63 ) =0 . , width_array ( : 63 ) =0 . 

while the couplings are stored in arrays which are indexed by the quantum num- 
bers of the lines meeting at the vertex 

real (kind=double) , public, target : : & 

g_w_lep ( l_mode : h_mode, l_mode : h_mode, gen_0:gen_2, & 
l_mode : h_mode, gen_0:gen_2, l_chir : r_chir ) = 0., & 

g_w_quark (l_mode :h_mode, l_mode : h_mode, gen_0:gen_2, & 

l_mode : h_mode, gen_0:gen_2, l_chir : r_chir ) = 0. 
real (kind=double) , public, target : : & 

g_z_lep (l_mode :h_mode, l_mode : lh_mode, gen_0:gen_2, & 
iso_up : iso_down, l_chir : r_chir ) = 0.,& 

g_z_quark ( l_mode : h_mode, l_mode : lh_mode, gen_0:gen_2, & 
iso_up : iso_down, l_chir : r_chir ) = 0. 
real (kind=double) , public, target : : & 

g_wwz ( l_mode : lh_mode , l_mode : h_mode ) = . 
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real (kind=double) , public, target : : & 

g_wwww (0 : 4) = (/ (0 . , i = l, 5) /) , & 

g_wwz z ( l_mode : lh_mode, l_mode : lh_mode) = 0., & 
g_wwza ( l_mode : lh_mode, l_mode : h_mode) = 0. 

(the coupling between four is an exception and is enumerated by the number 
of W at the vertex). All other couplings are fixed by electromagnetic or SU(3)c* 
gauge invariance. For convenience when calculating the couplings, the wavef unc- 
tions are organized in a similar way 

real (kind=double) : : & 

wf unct_w ( l_mode : h_mode, lat_0 : lat_l ) = 0., & 
wf unct_z ( l_mode : h_mode, lat_0 : lat_2 ) = 0. 
real (kind=double) : : & 

wf unct_lep_l ( l_mode : h_mode, gen_0:gen_2, & 

iso_up : iso_down, lat_0 : lat_l ) = 0.,& 
wf unct_lep_r ( l_mode : h_mode, gen_0:gen_2, & 

iso_up : iso_down, lat_l : lat_2 ) = 0.,& 
wf unct_quark_l (l_mode : h_mode, gen_0:gen_2, & 

iso_up : iso_down, lat_0 : lat_l ) = 0.,& 
wf unct_quark_r (l_mode : h_mode, gen_0:gen_2, & 

iso_up : iso_down, lat_l : lat_2 ) = 0. 

The parameters of the Three-Site Lagrangian are defined as 

real (kind=double) , target :: sigma_vev=0 . , g0=0 . , gl=0., & 
g2 = . , x=0 . , lambda=0 . , t=0 . , eps_l, & 

eps_r ( f type_l : f type_q, gen_0:gen_2, iso_up : iso_down) = 0., e 

where t is redundant and defined for convenience only. sigma_vev must be di- 
vided by \/2 to obtain the symmetry breaking scale v defined in chapter 13.31 and 
lambda has to be multiplied by 2 to obtain A (due to a difference in convention). 

The Standard Model parameters are defined as variables which have to be set 
prior to calling the diagonalization routines 

real (kind=double) , public : : & 

me_pdg=0 ._double, mmu_pdg=0 . 1 6_double, mtau_pdg=l . 7 8_double, & 
muq_pdg=0 ._double, mdq_pdg=0 ._double, mcq_pdg=l . 2 5_double, & 
msq_pdg=0 . 95_double, mtq_pdg=17 4 ._double, mbq_pdg=4 . 2_double, & 
mw_pdg=80 . 403_double, mz_pdg=91 . 188_double, & 
e_pdg=0 . 31332 9_double, ww_pdg=2 . 048_double, & 
wz_pdg=2 . 4 4 3_double, wt_pdg=l . 52 3_double 

In addition, there is an error reporting system built into the module which traps 
invalid function calls or points in parameter space at which the analytic formulae 
contain complex roots. Each function pushes its name via a call to estack_push on 
a stack when it begins execution and removes it again via estack_pop at the end. A 
subprogram panic is called if an error occurs, which then prints an error message 
and the contents of the error stack to help localizing the origin of the failure. A 
version of the square root msqrt is provided which traps complex roots. 
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Apart from the Standard Model parameters and the mass / width arrays, most 
variables, subprograms and functions are declared as private and prefixed with 
threeshl_ for export (the variables are exported by pointers). 

Prior to using the library, the subprogram in it must be called 

subroutine init 

# ifdef NLOW 

call nlow_init 

# endif 

call set_names 
call init_pointers 
end subroutine init 

which initializes a couple of internal quantities (if the O (a s ) corrections have been 
enabled at compile time, the corresponding library is also initialized). Afterwards, 
the Standard Model variables can be redefined if desired and the Model is then 
initialized by a call to 

subroutine pdg_init_wgap_bmass (mhw, bmass, el) 

real (kind=double) , intent (in) :: mhw, bmass 

real (kind=double) , optional, intent (in) : : el 

character (len=slength) , parameter :: f name="pdbJnit_wgap_bmass" 

call errstack_push ( f name) 

if ((mhw < mw_pdg) .or. (bmass < 0.)) & 

call panic (err_invalid_parameters, 0) 
call gauge_cpl_f rom_sm_wgap (mw_pdg, mz_pdg, e_pdg, mhw) 
if (present (el) ) then 

eps_l = el 
else 

eps_l = eps_l_of_x (x) 
end if 

lambda = bmass / sqrt (2 ._double) / sigma_vev 

call translate_f ermion_masses ( (/me_pdg, mmu_pdg, mtau_pdg/), & 
(/muq_pdg, mcq_pdg, mtq_pdg/), (/mdq_pdg, msq_pdg, mbq_pdg/)) 
call diagonalize 
call errstack_pop 
end subroutine pdg_init_wgap_bmass 

The call to gauge_cpl_f rom_sm_wgap first calculates the gauge couplings from the 
W, W and Z masses and the electromagnetic gauge coupling e. If ei is not given, 
ideal derealization is assumed and is calculated from x. A is then set from M\> u xk 
and translate_f ermion_masses is called to calculate the er from the SM fermion 
masses. Finally, diagonalize is called. 

For more information on the guts of diagonalize refer to the woven source. In 
a nutshell, this subprogram performs the following steps 

• Calculate the gauge boson wavef unctions and masses. 

• Calculate the fermion wavef unctions and masses. 



• Calculate the couplings. 
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In particular, the masses of the Standard Model particles are not taken from the 
input parameters, but are also calculated from the mass matrices in order to provide 
a simple consistency check. 

An error during the calculation usually triggers the termination of the program. 
This behavior can be modified by setting the flag threeshl_quit_on_panic to 
. false ., in which case the library returns to the caller, printing an error message 
to STDOUT and setting the flag threeshl_error. 

After the couplings have been calculated, the user can call the subprogram 
calculate_widths in order to calculate the widths and fill the width array (only 
the widths of the KK particles are calculated, the t, W and Z widths are taken from 
the Standard Model parameters, and the other Standard Model particles are set to 
zero width). If activated at compile time, the O (a s ) corrections to the heavy quark 
widths can be included. This can be controlled via the global flag threeshl_use_nlow 
which is set to .true . per default. In the calculation, a s is evaluated at the heavy 
quark mass. 

It is also possible to initialize the model in the massless limit described in chap- 
ter !5.3l To this end, the subprogram init_ward must be called 

subroutine init_ward (mx, ct, ph) 

real (kind=double) , intent (in) :: mx, ct, ph 

integer : : gen, iso 

character (len=slength) , parameter :: fname = "iniLward" 

call errstack_push ( fname) 

if ( (ct .le. -l._double) .or. & 

(ct . ge . l._double) .or. (mx == 0.)) & 
call panic (err_invalid_parameters, 0) 
mass_array = 0. 

t = msqrt (1 ._double/ct**2 - l._double) 
x = mx 
e = e_pdg 

gO = e * msqrt ( l._double + x**2 + 1 ._double/t**2 ) 

gl = gO / x 

g2 = gO * t 

sigma_vev = . 

lambda = . 

eps_l = . 

eps_r = . 

wf unct_w ( l_mode, :) = (/cos(ph), sin(ph)/) 

wf unct_w (h_mode, :) = (/-sin(ph), cos(ph)/) 

wf unct_z ( l_mode, :) = ( /-g2/2 ._double/gl - gl/g2 , & 

g2/2 ._double/g0 - gO/2 ._double/g2 , gl/gO + gO/2 ._double/gl /) 
wf unct_z ( l_mode, :) = wf unct_z ( l_mode, :) / & 
msqrt (wfunct_z ( l_mode, lat_0)**2 + & 
wf unct_z ( l_mode, lat_l)**2 + & 
wf unct_z ( l_mode, lat_2)**2) 
wfunct_z (h_mode, : ) = & 

( /-gO/2 ._double, gl, -g2/2 ._double/) / & 
msqrt (g0**2/4 ._double + gl**2 + g2**2/4 ._double) 
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do gen = gen_0, gen_2 ; do iso = iso_up, iso_down 

wfunct_lep_l (l_mode, gen, iso, :) = (/cos(ph), sin(ph)/) 
wf unct_lep_l (h_mode, gen, iso, :) = (/-sin(ph), cos(ph)/) 
wfunct_lep_r (l_mode, gen, iso, :) = (/sin(ph), cos(ph)/) 
wf unct_lep_r (h_mode, gen , iso, :) = (/cos(ph), -sin(ph)/) 
wf unct_quark_l ( l_mode, gen, iso, :) = (/cos(ph), sin(ph)/) 
wf unct_quark_l (h_mode, gen, iso, :) = (/-sin(ph), cos(ph)/) 
wf unct_quark_r ( l_mode, gen, iso, :) = (/sin(ph), cos(ph)/) 
wf unct_quark_r (h_mode, gen ,iso, :) = ( /cos (ph) , -sin (ph) / ) 

end do; end do 

call calculate_couplings 

call errstack_pop 
end subroutine init_ward 

This first sets all masses to zero and calculates the wavefunctions and parameters 
according to chapter 15.31 Then, the couplings are calculated via the overlap of the 
wavefunctions just like in the physical case. 

After using the library, the function finalize should be called. 

O (a s ) corrections to the fermion widths: module nlowidth 

The module nlowidth implements the numerical calculation of the analytical ex- 
pressions given in app. IC.2I As this might be useful in other contexts than the 
Three-Site Model, this module is kept separate from the noweb source of the rest 
of the package and doesn't contain any dependencies on the other packages. The 
analytical expressions are included from external files which have been autogener- 
ated from Mathematica, and a FORTRAN 90 wrapper around the Loop Tools library 
(which is FORTRAN 77 code) is supplied. As there is no woven documentation for 
this module, it is completely reprinted and documented here 

module nlowidth 
use ltglue 

The Loop Tools library outputs status and error messages to STDOUT. Unfortu- 
nately there is no way to change this behavior which disturbs the operation of pro- 
grams like spektrumorthreeshl. To this end, the module calls libc functions to 
redirect these messages to STDERR (a trick shamelessly stolen from the Mathemat- 
ica wrapper around Loop Tools). This relies on FORTRAN 2003 features and can be 
deactivated at compile time via the preprocessor variable no silencer . 

#ifndef NO_S I LENCER 

use, intrinsic : : iso_c_binding 
use, intrinsic : : iso_f ortran_env 
#endif 

implicit none 
private 

All public qualifiers are prefixed with nlow_, everything else is private. 
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public :: nlow_dz, nlow_nlo, nlow_rrad, nlow_lo, nlow_width, & 

nlow_wrel, nlow_init, nlow_f inalize, nlow_alfas, nlow_alf as_mz, & 
nlow_mt, nlow_mb, nlow_mc, nlow_mz, nlow_b2mode, & 
nlow_b2mode_series, nlow_b2mode_full, nlow_b2mode_auto, & 
nlow_b2auto_thresh 

The parameter fa used in fl67| for the definition of the three-particle phasespace 
integrals turns out to be numerically unstable in the massless limit for the final 
state fermion. Therefore, the module can either use the full analytic expression for 
fa or an expansion in 

_ m m 2 

^ 2 2 

ttiq — m\ 

to order O (e 19 ) (using the naming conventions from app.[C]). This choice can be 
controlled via nlow_b2mode with the default setting nlow_b2mode_auto using the 
full expression for fa above nlow_b2auto_thresh and the expansion below it. The 
remaining public parameters which can be modified are the c, b, t and Z masses and 
the value of a s at the Z pole which are required for the evolution of the coupling. 

integer, parameter :: nlow_b2mode_series=0 , nlow_b2mode_f ull=l , & 
nlow_b2mode_auto=2 

integer, parameter : : double=selected_real_kind& 
(precision ( 1 . ) + 1, range (1.) + 1) 

real (kind=double) , parameter :: pi=3 . 1415926535897932385_double 

real (kind=double) :: nlow_alfas_mz=0 . 1176_double, & 
nlow_mt=17 4 ._double, nlow_mb=4 . 2 , nlow_mc=l . 2 5 , & 
nlow_mz=91 . 188_double, nlow_b2auto_thresh=0 . l_double 

integer : : stdout_copy, nlow_b2mode=nlow_b2mode_auto 

The following functions are overloaded and can be called with or without 1x1%. If rri2 
is committed, then the limit rri2 — > is performed (see the actual implementation 
below for details); otherwise, the treatment of rri2 is determined from nl ow_b2mode. 

nlow_rrad calculates the 1—^3 decay width an external gluon, nlow_width 
the full O (a s ) width (including virtual + real corrections) and nlow_wrel returns 
the O (a s ) width normalized to the leading order result. 

interface nlow_rrad 

module procedure rrad_nobs 

module procedure rrad_bs 
end interface nlow rrad 



interface nlow_width 

module procedure width_nobs 
module procedure width_bs 

end interface nlow_width 

interface nlow_wrel 

module procedure wrel_nobs 
module procedure wrel_bs 

end interface nlow_wrel 
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This is the C interface necessary for redirecting the Loop Tools output. 

#ifndef NO_SILENCER 

interface 

function dup (s) result (t) bind (c, name="dup") 

import : : c_int 

integer (kind=c_int ) , value :: s 
integer (kind=c_int) : : t 
end function dup 

function dup2 (s, t) result (r) bind (c, name="dup2") 

import : : c_int 

integer (kind=c_int ) , value : : s, t 

integer (kind=c_int) : : r 
end function dup2 
end interface 
#endif 

contains 

redirect and end_redirect are responsible for the actual redirection. 

subroutine redirect 
integer : : dummy 

# ifndef NO_S I LENCER 

flush (output_unit, iostat=dummy ) 
flush (error_unit, iostat=dummy ) 
dummy = dup2 (2, 1) 

# endif 

end subroutine redirect 

subroutine end_redirect 
integer : : dummy 

# ifndef NO_S I LENCER 

flush (output_unit, iostat=dummy ) 
flush (error_unit, iostat=dummy ) 
dummy = dup2 ( stdout_copy, 1) 

# endif 

end subroutine end_redirect 

The following two functions should be called before and after using the library. 

subroutine nlow_init 

# ifndef NO_S I LENCER 

stdout_copy = dup (1) 

# endif 

call redirect 
call ffini 
call end_redirect 
end subroutine nlow_init 

subroutine nlow_f inalize 
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call redirect 
call ffexi 
call end_redirect 
end subroutine nlow_f inalize 

Calculates the normalization factor for 1 — > 2. mt, mb and raw are mo, m,2 and m\ 
(adhering to the conventions of app.O. 

function wdnorm (mt, mb, raw) result (res) 
real (kind=double) , intent (in) :: mt, mb, mw 
real (kind=double) : : res 

call redirect 

if ( (min (mt, mb, mw) < 0.) .or. (mt < mb + mw) ) & 
print *, "WARNING: u calculatingjDroken u width!" 

res = sqrt((mt**2 - (mb + mw) **2) * (mt**2 - (mb - mw)**2)) / & 

mt**3 / 32 ._double / pi 
call end_redirect 
end function wdnorm 

This function calculates the O (a s ) contribution to the residual of the quark prop- 
agator (|C.10|l . The actual code for AZ is included from dz . inc. m is the fermion 
mass, and 1 the infinitesimal gluon mass to 

function nlow_dz (m, 1) result (res) 
real (kind=double) , intent (in) : : m, 1 
real (kind=double) : : res 
complex ( kind=double ) :: dz 

call redirect 

call setlambda (1**2) 

# include "dz.inc" 

if (aimag (dz) /= O._double) print *, & 
"WARNINGiddzjioL.strictlyj'ealL", aimag (dz) 
res = real (dz) 
call end_redirect 
end function nlow_dz 

Calculates the O (a s ) 1PI correction to the 1 — > 2 width. Again, the actual code is 
included from an extra file, gv and ga are the vector and axial couplings and 1 is 
the gluon mass oj. 

function nlow_nlo (mt, mb, mw, gv, ga, 1) result (res) 
real (kind=double) , intent (in) :: mt, mb, mw, gv, ga, 1 
real (kind=double) : : res 
complex ( kind=double ) :: nlo 

call redirect 

call setlambda (1**2) 

# include "nlo.inc" 

if (aimag (nlo) /= O._double) print *, & 

"WARNING: u nlojioL J strictlyj'ealL", aimag (nlo) 
res = real (nlo) * wdnorm (mt, mb, mw) 
call end_redirect 
end function nlow_nlo 
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The next two private functions are the mi ^ resp. mi = branches of the 
public nlow_rrad and calculate the 1 — > 3 decay width for real gluon radiation. 

In the first case, the code for the analytical expressions of the phasespace in- 
tegrals (IC.llb is taken from II67II and included from external files. Depending on 
nlow_b2mode and mi, the series expansion of fa (see above) is used. The analyt- 
ical expression for the 1 — > 3 in terms of the phase space integrals is read in from 
rrad . inc. 

In the second case, an infrared singularity arises in the limit mi — > which 
cancels between virtual and real corrections. The function sets 7712 = 1CT 50 GeV in 
order to regularize this divergence and then evaluates a version of the 2—^3 width 
where only the logarithms in mi are retained. 

function rrad_nobs (mt, mb, mw, gv, ga, 1) result (res) 
real (kind=double) , intent (in) :: mt, mb, mw, gv, ga, 1 
real (kind=double) :: res, ep 

complex ( kind=double ) :: rrad, k, bO, bl, b2, into, intlOl, intlO, & 
intll, int200, int201, int211, intllO, b2_series 
call redirect 

# include "k.inc" 

# include "bO.inc" 

# include "bl.inc" 

ep = mb*mt/ (mt**2 - mw**2) 
select case (nlow_b2mode) 
case (nlow_b2mode_full) 

# include "b2.inc" 

case (nlow_b2mode_series) 

# include "b2_series.inc" 
b2 = b2_series 

case (nlow_b2mode_auto) 

if (ep < nlow_b2auto_thresh) then 

# include "b2_series.inc" 
b2 = b2_series 

else 

# include "b2.inc" 
end if 

case default 

print *, "WARNING:Jnvalid u nlow_b2mode; u using u nlow_b2mode_full!" 

# include "b2.inc" 
end select 



# 


include 


"intO.inc" 


# 


include 


"int101.inc" 


# 


include 


"intlO.inc" 


# 


include 


"intH.inc" 


# 


include 


"int200.inc" 


# 


include 


"int201 .inc" 


# 


include 


"int21 1 .inc" 


# 


include 


"intHO.inc" 


# 


include 


"rrad. inc" 



if (aimag (rrad) /= O._double) print *, & 
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"WARNING: jracLnot^strictlyj-ealL", aimag (rrad) 
res = real (rrad) 
call end_redirect 
end function rrad_nobs 

function rrad_bs (mmt, mmw, gv, ga, 1) result (res) 
real (kind=double) , intent (in) :: mmt, mmw, gv, ga, 1 
real (kind=double) :: res 

complex (kind=double) :: rrad_bsoft, mb, mt, mw 
mt = mmt ; mw = mmw 
mb = (1 .E-50_double, O._double) 
call redirect 

# include "rracLbsoft.inc" 

if (aimag (rrad_bsoft) /= O._double) print *, & 

"WARNING: jradJDSofLnot^strictlyj'ealL", aimag (rrad_bsoft) 
res = real (rrad_bsoft) 
call end_redirect 
end function rrad_bs 

nlow_lo calculates the leading order 1—^2 width. 

function nlow_lo (mt, mb, mw, gv, ga) result (res) 
real (kind=double) , intent (in) :: mt, mb, mw, gv, ga 
real (kind=double) :: res, lo 

# include "lo.inc" 

res = lo * wdnorm (mt, mb, mw) 
end function nlow_lo 

The following two functions are drivers which assemble the NLO width from the 
leading order result plus the real and virtual corrections. In the m 2 = case, m 2 is 
set to the infrared regulator m 2 = 10~ 50 GeV in order to regularize the divergence 
which cancels between virtual and real corrections, alf as is the value for a s to be 
used for the calculation (see below for a function which evolves a s to the desired 
scale). 

function width_bs (mt, mw, gv, ga, alfas, 1) result (res) 
real (kind=double) , intent (in) :: mt, mw, gv, ga, 1, alfas 
real (kind=double) :: res, low, mb 
mb = 1 . E-5 0_double 

low = nlow_lo(mt, mb, mw, gv, ga) 
res = low + alfas * & 

( (nlow_dz (mt, 1) + nlow_dz (mb, 1)) * low + & 

nlow_nlo (mt, mb, mw, gv, ga, 1) + & 

nlow_rrad (mt, mw, gv, ga, 1)) 
end function width_bs 

function width_nobs (mt, mb, mw, gv, ga, alfas, 1) result (res) 
real (kind=double) , intent (in) :: mt, mb, mw, gv, ga, 1, alfas 
real (kind=double) :: res, low 

low = nlow_lo(mt, mb, mw, gv, ga) 

res = low + alfas * & 
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( (nlow_dz (mt, 1) + nlow_dz (mb, 1)) * low + & 
nlow_nlo (mt, mb, mw, gv, ga, 1) + & 
nlow_rrad (mt, mb, mw, gv, ga, 1)) 
end function width_nobs 

These two function return the NLO with normalized to the LO result. 

function wrel_nobs (mt, mb, mw, gv, ga, alfas, 1) & 
result (res) 

real (kind=double) , intent (in) :: mt, mb, mw, gv, ga, 1, alfas 
real (kind=double) : : res 

res = nlow_width (mt, mb, mw, gv, ga, alfas, 1) & 
/ nlow_lo (mt, mb, mw, gv, ga) 
end function wrel_nobs 

function wrel_bs (mt, mw, gv, ga, alfas, 1) result (res) 
real (kind=double) , intent (in) :: mt, mw, gv, ga, 1, alfas 
real (kind=double) : : res 

res = nlow_width (mt, mw, gv, ga, alfas, 1) & 
/ nlow_lo (mt, 1 . E-5 0_double, mw, gv, ga) 
end function wrel_bs 

nlow_alfas evolves a s to a given scale. The mass thresholds and the value a s at 
the Z pole are global parameters which can be adjusted (see above). 

function nlow_alfas (1) result (res) 
real (kind=double) , intent (in) : : 1 
real (kind=double) :: res, a 

call redirect 

a = nlow_alf as_mz 

if (1 < 0) then 

print *, "WARNING^not^evolving^alfasJo^negative^scale!" 
elseif (1 > nlow_mt) then 

a = evolve (nlow_mt, nlow_mz, a, 5) 

a = evolve (1, nlow_mt, a, 6) 
elseif (1 > nlow_mb) then 

a = evolve (1, nlow_mz, a, 5) 
elseif ( 1 > nlow_mc) then 

a = evolve (nlow_mb, nlow_mz, a ,5) 

a = evolve (1, nlow_mb, a, 4) 
else 

a = evolve (nlow_mb, nlow_mz, a, 5) 
a = evolve (nlow_mc, nlow_mb, a, 4) 
a = evolve (1, nlow_mc, a, 3) 
end if 

res = a 

call end_redirect 
contains 



function evolve (11, 10, aO, n) result (resl) 
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real (kind=double) , intent (in) :: 11, 10, aO 
integer, intent (in) :: n 
real (kind=double) :: resl 
resl = l._double / & 

( (ll._double - 2 ._double*real (n, kind=double) /3 ._double) & 
/ 2._double /pi * log (11 / 10) + l._double / aO) 
end function evolve 
end function nlow_alfas 

end module nlowidth 



Interfacing O'Mega: module tglue 

The tglue module provides the interface to the matrix element code generated by 
O'Mega. This is necessary for several reasons: 

• Several couplings to the photon and to the gluons are not defined inthreeshl 
but are required for the matrix element code to function. 

• The couplings between three gauge bosons must be defined with an addi- 
tional factor i for usage in O'Mega. 

• The W^ff type couplings calculated in threeshl are taken w.r.t. to the Wy 2 
states and need to be divided by an additional factor of \[2 in order to obtain 
the couplings to the W ±<K '\ 

• The couplings between gauge bosons and fermions are decomposed w.r.t. to 
the left- and right-handed projectors, while the O'Mega implementation uses 
VA-type couplings. 

To this end, additional couplings are defined in this module 

complex (kind=double) , public :: & 

g_a_lep, g_a_quark ( iso_up : iso_down) , g_aaww, & 
ig_aww, ig_wwz ( l_mode : lh_mode, l_mode : h_mode) 
complex (kind=double) , public :: & 

g_w_lep_va (1:2, l_mode : h_mode, l_mode : h_mode, & 

gen_0:gen_2, l_mode : h_mode, gen_0 : gen_2 ) , & 
g_w_quark_va (1:2, l_mode : h_mode, l_mode : h_mode, & 

gen_0:gen_2, l_mode : h_mode, gen_0 : gen_2 ) , & 
g_z_lep_va (1:2, l_mode : h_mode, l_mode : lh_mode, & 

gen_0:gen_2, iso_up : iso_down) , & 
g_z_quark_va (1:2, l_mode : h_mode, l_mode : lh_mode, & 

gen_0:gen_2, iso_up : iso_down) 
real (kind=double) , public :: g_s=l . 2 1 8_double 

complex (kind=double) , public :: ig_s_norm, g_s_norm, g_s_norm2 

In order to initialize these couplings, the functions tglue_init (and 
tglue_init_ward for the massless limit) are provided which boast the same 



CHAPTER D. CODE 



153 



interfaces as threeshl_pdg_init_wgap_bmass and threeshl_ward (see above). 
The tglue initialization automatically initializes the threeshl module (and the 
tscript module, see below). Similarly, a tglue_f inalize is provided which calls 
all other initialization routines. 

An example of how to use tscript: program threeshleval 

Instead of discussing the tscript module in detail, the source treeshleval is 
reprinted, a small program which takes Mb u ik and optionally (for nonideal 
derealization) as well as a tscript function definition as command line arguments 
and returns the evaluated function. 

program threeshleval 

use threeshl 

use tdefs 

use tscript 

use tglue 

implicit none 

real (kind=double) :: mhw=-l., mbulk=-l., eps_l 
real (kind=double) , pointer : : value 
logical :: el_present = .false, 
character (len=slength) :: fun 

The tscript module defines the tscript_tokenize_ob ject data type which is 
initialized by a call to tscript_create_tob ject with the character string which 
is to be interpreted as a function definition. 

type (tscript_tokenize_ob ject) :: tobject 

First, the program checks if it has been correctly called. If this is not the case, it 
prints the syntax and quits (see below). Otherwise, the variables are set from the 
command line arguments. 

if (.not. ( (iargc () >= 3) .and. (iargc () <= 5))) call print_usage 
call getarg (1, fun); mhw = str_to_double (fun) 
call getarg (2, fun); mbulk = str_to_double (fun) 
call getarg (3, fun) 
if (iargc () == 4) then 

if (trim (fun) == "no_nlow") then 
threeshl_use_nlow = .false. 

else 

eps_l = str_to_double (fun) 
el_present = .true, 
end if 

call getarg (4, fun) 
end if 

if (iargc () ==5) then 

eps_l = str_to_double (fun) 

call getarg (4, fun) 

if (trim (fun) == "no_nlow") then 
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three shl_use_n low = .false, 
else 

call print_usage 
end if 

call getarg (5, fun) 
end If 

In order to parse the function string, a tscript_tokenize_ob ject must be cre- 
ated. This is then decoded by a call to tscript_decode_f spec in order to obtain a 
pointer to the corresponding quantity. If an error occurred during the parsing, the 
syntax is printed (can be switched off via the flag tscript_show_syntax) and the 
program is terminated if threesh_quit_on_panic is set. 

tobject = tscript_create_tob ject ( fun) 
value => tscript_decode_f spec (tobject) 

The model is initialized by a call to tglue_init and the function value is printed 
to STDOUT. After calling tglue_f inalize, the program quits. 

if (el_present) then 

call tglue_init (mhw, mbulk, el=eps_l) 
else 

call tglue_init (mhw, mbulk) 
end if 

print '(F25.15)', value 

call tglue_f inalize 

contains 

This subprogram prints the syntax and then in turn calls tscript_print_syntax 
which shows the tscript syntax and quits if threeshl_quit_on_panic is set. 

subroutine print_usage 

print * , "usage: JhreeshlevaLmhwjmbulk4epsJ]4no_nlow] Junction" 
print * 

print * , "mhw^and jribulk^are Jhe^mass^of Jhe Jneavy JA/^and Jhe Jxjlkjnass;" 
print *, "if^epsJJs^not^present.JtJsjdetermined^viaJdeal" 
print *, "derealization. J\sJorJheJunction^specifier:" 
print * 

call tscript_print_syntax 
end subroutine print_usage 

This tries to interpret a character string as a floating point number. 

function str_to_double (str) result (res) 
character (len=slength) , intent (in) :: str 
real (kind=double) : : res 
integer : : errstat 

read (unit=str, fmt='(F50.50)', iostat=errstat) , res 

if (errstat .ne. 0) call print_usage 
end function str_to_double 

end program threeshleval 
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D.2 O'Mega Module 

In this section, only a part of the O'Caml source of the model module is presented. 
The reader interested in studying the full code is referred to the commented ocaml- 
web source. 

Options 

In order to allow for an easy setup of the various options discussed !5.3[ the O'Mega 
module is implemented as a functor which maps two option modules to the final 
model module. The signatures of these modules are 

module type Threeshl_ options = 
sig 

val include-ckm : bool 
val include-hf : bool 
val diet : bool 
end 

type qcd -implementation = Disabled \ Colflow 

module type Threeshl- colopt = 
sig 

val o : qcd -implementation 

end 

The meaning of the different options is: 

• include-ckm: Include flavor violating couplings. 

• include-hf: Include the heavy fermions. 

• diet: Setting this to true discards all couplings between the W' and the lep- 
tons or the first two quark generations (only implemented in the case of 
include -ckm= false). 

• Theeshl_ colopt. o: Method for the treatment of color. Setting this to Colflow in- 
cludes the gluons and the QCD couplings with the sign conventions choosen 
correctly for the colorizer module to work, Disabled excludes the gluons. 

The functor is then defined as 

module Threeshl' (Module-options : Threeshl -options) 
(Module- colopt : Threeshl- colopt) = 

Other options that the user might want to control can be set via command line 
options when the binary is invoked 

let default -width = ref Timelike 
let all-feynman = ref false 
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let options = Options . create [ 

"constant_width", Arg. Unit (fun _ — > default. width := Constant), 

"use^constant^width^ (also^in^t-channel) "; 
"custom_width", Arg. String (fun x — > default-width := Custom x), 

"use^custom^width"; 
"cancel-widths", Arg. Unit (fun _ — > default-width := Vanishing), 

"use^vanishing^width"; 
"all_f eynman", Arg. Unit (fun _ — > all-feynman := true), 

" assign^f eynman^gauge^propagators^to^all^gauge^bosonsXn" 

" "\t ( f or^checking^the^ward^identities ) ; " 

' "use^only^if ^you^^really^^knowXn" 

" " Xtwhat^you^are^doing"] 



Flavors 

The first duty of the module is the definition of a type flavor which enumerates the 
different particles. This is done by first defining quantum numbers that encode KK 
mode, generation, "charge" (differentiates between particles / antiparticles and is 
also defined for neutrinos) and isospin and then defining flavor through suitable 
constructors. 

type kkmode = Light \ Heavy 

type generation = GenO \ Genl \ Gen2 

type csign = Pos \ Neg 

type isospin = Iso-up \ Iso-down 

type fermion = 

| Lepton of (kkmode x csign x generation x isospin) 
| Quark of (kkmode x csign x generation x isospin) 

type boson = 

W of (kkmode x csign) 
| Z of kkmode 

I G 

type flavor = Fermion of fermion \ Boson of boson 

Two functions revmap and revmap2 are defined which apply a list of functions 
to a value or resp. each element of a list of values and return the result as a flat list. 

let revmap funs v = List. map (fun x — > x v) funs 

let revmap2 funs vals = ThoList.flatmap (revmap funs) vals 

Together with a couple of functions that map to the constructors 
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let lepton kk cs gen iso = Lepton (kk, cs, gen, iso) 

let quark kk cs gen iso = Quark (kk, cs, gen, iso) 

let w kk cs = W (kk, cs) 

let z kk = Z kk 

let flavor _of _/ x = Fermion x 

let flavor .of -b x = Boson x 

and several functions which loop a list of functions over the quantum numbers 

let loop-kk flist = revmap2 flist [Light; Heavy] 

let loop-cs flist = revmap2 flist [Pos; Neg] 

let loop -gen flist = revmap2 flist [GenO; Genl; Gen2] 

let loop -iso flist = revmap2 flist [Iso -up; Iso-down] 

let cloop-kk flist = match Module _options. include -hf ^ with 

| true ->■ loop-kk flist 

| false —7- revmap flist Light 

these two functions allow for the easy creation of a list of all particles which is 
required by O'Mega 

let all-leptons = loop-iso (loop_gen (loop-cs (cloop-kk [lepton] ))) 
let all-quarks = loop_iso( loop-gen (loop-cs (cloop-kk [quark] ))) 
let all-bosons = (loop-cs (loop_kk [w] )) @ [Z Light; Z Heavy] @ [A] @ 
(match Module. colopt.o\N\th Colflow -)• [G] \ _ -)• []) 

\e\ flavors () = (List. map flavor _of _/ (all_leptons @ all-quarks)) @ 
(List. map flavor _of _b all-bosons) 

(using clopp-kk instead of loop-kk automatically implements the option include -hf 
for the exclusion of the heavy fermions). 

Helpers 

Apart the enumeration of the particles, the module is required to provide some 
query functions which O'Mega can use in order to obtain information on the prop- 
erties of the particles. The enumeration of the particles by quantum number allows 
to do this in a rather convenient way as demonstrated by the examples of SU(3)c 
and Lorentz representation. 

let color = 

let quarkrep = function 

| (_, Pos, _, _) -> Color. SUN 3 

| (_, Neg, _, _) -> Color. SUN (-3) 
in function 

| Fermion (Quark x) — > quarkrep x 
| Boson G -> Color. AdjSUN 3 
| _ — > Color. Singlet 
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let lorentz = 

let spinor = function 

| (_, Pos, _, _) — > Spinor 

| (_, Neg, _, _) — > ConjSpinor 
in function 

| Fermion (Lepton x) \ Fermion (Quark x) —> spinor x 
| Boson (W _) | Boson (Z _) — > Massive-Vector 
| Boson A —7- Vector 
| Boson G —7- Vector 

Another type of helper which is very important for the interaction with WHIZARD 
is the function pdg which returns the PDG code assigned to a particle. 

let int. of -csign = function Pos — > 1 | Neg — 1 

let int_of_gen = function GenO -)■ 1 | Genl -)■ 2 | Gen^ -)■ 3 

let pdg = 

let iso-delta = function /so_ down — >■ | Iso-up — >■ 1 

in let gen-delta = function GenO — > | Genl — > 2 | Gen2 — > 4 

in let kk_delta = function Zig/it — > | i^ea?;y — > 9900 

in function 

| Fermion ( Lepton (kk, cs, gen, iso)) — > 

(int - of - csign cs) x (11 + (gen -delta gen) + (iso -delta iso) + (kk-delta kk)) 
| Fermion ( Quark (kk, cs, gen, iso)) —> 

(int-of -csign cs) x (1 + (gen_delta gen) + (iso-delta iso) + (kk-delta kk)) 
| Boson (W (kk, cs)) — > (int _ of _ csign cs) x (24 + (kk-delta kk)) 
| Boson (Z kk) — > 23 + (kk-delta kk) 
| .Boson A -)• 22 
| .Boson G ^ 21 

This function is designed such that it adheres to the Monte Carlo numbering scheme 
for the Standard Model particles and assigns the PDG code of their partner prefixed 
with 99 to their KK partners, e.g. 11/ - 11 to the e~/e+ and 9911/ - 9911 to the 

e /e + . 

Although not included in the interface, a very important function is the trans- 
lation of a flavor into the constant identifying it in the FORTRAN module (see app. 

ED. 

let bcdi-of -flavor = 
let prefix = function 

I Fermion (Lepton (Heavy, _, _, _)) | Fermion (Quark (Heavy, _, _, _)) 

| Boson (W (Heavy , _)) | Boson (Z Heavy) — > "h" 

y II II 

in let rump = function 

| Fermion (Lepton spec) — > (match spec with 
| (_, _, GenO, Iso-up) — > "nue" 
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| (_, _, GenO, Iso_down) — > "e" 

| (_, _, Genl , Iso-up) — > "numu" 

| (_, _, Genl, Iso-down) — > "mu" 

| (_, _, Gen2, Iso-up) — > "nutau" 

| (_, _, Gen2, Iso-down) — > "tau") 
| Fermion (Quark spec) — > (match spec with 

| (_, _, GenO, Iso_up) ->■ "u" 

| (_, _, GenO, Iso-down) — > "d" 

| (_, _, Genl, Iso_up) — > "c" 

| (_, _, Genl, Iso-down) — > "s" 

| (_, _, Gen2, Iso.up) -)• "t" 

| (_, _, Gen2, Iso-down) ->■ "b") 
| Boson (W.) "w" | Boson (Z J) "z" 

| Boson A —7- invalid_arg "Csmodelsl . bcd_of .flavor : ^no^bcd^f or^photon ! " 
| Boson G — > invalid. arg "Csmodelsl . bcd_of _f lavor : ^no^bcd^f or^gluon ! " 
in function x — > (prefix x) ~ (rump x) ~ "_bcd" 

Other functions which are required by O'Mega for example to translate particle 
flavors in character strings suitable for the backend or to communicate with the 
user on the command line are constructed in a likewise fashion. 

Couplings and vertices 

In order to be able to represent the coupling constants as closely to the conventions 
of the FORTRAN module as possible, another type referring to the combination of 
KK modes at a vertex rather than the individual modes is defined. 

type kk2 = Light2 \ Heavy 2 \ Light. Heavy 

let loop-kk2 flist = revmap2 flist [Light2; Heavy2; Light _Heavy] 

let cloop-kk2 flist = match Module .options, include -hf ^ with 

| true ->■ loop_kk2 flist 

I false ->• revmap flist Light2 

The coupling constants are then defined to mimic their FORTRAN counterparts as 
closely as possible. 

type constant = 

| G_aJ,ep \ G-a_quark of isospin 

| G-aww | G-aaww 

| G-W-lep of (kkmode x kkmode x generation x kkmode x generation) 

| G -w -quark of (kkmode x kkmode x generation x kkmode x generation) 

| G-Z_lep of (kkmode x kk2 x generation x isospin) 

| G-Z-quark of (kkmode x kk2 x generation x isospin) 

| G-wwz of (kk2 x kkmode) 

| G-wwzz Of (kk2 x kk2) 

| G-iuwza of (kk2 x kkmode) 
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| G-wwww of int 
| Gs 
| IGs 

| Gs2 

Using the loop-xx function allows for a compact definition of the vertex lists, 
e.g. for the 7// and Zff type vertices 

let vertices -all = 
let vgen kk gen = 

((Fermion (Lepton (kk, Neg, gen, Iso-down)), Boson A, Fermion (Lepton (kk, Pos, gen, 
Iso-down))), FBF(\, Psibar, V, Psi), G-a_lep) 
in loop -gen (cloop-kk [vgen]) 

let vertices -aqq = 
let vgen kk gen iso = 

((Fermion (Quark (kk, Neg, gen, iso)), Boson A, Fermion (Quark (kk, Pos, gen, 
iso))), FBF(1, Psibar, V, Psi), G-a-quark iso) 
in loop-iso (loop_gen (cloop-kk [vgen])) 

let vertices-zll = 

let vgen kkz kk-f kk-fbar gen iso = 

((Fermion (Lepton (kk-fbar, Neg, gen, iso)), Boson (Z kkz), 
Fermion (Lepton (kk_f, Pos, gen, iso))), 
FBF (1, Psibar, VA2, Psi), 

G-Z-lep (kkz, get-kk2 (kk-f, kk-fbar), gen, iso)) 
in loop-iso (loop_gen (cloop-kk (cloop-kk (loop-kk [vgen] )))) 

let vertices-zqq = 

let vgen kkz kk-f kk-fbar gen iso = 

((Fermion (Quark (kk-fbar, Neg, gen, iso)), Boson (Z kkz), 

Fermion (Quark (kk-f, Pos, gen, iso))), 

FBF (1, Psibar, VA2, Psi), 

G-Z-quark (kkz, get-kk2 (kk-f , kk-fbar), gen, iso)) 
in loop-iso (loop-gen (cloop-kk (cloop-kk (loop-kk [vgen] )))) 

The Wff vertex lists are defined in different versions in order to implement the 
different options , e.g. for Wqq type vertices 

let vertices - wqq _no -ckm = 

let vgen kkw kk-f kk-fbar iso-f gen = 

((Fermion (Quark (kk-fbar, Neg, gen, conj-iso iso_f)), 

Boson (W (kkw, (match iso-f with Iso-up -)• Neg \ _ ->■ Pos))), 

Fermion (Quark (kk-f, Pos, gen, iso-f))), 

FBF (1, Psibar, VA2, Psi), 

G-W -quark (kkw, (match iso_f with Iso_ up — > kk_f \ _ — > kk-fbar), gen, 
(match iso-f with Iso- up ->■ kk-fbar \ _ ->■ kk-f), gen)) 
in loop-gen (loop-iso (cloop-kk (cloop-kk (loop-kk [vgen] )))) 
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let vertices _wqq_no-ckm -diet = 
let filter = function 

I ((Fermion (Quark (Light, _, gen, _)), Boson (W (Heavy, _)), 
Fermion (Quark (Light, _, _, _))), _, _) — > 
(match gen with Gen2 -)• true I _ ->■ false) 
| _ ->• true 
in List. filter filter vertices _wqq-no-ckm 

let vertices _wqq = 

let wgen fcfcw gen-f kk-fbar gen-fbar iso_f = 

((Fermion (Quark (kk-fbar, Neg, gen-fbar, conj-iso iso-f)), 

Boson (W (kkw , (match iso-f with Iso_up ->■ iVeg | _ ->■ Pos))), 
Fermion (Quark (kk-f, Pos, gen-f , iso-f))), 
FBF (1, Psibar, VA2, Psi), 
G-W- quark (match iso-f with 

| Iso_up — > (kkw, kk-f, gen_f, kk_fbar, gen_fbar) 
| Iso-down — > (kkw, kk-fbar, gen-fbar, kk-f, gen-f))) 
in loop -iso (loop_gen (cloop-kk (loop_gen (cloop-kk (loop_kk [vgen] ))))) 

The actual vertex list passed to O'Mega is then assembled according to the selected 
options. 

let vertices () = (vertices -all @ vertices-aqq @ 
(match Module -options, diet with 
| false — > vertices-wll 
| true ->■ vertices -wll -diet) @ 
(match (Module -options. include -ckm, M odule_ options. diet) with 
| (true, false) — > vertices-wqq 
| (false, false) — > vertices _wqq_no_ckm 
| (false, true) — > vertices -wqq-no-ckm -diet 
| (true, true) ->■ raise (Failure 

("Modules4 . Threeshl . vertices : ^CKM^matrix^together^with" " 
"^option^diet^is^not^implemented^yet ! "))) @ 
vertices -zll @ vertices -zqq @ vertices -aww @ vertices -zww @ 
(match Module -colopt.o with 

| Colflow — > vertices-gqq @ vertices-ggg 

I - -> []), 

vertices _aaww @ vertices _wwzz @ vertices -wwza @ vertices -wwww @ 
(match Module-colopt.o with 
| Colflow —7- vertices -gggg 

!-->[]) 
,[]) 
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D.3 WHIZARD Quirks 

In WHIZARD, models are described by a model file and a piece of FORTRAN glue. 
The model file defines free and derived parameters, the particles and a vertex list. 
The FORTRAN glue is called prior to the evaluation of the first matrix element and 
has to take care of setting up and providing all parameters and couplings that are 
required by the matrix element generation code to function. In principle, writing 
the model file is a straightforward task, and the FORTRAN glue would be the per- 
fect place to perform the tglue initialization. However, in the case of the Three-Site 
Model, implementing the masses of the particles in this fashion turns out to be a 
problem. 

In the particle definition, WHIZARD requires the mass to be either zero or a 
free / derived parameter. In the Three-Site Model however, the masses are com- 
plicated functions which are calculated by the threeshl module and stored in the 
mass_array array. It is possible to assign the members of this array to derived pa- 
rameters like 

derived mhz mass.array (hz_bcd) 

However, the members of this array get only initialized upon calling the tglue ini- 
tialization which, if this would be called from the FORTRAN glue, would be only 
after the derived parameters have been evaluated. This way, the correct masses 
would not propagate into the phasespace code which therefore would fail. 

The only way to get around this without modifying WHIZARD would be to 
define the masses as free parameters and calculate them with an external program 
which writes out a WHIZARD input file. However, this would be rather cum- 
bersome and error-prone. Therefore, the version of WHIZARD which contains 
the Three-Site Model has been modified to generate code which calls tglue_init 
before the derived parameters are calculated. The amplitude is modified not to 
use the usual FORTRAN glue (which is unnecessary for the Three-Site Model and 
therefore omitted) but instead use threeshl and tglue. 

Other minor changes include the build system which is modified to properly 
use the external library and an external O'Mega tree (WHIZARD includes its own 
version of O'Mega). However, these are minor issues which could be avoided, 
while the problem discussed in the last paragraph can only be solved cleanly by a 
change to the WHIZARD infrastructure. 

As a result of these issues, the Three-Site Model is currently not included in the 
official distribution of the 1.9x branch of WHIZARD but only available in a mod- 
ified WHIZARD package which can be downloaded from the URL quoted at the 
beginning of this chapter. The upcoming WHIZARD 2.0, however, will remove 
these limitations and will contain the implementation as part of the official distri- 
bution. 
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